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1.  INIRUUUC110N 


THE  PROBLEM 

Pavement  surfaces  for  runways,  taxiways,  end  aprons  of  modern  airports 
play  an  important  part  in  the  overall  performance  of  the  facility.  Pavement 
design,  construction,  and  maintenance  are  important  in  providing  acceptable 
surfaces  at  reasonable  overall  cost.  Experience  has  shown  that  pavements  on 
expansive  soils  often  require  rehabilitation  at  an  early  stage  of  pavement 
life. 

ExvonoLvo  soils  are  clayey  soils  which  exhibit  significant  volume  changes 
as  a  result  of  soil  moisture  variations.  Problems  with  airport  pavements  are 
caused  by  this  kind  of  volume  change  in  the  underlying  soil.  These  changes 
produce  unacceptable  differential  movements  in  the  pavement  surface,  causing 
surface  roughness  and  cracking.  An  example  of  a  pavement  surface  distorted  by 
expansive  soil  behavior  is  shown  in  Figure  1-1.  The  pavement  shown  is  capable 
of  supporting  the  load  it  was  designed  to  carry.  However,  as  an  aircraft 
travels  across  the  surface,  the  roughness  excites  the  airplane,  inducing 
unacceptable  vertical  accelerations  in  the  aircraft.  When  these  accelerations 
reach  a  certain  level,  pilots  are  unable  to  safely  negotiate  the  aircraft  on 
the  pavement,  passengers  are  uncomfortable,  and  damage  to  the  aircraft  or 
cargo  may  occur. 

The  behavior  of  expansive  soils  influences  the  performance  of  pavement 
structures  founded  on  these  soils.  Therefore,  understanding  the  properties 
that  make  a  soil  expansive  or  not  expansive  aids  in  evaluating  specific  soils 
for  pavement  design.  The  objectives  of  pavement  design  are  to  incorporate 
methods  in  the  construction  process  that  prevent  or  reduce  the  effects  of  soil 
volume  changes. 

OBJECTIVES  AND  SCOPE 

The  Systems  Research  and  Development  Service  of  the  Federal  Aviation 
Administration  initiated  a  study  of  expansive  soils  in  1975.  In  Phase  I  of 
the  study,  the  New  Mexico  Engineering  Research  Institute  (NMERI)  used  the 
technical  literature  to  update  airport  pavement  design  procedures.  The  pro¬ 
posed  design  method  which  resulted  from  the  literature  review  is  shown  in 
Figure  1-2.  NMERI  began  a  program  to  develop  this  method  in  February  1977. 
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FIGURE  1-1.  PAVEMENT  SURFACE  DISTORTED  BY  SOIL  BEHAVIOR. 


INDICATOR 

TEST 


FIGURE  1-2.  RECOMMENDED  DESIGN  PROCEDURE 


Phase  II  was  a  laboratory  study  of  soil  characterization  methods.  Phase 
III  was  a  field  study  of  airport  pavement  and  in  situ  soil  behavior  at  three 
sites  with  distinctly  different  climates.  Phase  IV  consisted  of  analysis  and 
design  recommendations.  The  results  of  Phases  I  and  II  are  reported  in  McKeen 
(1976)  and  McKeen  and  Nielsen  (1978).  Phase  III  and  Phase  IV  are  the  subject 
of  this  report. 


McKeen,  R.  G. ,  (1976)  Design  and  Construction  of  Airport  Pavements  on  Expan¬ 
sive  Soils ,  FAA-RD-76-66,  Federal  Aviation  Administration,  Washington, 
DC. 

McKeen,  R.  G.  and  Nielsen,  J.  P.,  (1978)  Characterizing  Expansive  Soils  for 
Airport  Pavement  Design ,  FAA-RD-78-59 ,  Federal  Aviation  Administration, 
Washington,  DC. 
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II.  BACKGROUND:  STATE  OF  THE  ART 


During  the  past  twenty  years  significant  research  efforts  around  the 
world  have  been  devoted  to  the  study  of  expansive  soil  interaction  with  struc¬ 
tures.  Costly  structures  and  equipment  are  damaged  by  expansive  soil  behavior, 
and  the  standard  methods  in  soil  mechanics  are  not  capable  of  describing  the 
phenomena  adequately  for  proper  design.  The  results  of  much  of  this  research 
have  been  available  in  the  technical  literature  for  several  years.  It  remains 
to  apply  these  evolving  theories  to  practical  problems  and  present  design 
methods  that  are  implementable. 

SUCTION 

Expansive  soils  cause  problems  beneath  structures  because  they  change 
volume.  The  soils  may  swell  with  great  force,  lifting  a  structure,  or  shrink 
away  from  it,  removing  support  altogether.  Only  two  causes  of  these  volume 
changes  are  involved  in  most  problems:  changes  in  load  and  changes  in  mois¬ 
ture  condition.  The  response  of  soils  to  loads  has  been  a  part  of  soil 
mechanics  since  its  beginning.  In  contrast,  the  volume  changes  accompanying 
variation  in  moisture  condition  have  only  recently  been  studied.  Soil  volume 
changes  caused  by  moisture  variation  are  the  predominant  factor  involved  in 
pavement  design  on  expansive  soils. 

The  potential  of  a  soil  to  undergo  volume  change  because  of  moisture 
variation  is  determined  by  the  type  of  clay  contained  in  the  soil  and  the 
amount  present.  Soils  that  do  not  contain  clay  (particles  less  than  0.002  pm 
in  size)  do  not  exhibit  expansive  soil  behavior.  In  addition,  the  type  of 
clay  minerals  in  the  soil  is  important.  In  engineering  work,  three  groups  of 
clay  minerals  are  identified.  In  descending  order  of  volumetric  activity, 
these  are  the  smectite,  illite,  and  kaolinite  groups.  These  minerals  are  com¬ 
posed  of  layers  of  silica  tetrahedra  and  gibbsite  stacked  in  various  config¬ 
urations  (Figure  2-1).  Naturally  occurring  soils  are  invariably  composed  of 
mixtures  of  these  groups  plus  other,  less  common,  varieties. 

Volume  changes  of  these  minerals  result  from  adsorbed  water  layers 
attached  to  the  clay  mineral  surfaces.  Water  molecules  are  dipolar  and  are 
attracted  to  the  surfaces  by  unbalanced  negative  charges.  The  imbalance  is 
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created  by  isomorphous  substitution  of  metallic  ions  in  the  clay  mineral 
lattice  systems.  Large  surface  areas  and  charge  densities  require  large  quan¬ 
tities  of  water  to  satisfy  the  imbalance.  As  large  amounts  of  water  are 
adsorbed,  the  clay  platelets  are  forced  apart,  occupying  a  larger  volume. 

Soil  moisture  suction  is  defined  as  a  macroscopic  property  of  the  soil 
which  indicates  the  intensity  with  which  a  soil  will  attract  water.  Suction 
is  normally  defined  as  a  negative  gage  pressure  (energy  i  volume  =  units  of 
force  per  unit  area),  but  it  should  not  be  confused  with  pore  water  pressure. 
Pore  water  pressure  is  normally  associated  with  the  density  of  a  liquid, 
distance  from  a  free  water  surface,  and  surface  tension  forces.  Suction 
results  from  (a)  the  interplay  of  attraction  and  repulsion  forces  of  charged 
clay  particles  and  polar  water  molecules,  (b)  surface  tension  forces  of  water, 
(c)  solution  potentials  caused  by  dissolved  ions,  and  (d)  gravity  potential. 

The  representation  of  suction,  the  sum  of  all  these  forces,  as  an  equivalent 
height  of  water  has  been  called  the  capillary  model.  This  model  was  contro¬ 
versial  until  1960  when  substantial  agreement  was  reached  on  its  validity 
(Aitchison,  1961 ) . 

In  engineering  problems,  soil  moisture  suction  is  considered  to  be  com¬ 
posed  of  matrix  suction  and  osmotic  (or  solute)  suction.  Matrix  action  is 
the  negative  pressure  which  will  hold  soil  water  in  equilibrium  through  a 
porous  membrane  with  the  same  soil  water  within  a  sample  of  the  soil.  Matrix 
suction  results  from  surface  adsorption  and  capillary  forces.  Osmotic  or 
solute  suction  is  a  negative  gage  pressure  which  will  hold  pure  water  in  equi¬ 
librium  with  soil  water  through  a  membrane  which  allows  only  water  molecules 
to  pass.  Osmotic  suction  results  from  variation  in  ion  concentration  in  the 
pore  fluid.  Matrix  suction  is  directly  influenced  by  mechanical  loads  and 
water  content  whereas  osmotic  suction  is  not. 

Two  independent  stress  variables  have  been  used  to  describe  the  state  of 
stress  in  unsaturated  soils.  The  preferable  stress  state  variables  are 
(o  -  U  )  and  (U  -  U)  where  o  =  total  stress,  U  -  pore  air  pressure,  and 
U  =  P°re  water  pressure.  The  term  (a  -  U  )  is  called  the  total  stress  tom, 

W  d 

and  (Ua  -  Uw)  is  called  the  matrix  suction  term.  This  combination  of  stress 

Aitchison,  G.  D.,  (1961)  "Relationships  of  Moisture  Stress  and  Effective  Stress 

Functions  in  Unsaturated  Soils,"  in  Pom  Pressure  and  Suction  in  Soils, 

Butterworths ,  London. 
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state  variables  is  most  satisfactory  because  the  effects  of  environmental 
variables  can  be  readily  separated  in  terms  of  stress  changes.  This  approach 
assumes  that  1)  is  approximately  atmospheric  and  that  the  osmotic  component  of 

a 

suction  remains  constant.  These  assumptions  are  adequate  for  many  engineering 
problems . 

The  logarithmic  unit  pF  is  used  to  measure  suction  in  the  present  paper, 
along  with  units  of  pressure  in  kPa.  pF  is  defined  as 

pF  =  1 og  i  o  (pressure  in  cm  of  water)  (2-1) 

Figure  2-2  illustrates  the  suction  phenomenon.  The  vertical  axis  repre¬ 
sents  suction  (free  energy)  of  the  water.  Also  on  the  vertical  axis  is  a 
scale  representing  molecular  layers  of  adsorbed  water,  taken  from  Michurin  and 
Lytayev  (1967).  An  explicit  relationship  between  water  layers  and  suction  is 
indicated.  Data  for  several  materials  were  used  to  illustrate  that  all  porous 
hygroscopic  materials  exhibit  similar  adsorption  of  water.  Curves  for  Portland 
cement  concrete  (McKeen,  1969),  lumber  (ASCE,  1975),  filter  paper  (McQueen  and 
Miller,  1977),  clay  minerals,  and  real  soils  are  illustrated.  For  practical 
use,  these  relationships  may  be  represented  as  linear.  Note  the  influence  of 
the  surface  areas  indicated  in  Figure  2-2  for  the  clays  and  soils. 

The  range  shown  in  Figure  2-2  is  the  adsorbed  water  range  (2.5  to  6.25  pF). 
Below  2.5  pF,  additional  large  quantities  of  water  are  held  in  the  capillary 
range.  Above  6.25  pF,  water  molecules  intimately  associated  with  crystalline 
structures  of  minerals  are  mobilized.  The  volume  changes  of  interest  in  study¬ 
ing  pavement  behavior  on  expansive  soils  occur  in  the  range  of  adsorbed  water 
from  2.5  to  5.5  pF  (McKeen  and  Nielsen,  1978). 


ASCE, 


(1  975)  lv,  v, i  ttni-  ■  tui’ cv;. 

Division,  1975. 


fury  ,  Structural 


McKeen,  R.  G.,  (1969)  "Shrinkage-Cracking  Characteristics  of  Structural  Light¬ 
weight  Concrete,"  Thesis,  Texas  A  &  M  University,  College  Station,  Texas. 

McQueen,  I.  S.,  and  Miller,  R.  F.,  (1977)  "The  Wide-Range  (Filter  Paper)  Method 
for  Measuring  the  Forces  of  Moisture  Retention,"  personal  correspondence. 

Michurin,  B.  N.,  and  Lytayev,  I.  A.,  (1967)  "Relationship  between  Moisture  Con¬ 
tent,  Moisture  Tension,  and  Specific  Surface  Area  in  Soil,"  ;.>v,  .'  .7 

’M;*t ’ ,  Vol  .  8,  pp.  1093-1103. 
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FIGURE  2-2.  WATER  CONTENT-FREE  ENERGY  RELATIONSHIPS 


SOIL  RESPONSE  TO  SUCTION  CHANGE 


The  response  of  soil  volume  to  suction  changes  may  be  studied  in  much 
the  same  manner  as  soil  consolidation  in  response  to  applied  load  changes 
(Aitchison  and  VJoodburn,  1969;  McKeen  and  Nielsen,  1978;  Escario,  1969;  Lytton, 
1977;  Compton,  1970;  Fredlund  and  Morgenstern,  1977).  However,  some  basic 
differences  must  be  recognized.  Figure  2-3  shows  consolidation  data  for  tests 
on  a  sample  of  expansive  soil.  Curve  A  is  the  result  of  a  conventional  one¬ 
dimensional  consolidation  test.  Curve  B  is  the  result  of  a  hydrostatic  con¬ 
solidation  test.  Curve  C  represents  consolidation  resulting  from  suction 
induced  by  drying  samples  to  various  moisture  conditions.  In  all  cases  the 
stresses  are  effective  stresses.  The  effective  stresses  for  the  suction  curve 
were  computed  using  an  equation  presented  by  Aitchison  (1961): 

1  p" 

o"  =  T50  (S  *  p,,)  +  Z  °-3p"  ‘  AS  (2-2) 

o 

o"  =  effective  stress  component  caused  by  p" 

S  =  degree  of  saturation 

AS  =  100  -  S 

p"  =  pressure  deficiency,  kg/cm2 

An  important  difference  between  suction  (Curve  C)  and  the  conventional 
test  (Curve  A)  is  that  suction  changes  are  not  one-dimensional.  The  suction 
forces  on  a  soil  element  are  body  forces  rather  than  externally  applied  forces 


Aitchison,  G.  D.,  and  Woodburn,  J.  A.,  (1969)  "Soil  Suction  in  Foundation 

Design,"  Proceedings,  7th  International  Conference  on  Soil  Mechanics  and 
Foundation  Engineering,  Mexico  City,  Vol .  21,  pp.1-8. 

Compton,  Phil  V.,  (1970)  A  Study  of  the  Swelling  Behavior  of  an  Expansive 

Clay  as  Influence  l  by  the  Clay  Micros  true tura t  Soil  Suction,  dnl  Exterruil 
Loading,  AFWL-TR-/0-26,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force 
Base,  New  Mexico. 

Escario,  Ventura,  (1969)  "Swelling  of  Soils  in  Contact  with  Water  at  a  Negative 
Pressure,"  Proceedings,  2nd  International  Research  and  Engineering  Confer¬ 
ence  on  Expansive  Clay  Soils,  College  Station,  Texas,  pp.  207-217. 

Fredlund,  D.  G.,  and  Morgenstern,  N.  R.,  (1977)  "Stress  State  Variables  for 
Unsaturated  Soils,"  Journal  of  the  Geotechnical  Enainceri  n  ■■  ’division , 

ASCE,  Vol.  103,  No.  GT5,  pp.  447-466. 

Lytton,  R.  L.,  (1977)  "The  Characterization  of  Expansive  Soils  in  Engineering," 
presentation  at  the  Symposium  on  Water  Movement  and  Equilibrium  in 
Swelling  Soils,  American  Geophysical  Union,  San  Francisco,  California, 
(unpubl i shed) 


CURVE  A  (LOADING) 


CURVE  A  (UNLOADING) 


at  the  surface  of  the  element.  Suction  forces  yield  a  single  curve  for  a 
soil,  regardless  of  whether  it  is  wetting  or  drying.  When  the  effective 
stress  is  reduced,  the  relieved  body  forces  remold  the  soil.  In  the  case  of 
applied  loads  at  the  surface,  this  remolding  does  not  occur;  therefore  the 
rebound  curves  from  mechanical  loading  are  different  from  the  loading  curves. 

Because  suction  stresses  are  largely  effective  stresses  applied  directly 
to  the  soil  skeleton  as  body  forces,  volume  change  studies  may  be  made  using 
total  suction  values.  In  practical  engineering  problems,  it  is  the  slope  of 
the  volume  change  versus  suction  curve  that  quantifies  soil  response  to  mois¬ 
ture  changes.  Thus  the  suction  compression  index,  y^,  is  defined  as  the  slope 
of  the  volume-total  suction  curve  (Lytton,  1977).  No  constraint  on  the  manner 
in  which  the  volume  change  is  determined  is  specified  here.  The  importance  of 
this  aspect  is  discussed  later. 


Yh 


AV/V. 

log  10  ^f 

h. 


(2-3) 


where 

Yh  =  suction  compression  index 
AV/V.j  =  change  of  volume  with  respect  to  initial  volume 
hp  h^  =  final  and  initial  total  suction,  arithmetic  units 

Values  of  y^  may  be  determined  by  several  testing  methods.  A  measured 
volume  change  and  suction  change  are  required,  and  these  data  must  cover  the 
range  of  moisture  suction  applicable  to  the  problem  under  study.  The  volume 
measurement  may  be  made  in  a  conventional  oedometer  which  constrains  the 
sample  to  one-dimensional  volume  change.  Another  test,  the  COLE  (coefficient 
of  linear  extensibility)  test,  measures  unrestrained  three-dimensional  volume 
change  (McKeen  and  Nielsen,  1978).  Natural  soils  will  experience  lateral  res¬ 
traint  conditions  that  vary  between  the  two  conditions  represented  in  these 
tests.  A  dry,  cracked  clay  soil  profile  will  approach  the  COLE  condition;  an 
uncracked  clay  will  approach  the  fully  restrained  condition.  The  selection  of 
a  proper  y^  value  requires  some  knowledge  of  the  lateral  restraint  experienced 
by  the  soil . 
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Another  factor  to  consider  in  arriving  at  values  is  the  effect  of 
applied  loads  in  reducing  the  soil  swell.  The  compression  index  for  load 
changes  is  determined  as  follows  (Lytton,  1977): 


AV/Vi 

Y°  log  ff 
°1 

where 

Yo  =  total  stress  compression  index 
AV/V^  =  change  of  volume  with  respect  to  initial  volume 
Of,  a.j  =  final  and  initial  total  stress 


(2-4) 


The  two  compression  indexes,  Yh  and  ya*  may  be  combined  into  an  incre¬ 
mental  volume  strain  equation  for  the  soil  (Lytton  et  al,  1976;  McKeen  and 
Nielsen,  1978): 


AV 

T  “ 


-Y 


Ah 

h  hj 


-Y  — 

'a 


(2-5) 


Evaluation  of  these  expressions  reveals  the  surface,  describing  soil 
behavior  in  response  to  load  and  suction  changes.  This  surface  is  illustrated 
in  Figure  2-4. 


THE  GILGAI  PHENOMENON 

A  natural  clay  soil  left  to  interact  with  climatic  forces  will  develop 
an  irregular  surface  microrelief  pattern.  As  illustrated  in  Figure  2-5  (from 
Spotts,  1974),  the  pattern  develops  in  seven  stages.  First,  sedimentary  clay 
is  deposited.  Following  deposition,  drying  occurs,  producing  surface  shrinkage 
cracks.  Further  desiccation  results  in  the  formation  of  major  cracks  penetrat¬ 
ing  far  into  the  clay.  The  extent  of  penetration  and  spacing  are  determined 


Lytton,  R.  L. ,  Boggess,  R.  L. ,  and  Spotts,  J.  W.,  (1976)  "Characteristics  of 
Expansive  Clay  Roughness  of  Pavements,"  swelling  Soils,  Transportation 
Research  Record  No.  568,  Washington,  DC,  pp.  9-23. 

Spotts,  J.  W.,  (1974)  "The  Role  of  Water  in  Gilgai  Formation,"  Dissertation, 
Texas  ASM  University,  College  Station,  Texas. 


SUCTION-VOLUME  SURFACE  FOR  SWELLING  SOILS. 
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by  a  number  of  factors  of  which  the  degree  and  rate  of  drying  and  the  clay 
type  are  the  most  significant.  Subsequent  development  of  microrelief  occurs 
as  the  clay  is  wetted  and  dried  by  cyclic  variation  of  the  climate.  Because 
water  penetrates  deeper  at  the  major  cracks,  greater  swell  occurs.  As  the 
next  cycle  of  drying  proceeds,  the  higher  areas  dry  first,  causing  the  major 
cracks  to  reappear  in  approximately  the  same  locations.  Thus,  over  a  period 
of  several  years,  an  undulating  surface  pattern  called  yiTyi:  develops  on  the 
soil  surface.  Characteristics  of  some  gilgai  surfaces  found  in  the  United 
States  are  shown  in  Table  2-1. 

TABLE  2-1.  GILGAI  CHARACTERISTICS 


Soil  -  Location 

Dominant  wavelenqths, 
m  (ft) 

Maximum 
amp! itude, 
m  (ft) 

Burl eson-Texas 

4.5  (15)  to  9.0  (30) 

0.12  (0.4) 

Houston  Black-Texas 

4.0  (13)  to  9.0  (30) 

0.16  (0.5) 

Puerco-New  Mexico 

1.5  (5) 

0.60  (2.0) 

Oenhawken-Texas 

4.5  (15)  to  9.0  (30) 

0.46  (1.5) 

The  gilgai  structure  extends  through  the  depth  of  soil  that  interacts  with 
the  climate,  the  active  zone  of  the  soil  profile.  Gilgai  structural  develop¬ 
ment  begins  when  the  soil  is  deposited.  The  depth  of  actual  gilgai  occurrence 
depends  on  the  most  extreme  climate  in  the  soil  history  and  may  not  represent 
current  active  zone  depth.  Also,  because  of  successive  deposition  cycles, 
some  gilgai  may  be  buried  beneath  more  recent  deposits  (O'Neill  and  Ghazzaly, 
1977).  Gilgai  structure  in  soil  profiles  causes  soil  response  to  environmental 
change  to  vary  depending  on  position  within  the  structure.  The  variation  is 
produced  by  soil  leaching,  deposition  of  dissolved  salts  by  evaporating  water, 
development  of  internal  drainage  paths,  and  infiltration  of  surface  material 
via  cracks.  Gilgai  structure  is  most  evident  where  soil  surfaces  are  left 
undisturbed  for  long  periods  of  time  (Figure  2-6).  Thus,  an  expansive  soil 
profile  will  respond  differentially  from  point  to  point  when  subjected  to  an 
environmental  change. 


O'Neill,  Michael  W.,  and  Ghazzaly,  Osman  I.,  (1977)  "Swell  Potential  Related 
to  Building  Performance,"  Joumix!  <>f  oot.  •hnL\:'  h!min,.o»:n : 
cion ,  ASCE,  Vcl  .  103,  No.  GT  12,  pp'.  1363-1379. 
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EFFECTS  OF  PAVEMENT  ON  SOILS 

When  an  airport  pavement  is  constructed  on  an  expansive  clay  soil,  the 
environment  to  which  the  soil  is  exposed  changes  dramatically.  The  pavement 
structure  serves  as  a  barrier  to  moisture  exchange  with  the  atmosphere.  This 
barrier  usually  causes  the  active  zone  soil  to  exchange  moisture  with  the  soil 
beneath  it  until  an  equilibrium  condition  is  approached.  The  ability  of  mois¬ 
ture  to  move  through  a  soil  profile  via  fissures,  seams,  and  so  on,  may  inter¬ 
rupt  the  equilibration  process  or  prevent  it  entirely.  If  an  area  contains 
numerous  routes  for  water  movement  in  the  soil,  a  large  area  may  heave,  alter¬ 
ing  the  pavement  profile  by  displacing  sections  of  pavement.  Examples  of  this 
alteration  are  given  in  the  elevation  profiles  in  Appendix  A.  Once  the  soil 
attains  an  equilibrium  condition  beneath  the  pavement,  the  edges  will  still  be 
subjected  to  seasonal  variations. 

The  initial  equilibration  of  the  soil  may  produce  differential  vertical 
movement  in  the  pavement  profile.  Additional  disturbance  of  the  equilibrated 
soil  will  occur  if  cracks,  joints,  or  other  pavement  surface  characteristics 
permit  infiltration  of  surface  water  into  the  subgrade.  Usually  this  Kind 
of  distress  is  localized  about  a  specific  moisture  source.  Both  initial 
equi 1 ibration  and  subsequent  infiltration  involve  suction  changes,  and  both 
result  in  unacceptable  profile  characteristics.  These  two  mechanisms  are  the 
major  cause  of  difficulties  in  the  maintenance  of  airport  pavements  on  expan¬ 
sive  soils. 

A  significant  factor  contributing  to  pavement  profile  deterioration  is 
the  pavement  response  to  aircraft  loads.  The  loads  produce  permanent  deforma¬ 
tion  in  the  subgrade  which  accumulates  with  repetition.  In  addition,  any 
differential  support  beneath  the  pavement  structure  caused  by  differential 
movement  in  the  soil  will  be  enhanced  by  the  aircraft  loads.  Thus  the  two 
processes  of  pavement  deterioration  complement  one  another. 

SUMMARY 

The  problem  at  hand  is  to  develop  design  methods  to  prevent  the  rapid 
development  of  excessively  rough  surfaces  on  airport  pavements  built  on  expan¬ 
sive  clays.  The  roughness  exhibited  by  such  pavements  is  the  result  of  normal 
soil  behavior  in  response  to  the  conditions  imposed  by  pavement  construction. 
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Once  the  pavement  imposes  a  change  in  the  environment,  considerable  time  may 
be  required  to  reestablish  moisture  equilibrium  in  the  soil.  The  literature 
in  this  field  indicates  that  soil  characterization  using  suction  is  the  best 
approach  to  understanding  equilibrium. 

Three  mechanisms  may  cause  distress  when  airport  pavements  are  founded  on 
expansive  soils.  First,  normal  soil  equilibration  with  the  boundary  condi¬ 
tions  imposed  by  pavement  construction  may  produce  differential  surface  move¬ 
ment.  Second,  soil  beneath  pavement  edges  swells  and  shrinks  in  response  to 
normal  climate  cycling.  This  behavior  is  strongly  affected  by  shoulder  and 
drainage  design.  Third,  the  infiltration  routes  through  the  pavement  (e.g., 
cracks  and  joints)  result  in  localized  heave  near  the  point  of  infiltration. 
The  first  two  mechanisms  are  considered  design  problems  while  the  third  is 
primarily  a  maintenance  problem. 
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III.  SITES  STUDIED 


Three  sites  were  selected  for  the  Phase  III  field  study  portion  of  this 
investigation:  Gallup  Municipal  Airport,  New  Mexico;  Dallas/Fort  Worth 
Regional  Airport,  Texas;  and  Jackson  Municipal  Airport,  Mississippi.  The 
bases  of  the  selections  were  (1)  previous  expansive  soil  problems  in  the  area, 
(2)  existing  pavement  structures,  and  (3)  variation  in  climate  types  among  the 
three  sites.  The  following  paragraphs  briefly  describe  the  geology,  climate, 
and  soil  properties  of  the  three  sites. 

Geological  information  was  extracted  from  the  literature.  Soil  properties 
were  determined  from  reports  written  during  the  construction  of  pavements  at 
the  airports  and  from  borings  made  at  the  sites.  Specific  references  are 
given  in  the  discussion  of  the  sites.  Borings  were  made  to  a  depth  of  6  m 
(20  ft)  at  each  site.  Most  of  the  data  discussed  here  were  taken  from  one 
borehole  at  each  site.  Samples  taken  from  the  boreholes  were  subjected  to 
standard  engineering  tests  in  the  laboratory.  Properties  investigated  were 
clay  content  (C),  liquid  limit  (LL),  plasticity  index  (PI),  and  water  content 
(w). 

The  climate  for  each  site  is  categorized  using  a  procedure  developed  by 
Thornthwaite  (1948).  The  climate  cycle  is  described  by  balancing  the  rain¬ 
fall  (R),  potential  evapotranspiration  (PE),  and  soil  water  holding  capacity 
(HC).  In  this  method  the  cycle  is  divided  into  periods  as  follows: 

Depletion  (DP)  PE  -  R,  soil  drying 

Deficit  (OF)  PE  R,  soil  moisture  depleted 

Recharge  (RC)  PE  <  R,  soil  wetting 

Surplus  (S)  PE  <  R,  water  drains  away 

The  Thornthwaite  Index  (TI)  is  derived  as  follows: 

TT  100S  -  60DF 
1  PE 

The  TI  provides  an  indication  of  PE  and  R  balance  on  an  annual  basis. 
Negative  values  indicate  dry  climates,  and  positive  values  mean  wetter  climates. 

Thornthwaite,  C.  W.,  (1948)  "An  Approach  Toward  a  Rational  Classification  of 

Climate,"  •  leo.ji\iyhu\il  ,  Vol .  38,  No.  1,  pp.  54-94. 
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GALLUP  AIRPORT  SITE 


The  Giilluj  airport  site  (GAL)  is  located  in  the  city  of  Gallup,  New 
Mexico,  in  the  northwo; tern  part  of  the  state.  The  surface  soil  consists  of 
a  lacustrine  deposit  of  recent  origin.  The  clay  layer  is  estimated  to  be 
approximately  91  m  (300  ft j  hick.  Very  little  information  is  available  on 
the  specific  deposit  on  which  the  airport  is  constructed.  However,  the  mater 
ial  is  derived  by  erosion  from  Manms  Shale  in  the  surrounding  area.  An  in- 
progress  study  of  Mancos-der ived  soils  over  a  four-state  area  has  yielded  some 
information  applicable  to  the  site  (Potter  et  al,  1979): 

As  a  result  of  its  formation  in  a  marine  environment,  the 
Mancos  has  a  high  content  of  sodium  and  sulfate.  The  influence 
of  variable  amounts  of  organic  matter  and  calcium  carbonate 
affect  the  weathering  pattern  of  exposed  shale  and  the  physical 
and  chemical  make-up  of  shale-derived  soils.  Mancos  soils  have 
a  high  clay  content,  generally  including  a  high  percentage  of 
expandable  clays.  The  soils  swell  on  wetting  and  often  exhibit 
poor  water  penetration.  High  concentrations  of  sodium  in  the 
crystal  lattice  exaggerate  this  condition. 

The  X-ray  studies  of  clays  from  shale  samples  show  that 
the  typical  mineralogy  consists  of  small  amounts  of  kao finite 
and  illite,  and  larger  amounts  of  interstratified  illite-smec- 
Lite,  the  smectite  composed  principally  of  montmorrillonite. 

Except  for  a  few  localities  near  igneous  intrusions,  the  inter¬ 
stratified  species  is  approximately  5 Ob  illite  [I]--50'T  smectic 
[S].  The  range  observed  is  601  -  40S  to  401  -  60S,  with 
majority  near  50  -  50. 

The  main  modification  to  the  materials  at  the  Gallup  site  has  been  the 
leaching  associated  with  erosion  and  lacustrine  deposition.  A  boring  log 
is  shown  in  Table  3-1  along  with  information  from  other  investigations. 

Initial  suction  tests  reveal  a  very  uniform  profile,  indicating  that  the  mois¬ 
ture  is  in  equilibrium. 

The  present  runway  and  taxiway  pavements  were  constructed  of  229  mm 

(9  in)  of  compacted  subyrade,  152  mm  (6  in)  of  lime-stabilized  subgrade, 

229  mm  (9  in)  of  cement- treated  base,  prime  coat,  and  51  mm  (2  in)  of  bitumi¬ 
nous  surface  course;  critical  areas  had  a  76-mm  (3-in)  surface  course  and  a 
279-mm  (11-in)  base.  The  cement- treated  base  consists  of  pulvarized  reclaimed 


Potter,  Loren  0.,  Reynolds,  Robert  C.,  and  Louderbough,  Ellen  T.,  (1979) 

/y,'>;v»" .•  .»_<*  in:  '.'i'.-  ■>?  \  vTv.ru  ■'//•  n ,  final 

Report,  U.S.  forest  Service,  Contract  No.  16-729-GR,  University  of  New 
Mexico,  february  1979. 


TABLE  3-1.  GALLUP  BORING  LOGS. 


Construction 
tes  t i ng 


Red -brown  clay 
very  stiff  to 
hard 


NMERI  Boring  No.  1  Laboratory  test  data 


c 

LL 

PI 

28 

49 

23 

81 

95 

58 

75 

51 

82 

85 

87 


53 

72 

75  79  49 


C  =  clay  content 
LL  =  1 iquid  1 imit 
PI  =  plasticity  index 
w  =  water  content 


asphaltic  pavement,  red  dog  (a  local  burned  clay),  old  base  materials,  and 
approximately  4  percent  cement.  The  pavements  were  repaired  in  1972  with  deep 
patching  and  an  overlay  of  114  mm  (4.5  in)  of  dense-graded  bituminous  surface. 
An  open-graded  porous  friction  course  (PFC)  was  then  placed.  Shoulders  were 
constructed  on  both  the  runway  and  taxiway.  The  shoulders  were  3  m  (10  ft) 
wide  and  consisted  of  127  mm  (5  in)  of  crusher  fines  and  a  double  surface 
treatment  with  cover  aggregate. 

Weather  at  the  site  is  dry  with  cold  winters.  The  precipitation  and 
potential  evapotranspi ration  are  shown  in  Figure  3-la.  Note  the  absence  of 
surplus  at  anytime  during  the  year.  A  large  moisture  deficit  is  apparent 
during  the  summer  months.  These  factors  combine  to  produce  very  dry  soil 
conditions.  The  actual  weather  during  the  experiment,  shown  in  Figure  3-lb, 
included  an  extremely  wet  winter  followed  by  normal  weather. 


•JACKSON  AIRPORT  SITE 

The  Jackson  airport  site  (JSN)  is  located  southeast  of  the  city  of 
Jackson,  Mississippi.  Within  the  airport  boundaries,  surface  soils  consist 
of  weathered  Yazoo  clay  and  alluvium.  Boring  logs  for  this  test  program 
(NMERI  Nos.  1  and  2)  and  a  log  of  a  nearby  boring  made  during  construction  of 
the  airport  are  shown  in  Table  3-2. 

Baughman  described  the  Yazoo  clay  as  follows  (Baughman,  1971): 

The  unweathered  Yazoo  clay  is  a  more  or  less  homogeneous 
unit  of  blue-green,  blue-gray,  and  gray  calcareous  clay  with 
occasional  pyrite,  many  thin  soft  fossils  and  fossi 1 i ferous 
zones,  foraminifera ,  occasionally  slightly  sandy  and  usually 
limy  near  the  base.  Some  thin  beds  of  soft  marly  limestone 
are  often  present  in  the  upper  portion  of  the  formation. 

Yazoo  clay  is  greenish  yellow-tan.  On  the  surface  it  is 
usually  blocky  and  crumbly  when  dry  with  much  black  staining 
which  is  due  to  the  presence  of  manganese.  Calcarious  nodules 
on  weathered  surfaces  of  the  clay  are  common  in  some  areas. 

Pyrite  is  not  found  in  weathered  Yazoo.  Gypsum  is  present  in 
the  weathered  zone  in  the  form  of  selenite  crystals.  These 
crystals  are  a  weathering  phenomenon  and  do  not  appear  below 
the  weathered  zone,  which  is  normally  27  to  31  feet. 


Baughman,  Wilbur  T.,  (1971)  "Rankin  County  Geology  and  Mineral  Resources," 
Mississippi  Geological,  Economic,  and  Topographic  Survey, 
Jackson,  Mississippi . 
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FIGURE  3-1.  THORNTHWAITE  DATA  FOR  GALLUP. 


TABLE  3-2.  JACKSON  BORING  LOGS 


Depth , 
m  (ft) 
Or 


1 

(5) 


3 

00) 


05) 

5)- 


6(- 

(20) 


Construction 

testing 


Pavement,  base, 
and  fill 


Stiff  brown  and  gray 
silty  clay  (weathered 
Yazoo) 


Stiff  tan  and  gray 
(weathered  Yazoo) 


NMERI  Borings 
No.  1  and  2 


Gray  silty  clay 
(alluvium  fill ) 


Very  stiff  gray 
and  tan  clay 
(weathered  Yazoo) 


Laboratory  test  data 
C  LL  PI  w 

20  36  14  26 


40  64  42  25 


Tan  and  light  gray 
clay  (weathered  Yazoo) 


55  76  47  28 


Light  gray  and  tan 
clay  with  oxidized 
seams  (weathered 
Yazoo) 


72  114  77  39 


C  =  clay  content 
LL  =  liquid  limit 
PI  -  plasticity  index 
w  =  water  content 
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The  Yazoo  clay  varies  from  approximately  340  to  460  feet 
thick  in  Rankin  County.  The  average  thickness  of  the  Yazoo 
in  Rankin  County  is  approximately  375  feet. 

Further  information  is  provided  from  the  results  of  a  Mississippi 
Highway  Department  study  (Teng  et  al,  1972).  Findings  significant  to  the 
present  work  were  as  follows: 

1.  Mineralogical  analysis:  smectite,  65  percent;  kaolinite,  20  percent; 
mica,  10  percent;  quartz,  5  percent. 

2.  Horizontal  and  vertical  particle  orientations  are  not  significantly 
different. 

3.  Because  of  small  seams  of  sand  in  the  clay,  no  moisture  change 
occurs  as  a  result  of  capillary  rise.  The  only  mechanism  by  which  moisture 
changes  occur  is  surface  runoff. 

4.  The  field  coefficient  of  permeability  is  extremely  low  for  remolded 
Yazoo  clay. 

5.  The  most  severe  pavement  maintenance  problem  is  the  occurrence  of 
fissures  caused  by  desiccation  and  their  subsequent  filling  with  water. 

The  present  pavement  at  the  Jackson  airport  consists  of  varying  amounts 
of  medium-dense  clayey,  fine  sand  fill,  generally  ranging  from  203  to  381  mm 
(8  to  15  in)  thick.  The  fill  is  overlain  by  152  mm  (6  in)  of  cement- treated 
base,  and  229  to  356  mm  (9  to  14  in)  of  portland  cement  concrete  pavement. 

The  east  runway  and  east  parallel  taxiway  have  been  overlaid  with  varying 
amounts  of  asphaltic  concrete,  152  to  203  mm  (6  to  8  in)  thick  in  most  areas. 

A  considerable  amount  of  distress  has  been  produced  in  the  pavement  by 
swelling  soil  behavior. 

Figure  3-2a  illustrates  a  normal  weather  cycle  for  the  Jackson  airport 
site.  The  dry  hot  summers  and  wet  winters  which  normally  occur  in  Jackson  are 
the  factors  producing  a  large  swing  in  soil  moisture  conditions.  However, 
because  of  an  unusually  wet  cycle,  as  shown  in  Figure  3-2b,  the  soil  was  sat¬ 
urated  (surplus  condition)  throughout  the  field  experiment. 

Teng,  T.  C.  Paul,  Mattox,  Robert  M.,  and  Clisby,  M.  Barrett,  (  1972)  rf 

/ha  iih>  Chvja  m  Hd.tt?  }  tv  Ihsaijn ,  Final  Report,  MSHD-RD-72-045, 

Vol .  1  and  Vol .  2,  Mississippi  State  Highway  Department,  Jackson, 
Mississippi . 
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OALLAS/FORT  WORTH  AIRPORT  SITE 


The  Oallas/Fort  Worth  airport  site  (DFU)  is  located  between  the  cities 
of  Dallas  and  Fort  Worth,  Texas.  The  site  is  in  the  Gulf  Coastal  Plain  physio¬ 
graphic  province.  Formations  outcropping  on  the  airport  are  the  Woodbine  and 
Eagleford.  All  work  in  the  present  study  involved  soils  derived  from  the  Upper 
Cretaceous  Eagleford  Formation.  The  Eagleford  Formation  has  been  subdivided 
into  a  number  of  members  at  various  locations.  At  the  study  site,  the  parent 
material  has  weathered  to  a  heavy  black  clay  at  the  surface  [0  to  1.2  m 
(4  ft)].  Below  this  material  is  a  tan  gray,  fissured,  heavy  clay  containing 
numerous  seams  filled  with  sand.  Typical  logs  of  borings  are  shown  in 
Table  3-3.  X-ray  diffraction  testing  indicates  that  the  clays  are  dominantly 
montmorril Ionite  with  smaller  amounts  of  i 11 i te  and  kaolinite.  The  area  is 
widely  known  for  problems  encountered  in  construction  and  performance  of 
engineering  structures  on  these  soils. 

The  present  pavement  at  the  DFW  Regional  Airport  was  constructed  in  the 
early  1 97 0 ‘ s  after  substantial  testing  of  materials  and  test  sections  and  a 
great  amount  of  analysis  and  study.  The  taxiway  section  studied  in  the 
present  work  consisted  of  229  mm  (9  in)  of  lime-stabilized  subgrade,  229  mm 
(9  in)  of  cement-treated  base,  and  381  mm  (15  in)  of  Portland  cement  concrete 
pavement.  The  lime-stabilized  subgrade  was  constructed  using  6  to  7  percent 
lime.  The  subqrade  was  covered  with  about  1.8  L/m2  (0.4  gal/yd:)  of  asphalt 
emulsion  to  form  a  seal.  The  sealed  subgrade  was  then  left  undisturbed  for  at 
least  75  days,  permitting  equilibrium  of  the  soil  moisture  conditions.  As  of 
the  present  report,  the  pavements  have  performed  well  and  exhibit  no  distress. 

The  normal  weather  cycle  for  the  OFW  Regional  Airport  is  shown  in 
Figure  3-3a.  The  cycle  shows  both  deficit  and  surplus  areas,  indicating  that 
the  normal  climatic  cycle  will  produce  large  variations  in  the  soil  moisture 
condition.  Departures  from  the  normal  cycle  will  result  in  prolonged  deficit 
or  surplus  conditions.  Figure  3-3b,  the  actual  weather  during  the  investiga¬ 
tion,  reveals  a  nearly  normal  weather  cycle  for  the  DFW  site  during  the  field 
experiment. 
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TABLE  3-3.  DALLAS/ FORT  WORTH  BORING  LOGS 


Airmail  Facility 
Depth,  Boring  TB3 

m  (ft)  (Nov  71) 


, 

Fill 

1 

1  | 
. 

Dark  brown 

"  | 

i 

clay 

(5) 

; 

Tan  brown 
clay 

Z\ 

1 

. 

' 

t 

3 

■ 

■ 

, 

(10) 

1 

1 

• 

4 

(15) 

- 

1 

Light  tan  and 
gray  clay  with 

r 

i 

sand  seams  and 
bentonite  seams 

5 

-  i 

! 

1 

; 

6 

! 

1 

1 

(20) 

i 

■ 

1 

1 

7 

/ 

Dark  blue  shale 

(30) 

... 

(weathered) 

8 

' 

9 

i 

! 

(40) 

Dark  blue  shale 

NMERI  Boring  No.  4 
(Oct  78) 


C 

LL 

PI 

w 

Heavy  dark 

56 

70 

44 

20 

brown  clay 

55 

76 

50 

28 

67 

73 

46 

25 

Weathered  shale 
tan  and  light  gray 
with  sand  seams 


64  68  43  28 


C  =  clay  content 
LL  =  liquid  limit 
PI  =  plasticity  index 
w  =  water  content 
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IV.  SUCTION  STUDIES 


The  present  state  of  the  art  (Lytton,  1977;  Fredlund  anu  Mortens  tern, 

1977;  Aitchison  and  Woodburn,  1969;  Fredlund  et  ?1,  1980;  Tor  characterizing 
expansive  soils  shows  that  volume  changes  result  directly  from  changes  in 
aDplied  load  and  changes  in  suction.  Soil  response  to  surtion  changes  is  the 
volume  change  of  interest  in  problems  with  pavements  on  expansive  soil  a.  The 
following  paragraphs  describe  the  methods  used  in  studying  suction,  some 
results  of  laboratory  and  field  measurements,  and  the  soil  properties  obtaineo 
with  suction-based  testing. 

i EST  MET  tons 

Thermocouple  P sychrometers 

The  idea  of  inferring  suction  values  from  thermocouple  asyc. geometer  (TcP) 
measurements  of  relative  humidity  was  introduced  by  Spanner  ( i 9bl ) .  Numerous 
researchers  have  worked  toward  the  improvement  oF  equipment  and  procedures 
for  use  of  suction  in  various  research  fields.  Through  the  expenditure  of  a 
great  effort  primarily  at  the  U.S.  Salinity  Laboratory  (Riverside,  California), 
Utah  State  University,  and  the  U.S.  Forest  Service  {Logan,  Utah),  many  vari¬ 
ables  have  been  studied.  The  diversity  of  applications  has  ranged  from  studies 
of  plant  growth  to  irrigation,  forest  hydrology,  soil  mechanical  behavior,  and 
dust  emissions  from  roadways.  Much  can  be  learned  from  the  wealth  of  docu¬ 
mented  studies  involving  these  instruments.  The  intention  here  is  to  famili¬ 
arize  the  reader  with  the  instrument  and  some  of  the  precautions  necessary  for 
proper  usage. 

Principles  involved— A  ;  rywhiowo  tcv  (Figure  4-1)  is  defined  as  two  simi¬ 
lar  thermometers  with  the  bulb  of  one  being  kept  wet  so  that  the  cooling  that 
results  from  evaporation  makes  it  register  a  lower  temperature  than  the  dry 
one;  the  difference  between  the  readings  represents  a  measure  of  the  dryness 
of  the  atmosphere. 

Fredlund,  0.  G.,  Hasan,  J.  U.,  and  Filson,  H.  L.,  (1980)  "The  Prediction  of 
Total  Heave,"  Ppo<?nedin,-;8 ,  4th  International  Conference  on  Expansive 
Soils,  Vol .  1,  ASCE,  pp.  1-17. 

Spanner,  0.  C.,  (1951)  "The  Peltier  Effect  and  Us  Use  in  the  Measurement  of 
Suction  Pressure,"  Joumtl  .-•/  />;••  vrr.’Ktrl  /«<.-,  Vol.  11,  op.  145  168, 
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CONSTANTAN 

REFERENCE  JUNCTION 
(DRY  BULB) 


FIGURE  4-1.  A  PEI  TIER  THERMOCOUPLE  PSYCHROMETER . 


The  difference  between  the  two  readings  is  called  the  wet  b.-.'b  d^'vession. 

For  a  known  wet  bulb  depression  and  dry  bulb  temperature,  the  ratio  of  the 
water  vapor  pressure  in  the  air  to  that  of  pure  free  water  at  the  same  tempera¬ 
ture  and  pressure  can  be  computed.  This  ratio  is  the  relative  humidity. 

Psychrometers  have  been  used  for  many  years  in  studying  the  environment. 
One  essential  feature  of  conventional  instruments  is  that  the  psychrometer  is 
ventilated  in  order  to  yield  valid  measurements.  Traditionally  the  ventilation 
rate  was  a  critical  factor  involved  in  psychrometry.  Spanner  (1951)  demon¬ 
strated  that,  by  reducing  the  size  of  the  wet  bulb,  the  psychrometer  could  be 
used  in  static  air. 

Once  the  ventilation  requirement  is  waived,  a  procedure  for  measuring 
wet  bulb  depression  is  necessary.  Spanner  employed  two  principles  of  thermo¬ 
electricity,  the  Peltier  effect  and  the  Seebeck  effect.  The  Seebeck  effect 
states  that  a  current  will  flow  around  a  loop  of  two  dissimilar  metals  if  the 
two  junctions  are  kept  at  different  temperatures .  The  Peltier  effect  states 
that,  if  a  small  current  is  passed  through  a  thermocouple  sensing  junction 
(wet  bulb)  in  the  proper  direction,  cooling  of  the  junction  occurs;  and  an 
exactly  equal  amount  of  heating  occurs  at  the  reference  junction  (dry  bulb). 
Once  the  temperature  of  the  sensing  junction  drops  below  the  dew  point  of  the 
atmosphere  to  which  it  is  exposed,  water  will  condense  on  the  junction.  The 
maximum  attainable  cooling  is  theoretically  about  1.8°C  (Verbrugye,  19/4).  At 
this  point,  the  temperature  of  the  sensing  junction  has  been  dropped,  and  the 
temperature  of  the  reference  junction  has  been  raised.  In  order  to  use  the 
technique,  Spanner  constructed  massive  heat  sinks  at  the  reference  junction. 

The  design  was  intended  to  reduce  the  temperature  change  at  the  reference  junc¬ 
tion  to  10  3  to  10~’4OC.  In  this  way  the  reference  junction  could  be  assumed 
not  to  change  and  to  be  equal  to  the  temperature  of  the  water  bath  in  which 
the  entire  apparatus  was  immersed.  Thus,  the  dry  bulb  temperature  is  known. 

Wnen  the  cooling  current  is  stopped,  water  condensed  on  the  sensing 
junction  begins  to  evaporate,  maintaining  a  lower  temperature  at  the  sensing 
junction  than  at  the  reference  junction.  The  wet  bulb  depression  maintained  is, 

Verbrugge,  J.  C.,  (1974)  "Contribution  A  La  Mesure  De  La  Succion  Ft  Oe  La 

Pression  Interstitiele  Dans  Les  Sols  Non  Sntures,"  Brussels  Free  Univer¬ 
sity,  School  of  Applied  Sciences,  Jacques  Verdeyen  Laboratory.  (USACCRREL, 
Draft  Translation  No.  5P9). 
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as  in  conventional  psychrometers ,  a  function  of  the  relative  humidity.  As 
long  as  this  temperature  difference  exists,  an  electromotive  force  (emf)  is 
produced  between  the  leads,  in  accordance  with  the  Seebeck  effect.  Thus,  the 
sensor  output  is  a  function  of  wet  bulb  depression  at  the  sensing  junction 
with  respect  to  the  reference  junction.  The  thermocouple  output  resulting 
from  the  Seebeck  effect  reflects  the  temperature  difference  between  the  sens¬ 
ing  and  reference  junctions  regardless  of  the  cause,  e.g.,  evaporative 
cooling,  initial  temperature  difference,  heating  caused  by  the  Peltier  cur¬ 
rent.  Determining  the  cause  of  the  temperature  difference  is  a  major  problem 
in  using  the  sensors  for  studying  water  potential  in  soil  systems. 

The  relationship  between  the  wet  bulb  depression  and  the  water  vapor 
pressure  in  the  air  is  affected  by  a  number  of  factors.  In  the  theoretical 
treatment  presented  by  Rawlins  (1966)  these  include 

1.  Radii  of  wet  junction,  equilibration  chamber,  and  thermocouple  wire 

2.  Thermal  conductivities  of  wire  and  air 

3.  Diffusity  of  water  vapor  in  air 

4.  Latent  heat  of  vaporization  of  water 

5.  Saturation  specific  humidity 

6.  Temperature 

Most  researchers  have  been  more  concerned  with  the  emf  output  as  a  func¬ 
tion  of  water  potential  than  with  the  wet  bulb  depression;  therefore  the 
Peltier  and  Seebeck  coefficients  of  the  metals  involved  must  enter  into  the 
relationships.  Although  the  theoretical  treatments  have  provided  insights 
into  the  response  of  the  instrument  to  various  design  factors,  they  are  of 
minor  interest  in  appl ications.  The  complexities  are  overcome  through  cali¬ 
bration  with  salt  solutions  of  known  water  potential.  However,  the  use  of 
calibration  does  not  eliminate  all  pitfalls  in  the  use  of  TCP. 

Present  usaqe--A  variety  of  instruments  are  presently  being  used  for 
studies  involving  TCPs.  These  are  of  two  types  generally.  One  (Spanner,  1951) 
depends  on  Peltier  cooling  to  condense  a  water  droplet  on  the  wet  junction; 


Rawlins,  S.  L.,  (1966)  "Theory  for  Thermocouple  Psychrometers  Used  to  Measure 
Water  Potential  in  Soil  and  Plant  Samples,"  Anri  cultural  Mctcw>loi)n , 

Vol .  3,  pp.  293-310. 
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the  other  (Richards  and  Ogata,  1958)  requires  manual  placement  of  a  water  drop 
on  the  wet  junction.  The  Peltier  cooled  instrument  was  best  suited  for  the 
present  study  because  the  other  cannot  be  used  for  field  measurements. 

Several  varieties  of  the  Peltier  instrument  exist.  In  some  cases,  addi¬ 
tional  terminals  are  installed  in  the  measuring  chamber  to  provide  a  means  of 
temperature  compensation  or  to  avoid  heating  at  the  reference  junction  when 
the  Peltier  current  is  flowing.  These  types  may  mask  errors  that  result  from 
temperature  gradients.  Another  type  has  a  single  sensing  junction  enclosed  in 
a  protective  chamber.  An  additional  temperature-sensing  thermocouple  is  pro¬ 
vided  as  well.  This  instrument  was  selected  for  the  present  work  as  the  best 
compromise  of  capability  and  simplicity.  Therefore,  this  variety,  shown  in 
Figure  4-1,  is  the  TCP  discussed  in  this  paper. 

The  TCP  may  be  operated  in  either  a  psychrometric  or  hygrometric  mode. 

In  both  modes,  a  drop  of  water  is  condensed  on  a  sensing  junction  by  passing 
a  cooling  current  through  it.  The  hygrometric  mode  uses  an  electric  circuit 
to  cool  the  junction  periodically  by  pulses  in  the  proper  direction.  The 
length  of  the  cooling  period  is  adjusted  according  to  the  TCP  output.  By 
entering  a  gain  in  the  system  (obtained  by  calibration),  the  circuit  can  main¬ 
tain  the  drop  of  water  and  thus  measure  the  emf  resulting  from  keeping  the 
wet  junction  at  the  dew  point.  Although  this  instrument  has  been  used  in 
field  environments,  it  was  not  suitable  for  the  present  study  because  of  very 
little  experience  in  field  applications. 

In  the  psychrometric  mode  the  sensing  junction  is  cooled  by  passing  a 
current  through  it  in  the  proper  direction.  Once  the  cooling  is  finished, 
the  current  is  stopped,  and  the  emf  output  measured.  After  a  certain  period 
(usually  2  to  5  s)  has  elapsed,  the  emf  is  recorded.  As  the  water  drop  evapo¬ 
rates,  the  wet  junction  is  held  at  the  dew  point.  When  evaporative  cooling 
drops  off,  the  emf  decays  back  to  a  zero  reading  if  there  is  no  temperature 
difference  between  the  reference  and  sensing  junctions.  By  comparing  this 
emf  to  the  calibration  curve,  a  value  of  water  potential  may  be  obtained. 

Richards,  L.  A.,  and  Ogata,  G.,  (1953)  "Thermocouple  for  Vapor  Pressure 

Measurement  in  Biological  and  Soil  Systems  at  High  Humidity,"  . 

128:1089-1090. 
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Readout  procedures--In  order  to  operate  TCPs,  several  procedures  are 
required.  A  cooling  current  must  be  applied  in  the  proper  direction  to  cool 
the  sensing  junction.  When  it  is  stopped,  a  voltmeter  in  the  range  of  1  to 
30  pV  is  required  to  measure  the  TCP  emf  output.  Most  researchers  agree  that, 
in  order  to  develop  meaningful  procedures  using  TCPs,  an  X-Y  plotter  or  strip 
chart  recorder  should  be  used  to  study  the  emf  versus  time  data  for  the  psy- 
chrometers.  Once  some  experience  is  gained,  the  necessity  of  plotting  these 
outputs  diminishes. 

For  this  study,  an  MJ-55  thermocouple  microvoltmeter  was  used,  a  Keithly 
Model  155  Microvoltmeter  modified  and  sold  by  Wescor,  Inc.  Modifications 
include  a  cooling  current  circuit  and  manual  switch,  a  0°C  thermocouple  refer¬ 
ence  junction,  and  connectors  for  use  with  Wescor  instruments.  The  emf-time 
data  were  plotted  using  a  Hewlett-Packard  Model  7046A  X-Y  recorder.  The  time 
was  recorded  at  10  s/in,  and  the  output  at  0.5  V/in.  The  MJ-55  provides  a 
±1-V  output  for  full-scale  positive  or  negative  deflection  at  the  terminal 
strip,  through  internal  amplifiers,  for  driving  the  plotter. 

The  TCPs  used  were  chromel -constantan  thermocouples  with  a  single  sensing 
junction.  The  Wescor,  Inc.,  model  PT-51  and  P51  instruments  were  used  with 
lead  wire  lengths  from  1.5  to  6.1  m  (5  to  20  ft).  The  difference  between 
these  instruments  is  the  addition  of  a  copper-constantan  thermocouple  for 
temperature  measurements  on  the  PT  model.  These  instruments  are  covered  with 
a  ceramic  bulb  that  forms  the  sample  chamber.  Another  type  of  TCP  used  was 
the  No.  74-13  produced  by  J.R.D.  Merrill  Specialty  Equipment.  These  are 
screen-covered  instruments  with  a  single  sensing  junction.  An  additional  con¬ 
stantan  wire  is  connected  to  one  reference  junction  to  provide  temperature 
measuring  capability.  They  are  covered  with  a  double  stainless  steel  wire 
mesh. 

The  instruments  were  connected  manually  to  the  MJ-55.  The  meter  was 
turned  on  and  zeroed.  (A  large  deflection  of  the  meter  indicates  the  lack 
of  an  equilibrium  condition  between  the  junctions.)  The  recorder  was  zeroed 
and  started.  The  cooling  current  was  applied  using  the  manual  switch  on  the 
MJ-55.  After  some  experimenting,  a  15-s  cooling  period  was  adopted  as  stand¬ 
ard. 
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It  is  important  to  select  a  uniform  procedure  for  reading  the  thermocouple 
output.  Plots  of  output  versus  time  were  recorded  arid  studied  to  select  a 
method.  Consistent  results  were  obtainei  using  a  2-  to  5-s  time  delay.  The 
delay  begins  when  the  cooling  current  is  shut  off  and  ends  when  the  reading 
is  made.  Although  a  variety  of  output  curve  forms  were  observed,  this  pro¬ 
cedure  provided  consistent  results.  Typical  output  plots  are  shown  in  Figure 
4-2  and  4-3  for  high-  and  low-range  measurements  respectively. 

Calibration  model i ng--The  present  state  of  the  art  for  using  TCP  sensors 
has  been  totally  dependent  on  instrument  calibration.  Despite  several  attempts 
to  theoretically  predict  the  output,  calibration  is  still  an  essential  part  of 
applications  to  real  measurement  problems.  The  emf  will  vary  as  the  relative 
humidity  changes.  Normally  the  TCP  output  emf  is  plotted  versus  water  poten¬ 
tial.  Figure  4-4  illustrates  a  typical  calibration  curve  using  several  cooling 
periods  (Bingham  et  al ,  1972).  Figure  4-5  illustrates  calibration  variation 
with  ambient  temperature  (Meyn  and  White,  1972). 

Because  both  water  potential  and  emf  output  vary  with  temperature,  it  is 
necessary  to  evaluate  the  influence  of  temperature  changes  for  all  applications 
in  which  it  is  not  closely  controlled.  Meyn  and  White  (1972)  have  proposed  a 
model  for  TCP  calibration  that  accounts  for  the  influence  of  temperature  as 
fol lows: 


Y  =  b0x  >  +  btXi  +  b2x:;  +  b 3x ] r  +  b4X2‘i  +  b t. x i x 2 


(4-1) 


where 


Y  =  predicted  water  potential,  bars 
xo  =  dummy  variable  (1) 
xi  =  emf,  iiV 
X2  =  temperature,  °C 

bo,  b i ,  b2,  b  j,  b ,  b;.  -  regression  constants 


Bingham,  G.  E.,  Johnson,  M.  N.,  and  Lemon,  E.  R.  ,  (1972  )  "Influence  of  Heat 
Sink  Design  on  Thermocouple  Psychrometer  Response,"  in  :  romei  ei:  ;n 

Alter  He  Lit  tom:  Research,  Utah  Agricultural  Experiment  Station,  Utah 
State  University,  Logan,  Utah. 

Meyn,  R.  L.  and  White,  R.  S.,  (1972)  "Calibration  of  Thermocouple  Psychroineters : 
A  Suggested  Procedure  for  Development  of  a  Reliable  Predictive  Model," 
in  I'syiihirome try  in  Wu ter  HelutLw  /t.  iseureh,  Utah  Agricultural  Experiment 
Station,  Utah  State  University,  Logan,  Utah,  pp.  56-64. 
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TIME,  s 

a.  SUCTION  =  4.75  pF  (5500  kPa) 


FIGURE  4-2.  TYPICAL  THERMOCOUPLE  PSYCHROMETER  OUTPUT,  HIGH  RANGE. 
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OUTPUT,  uV  OUTPUT 


TIME,  S 

b.  SUCTION  =  3.85  pF  (700  kPa) 


FIGURE  4-3.  TYPICAL  THERMOCOUPLE  PSYCHROMETER  OUTPUT,  LOW  RANGE. 
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OUTPUT,  yV 
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0  TO  20  30  40  50 


FIGURE  4-5.  CALIBRATION  VARIATION  WITH  AMBIENT  TEMPERATURE. 
(AFTER  MEYN  AND  WHITE,  1972) 
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Slack  (1975)  and  Riggle  (1978)  have  followed  the  same  approach.  The  results 
obtained  by  Riggle  and  Meyn  and  White,  along  with  those  of  the  present  study, 
are  shown  in  Table  4-1. 


TABLE  4-1.  REGRESSION  MODELS 


Riggle  (1978) 

Meyn  and  White  (1972) 

. 

NMERI 

(dummy)  b0 

2.292 

0.10617 

-16.4704 

(xi)  bi 

-4.138 

-4.6283 

-  4.0883 

(x2)  b2 

— 

0.0778 

1.6769 

:*i)2  b3 

-0.0557 

-0.0076 

0 

(x2)2  b„ 

-0.0075 

-0.0413 

( x i x 2 )  bs 
coefficient  (R2) 

0.1110 

0.1036 

0.0728 

of  determination 

0.9756 

0.9623 

0.9630 

In  these  experiments ,  a  model  involving  fewer  terms  was  developed  in  each 
case.  Note  that  the  calibration  is  valid  only  if  the  measurements  are  made  in 
exactly  the  same  manner  as  the  calibrations. 

Meyn  and  White  developed  their  model  for  use  with  a  large  number  of 
instruments  whereas  Riggle  calibrated  each  instrument  individually.  In  the 
NMERI  calibrations,  a  single  satisfactory  model  was  sought.  Development 
efforts  commenced  with  the  calibration  of  numerous  instruments  at  ambient  tem¬ 
perature  in  the  laboratory.  This  was  followed  by  additional  calibrations  in 
an  environmental  chamber  at  temperatures  of  roughly  10  and  15°C.  These  data 
were  then  evaluated  using  a  stepwise  multiple  regression  computer  program.  The 
form  of  the  model  was  identical  to  those  previously  used. 


Riggle,  F.  R.,  (1978)  "Soil  Water  Potential  Determinations  with  Thermocouple 
Psychrometers,"  Master  of  Agriculture  Integrating  Paper,  University  of 
Minnesota,  St.  Paul,  Minnesota. 

Slack,  D.  C.,  (1975)  "Modeling  the  Uptake  of  Soil  Water  by  Plants,"  Disserta¬ 
tion,  University  of  Kentucky,  Lexington,  Kentucky. 


Field  installations— At  the  Gallup  airport,  two  boreholes  were  made  to  a 
depth  of  6.1  m  (20  ft)  for  installation  of  the  instruments.  This  depth  was 
expected  to  include  the  soil  active  zone.  On  the  recommendations  of  an 
expert  on  use  of  the  instruments,  redundancy  was  included  to  provide  for 
failure  of  the  units.  Table  4-2  shows  the  instrumentation  included  in  each 
borehole.  Two  types  of  instruments  were  used.  The  instruments  with  the  model 
designation  beginning  with  P  were  obtained  from  Wescor,  Inc.,  in  Logan,  Utah. 
For  these  instruments,  the  P  indicates  psychrometer;  T  means  that  an  addi¬ 
tional  temperature-measuring  thermocouple  is  included;  and  numbers  are  the 
lead  wire  length  in  feet.  The  others,  numbered  in  the  200 ' s  in  Table  4-2, 
were  obtained  from  the  J.R.D.  Merrill  Specialty  Co.  in  Logan,  Utah.  These 
instruments  included  a  temperature-measuring  thermocouple  and  a  3-m  (9.8-ft) 
lead  wire. 

Initially  two  boreholes  were  drilled  at  the  DFW  airport  (borings  5  and  6). 
Because  of  the  poor  performance  of  these  instruments,  five  additional  psychro- 
meters  were  installed  in  another  site  (boring  8).  Table  4-3  shows  the  instru¬ 
mentation  used. 

Because  of  the  climate  at  the  Jackson  airport,  the  equilibrium  suction 
below  the  depth  of  the  active  zone  was  estimated  to  be  too  low  for  reliable 
performance  at  the  depths  used  for  DFW  and  GAL.  Therefore,  instruments  were 
installed  only  to  1.5  m  (5  ft)  for  measurements  in  the  active  zone  of  the 
soil  (Table  4-4). 

Instrument  performance  was  monitored  to  evaluate  the  redundancy  designed 
into  the  instrument  configuration.  The  percentage  of  instruments  doubled  was 
as  follows:  DFW  boring  5,  67;  DFW  boring  6,  57;  DFW  boring  8,  0;  GAL 
boring  1,  50;  GAL  boring  2,  50;  JSN  boring  1,  0.  The  GAL  installations  were 
designed  with  less  redundancy  because  of  the  dry  climate  and,  therefore, 
favorable  measurement  conditions.  More  redundancy  was  used  for  the  DFW  sites 
because  of  wetter  conditions  and,  therefore,  more  difficult  operating  environ¬ 
ment.  The  instruments  of  JSN  were  not  duplicated  because  of  the  ease  of 
replacing  instruments  to  a  depth  of  1,5  m  (5  ft).  The  total  percent  of  avail¬ 
able  instruments  at  each  site  was  as  follows:  DFW  boring  5,  83;  DFW  boring  6, 
36;  DFW  boring  8,  100;  GAL  boring  1,  69;  GAL  boring  2,  81;  JSN  boring  1,  100. 
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TABLE  4-2.  INSTRUMENTATION  AT  GAL 


Boring  1 

Boring  2 

Depth, 
m  (ft) 

Sensor 

number 

Type 

Depth , 
in  (ft) 

Sensor 

number 

Type 

0.61  (2.0) 

8 

PT51-05 

0.5  (1.5) 

3 

PT51-05 

0.61  (2.0) 

2 

PT51-05 

0.5  (1.5) 

4 

PT51-05 

1.2  (4.0) 

22 

PT51-10 

1.1  (3.5) 

24 

PT51-10 

1.2  (4.0) 

203 

74-13 

1.1  (3.5) 

204 

74-13 

1.8  (6.0) 

21 

PT51-10 

1.7  (5.5) 

23 

PT51-10 

1.8  (6.0) 

25 

PT51-10 

1.7  (5.5) 

27 

PT51-10 

2.1  (7.0) 

201 

74-13 

2.0  (6.5) 

203 

74-13 

2.4  (8.0) 

108 

PT51-16 

2.3  (7.5) 

109 

PT51-16 

3.0  (10.0) 

110 

P51-16 

2.9  (9.5) 

112 

P51-16 

3.0  (10.0) 

111 

PT51-16 

2.9  (9.5) 

107 

PT51-16 

3.7  (12.0) 

115 

PT51-16 

3.5  (11.5) 

113 

P51-16 

4.3  (14.0) 

116 

P51-16 

4.1  (13.5) 

114 

PT51-16 

4.9  (16.0) 

PT51-20 

4.7  (15.5) 

103 

PT51-20 

4.9  (16.0) 

104 

P51-20 

4.7  (15.5) 

105 

P51-20 

5.5  (18.0) 

P51-20 

5.2  (17.5) 

106 

PT51-20 

TABLE  4-4.  INSTRUMENTATION  AT  JACKSON  AIRPORT 


Boring  No. 

3 

Depth, 
m  (ft) 

Sensor 

number 

Type 

0.3  (1.0) 

1 

PT51-05 

0.6  (2.0) 

33 

PT51-10 

0.9  (3.0) 

130 

PT51-16 

1.2  (4.0) 

38 

PT51-10 

1.5  (5.0) 

39 

PT51-10 

The  number  available  refers  to  the  number  of  instruments  in  which  the  circuit 
was  intact.  However,  the  availability  of  an  instrument  did  not  guarantee  that 
a  measurement  was  made.  In  some  cases,  the  soil  was  too  wet  for  TCP  opera¬ 
tion.  When  the  thermocouples  corroded  or  otherwise  failed,  it  was  easily 
detected  by  measuring  circuit  resistance.  Results  of  measurements  using  these 
instruments  are  discussed  later  in  this  section. 

Filter  Paper  Method 

Researchers  in  the  field  of  soil  science  recognized  that  forces  attracting 
water  to  paper  could  be  used  to  estimate  the  magnitude  of  soil  attraction  for 
water.  Hansen  (1926)  used  blotting  paper  as  a  carrier  for  sugar  solutions  in 
work  at  the  University  of  Copenhagen.  Strips  saturated  with  four  different 
solutions  were  exposed  to  soils  in  closed  containers.  The  one  with  the  least 
moisture  change  at  equilibrium  was  interpreted  as  being  nearest  the  stress 
level  of  the  soil.  Stocker  (1930)  improved  the  accuracy  of  the  method  by 
employing  more  solutions.  Another  improvement  was  made  by  Gradmann  (1934)  who 
calibrated  the  blotting  paper  for  weight  versus  stress  after  soaking  in  salt 
solutions.  The  sensors  were  then  equilibrated  with  soil  samples,  and  water 


Gradmann,  H.,  (1934)  "Uber  die  Messung  von  Bodensaugwerten"  (About  the  measure¬ 
ment  of  soil  suction),  Jh.  Wise  Bot.  80:92-111. 

Hansen,  H.  C.,  (1926)  "The  Water  Retaining  Power  of  the  Soil,"  Journal  of 
Ecology,  Vol .  14:111-119. 

Stocker,  0.,  (1930)  "Uber  die  Messugn  von  Bodensaughkraften  and  ihren  Verhalt- 
nis  zu  den  Wurzelsaughkraften,"  Z.  Bot.  23:27-56. 
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contents  were  used  to  infer  moisture  stress.  With  some  minor  changes,  Eckhardt 
(1960)  used  this  method  in  France. 


The  first  use  of  paper  as  a  moisture  stress  sensor  without  a  hygroscopic 
salt  was  reported  by  Gardner  (1937).  He  used  an  ash-free  quantitative  filter 
paper  (Schleicher  and  Schuell,  No.  589,  White  Ribbon)  calibrated  for  water 
content  versus  stress.  From  this  beginning,  numerous  applications  have  fol¬ 
lowed,  leading  to  the  present  method  used  extensively  and  routinely  by  the 
Water  Resources  Division  of  the  U.S.  Geological  Survey  (USGS).  The  first  use 
of  this  method  for  civil  engineering  purposes  was  by  McKeen  (1977). 

Cal ibration--Numerous  calibrations  have  been  made  for  different  papers 
used  in  soil  moisture  studies.  The  calibration  techniques  have  generally 
varied  depending  on  the  range  of  interest.  Soil  samples  are  prepared  at  a 
known  moisture  stress  or  suction  level,  then  the  filter  paper  is  equilibrated 
with  it  for  calibration.  At  the  higher  pressures,  filter  papers  may  be 
equilibrated  on  a  pressure  membrane  or  above  salt  or  acid  solutions  of  known 
concentration.  In  this  way  a  full  calibration  curve  can  be  developed. 

McQueen  and  Miller  (1974)  found  that  a  two-part  relationship  best  fit 
their  calibration  data.  They  determined  that  the  upper  segment  (pF  >  2.5) 
represents  moisture  retained  as  film  adsorbed  to  particle  surfaces  and  that 
the  lower  segment  (pF  <  2.5)  represents  moisture  retained  by  capillary  or  sur¬ 
face  tension  forces.  It  is  no  coincidence  that  this  level  (2.5  pF)  of  moisture 
stress  corresponds  to  the  field  capacity  of  soils  as  found  in  the  agricultural 
literature.  In  an  extensive  and  thorough  study  by  McQueen  and  Miller  of  cali¬ 
bration  and  factors  influencing  it,  a  range  of  calibration  relationships  was 


Eckardt,  F.  E.,  (1960)  "Eco-physiological  Measuring  Techniques  Applied  to 
Research  on  Water  Relations  of  Plants  in  Arid  and  Semi-arid  Regions,  in 
Avid  Zone  Res. :  XV.  UNESCO. 

Gardner,  R.  (1937),  "A  Method  of  Measuring  the  Capillary  Tension  of  Soil  Mois¬ 
ture  over  a  Wide  Moisture  Range,"  Roil  Science,  Vol .  43:277-283. 

McKeen,  R.  G.,  (1977)  "Characterizing  Expansive  Soils  for  Design  Purposes," 
presentation  at  the  joint  meeting  of  Texas,  New  Mexico,  and  Mexico  sec¬ 
tions,  American  Society  of  Civil  Engineers,  Albuquerque,  New  Mexico, 
October  6-8. 

McQueen,  I.S.  and  Miller,  R.  F.,  (1974)  "Approximating  Soil  Moisture  Chara^ 
istics  from  Limited  Data:  Empirical  Evidence  and  Tentative  Mode1 
Resources  Research,  Vol.  10,  No.  3:521-527. 


defined  (Figure  4-6),  The  calibrations  for  the  USGS  work  and  for  the  NMER1 
measurements  are  also  shown  in  Figure  4-6.  The  NMERI  curve  falls  lower  than 
the  USGS  curve  because  the  NMERI  calibrations  used  salt  solutions  in  the  range 
of  3.0  to  4.7  pF.  The  calibration  was  intended  for  use  with  clay  soils  in 
this  range  of  suctions  and  was  therefore  deemed  satisfactory.  Subsequent  cali¬ 
brations  at  the  low  end  and  communications  with  the  developers  of  the  USGS 
method  substantiated  the  NMERI  calibration  relationship.* 

Procedure- -The  method  adopted  for  this  study  was  essentially  that  pub¬ 
lished  by  McQueen  and  Miller  (1968).  The  procedure  uses  filter  paper  discs 
5  cm  in  diameter.  In  this  study,  Schleicher  and  Schuell  No.  589  White  Ribbon 
was  used.  As  shown  by  McQueen  and  Miller  (1977),  virtually  any  uniform 
quality  paper  may  be  used.  Air-dried  filter  papers  were  exposed  to  salt  solu¬ 
tions  (potassium  chloride)  for  a  period  of  seven  days.  The  papers  absorbed  a 
small  amount  of  water  from  the  salt  solution.  The  papers  were  quickly  and 
carefully  transferred  into  small  aluminum  weighing  cans  and  weighed  to  the 
nearest  0.0001  g.  The  containers  were  opened  and  oven-dried  overnight.  The 
containers  were  then  closed  inside  the  oven,  removed,  and  placed  on  an  aluminum 
block  to  cool  for  30  s.  The  filter  paper  was  weighed  to  obtain  the  dry  weight. 
In  each  case  a  separate  tare  weight  was  taken  for  both  the  wet  and  dry  weigh¬ 
ings.  The  moisture  stress  of  the  salt  solution  was  plotted  versus  the  filter 
paper  water  content. 

Suction  measurements  on  soil  samples  were  made  using  filter  paper  from 
one  lot.  First,  the  sample  was  placed  in  a  moisture  can.  Then  two  filter 
papers  were  placed  directly  on  top  of  the  soil  sample.  In  most  of  the  samples 
tested  in  this  study,  suctions  were  beyond  the  range  of  capillary  water;  there¬ 
fore,  the  filter  paper  could  not  absorb  water  through  capillary  action.  The 
moisture  transferred  by  vapor  diffusion  and  thus  represents  a  total  suction 
value.  [For  soils  with  suctions  below  the  capillary  limit  (<2.5  pF),  precau¬ 
tions  are  required  to  prevent  capillary  flow.]  Following  equilibration,  both 
filter  papers  were  removed  and  placed  in  separate  moisture  containers.  Their 
respective  water  contents  were  determined.  The  water  content  of  the  bottom 


McQueen,  Irel  S.  and  Miller,  R.  G. ,  (1968)  "Calibration  and  Evaluation  of  a 
Wide-Range  Gravimetric  Method  for  Measuring  Moisture  Stress,"  :'oii 
Science,  V ol.  106,  No.  3:225-321. 

♦Personal  communication  from  R.  F.  Miller  to  the  author,  1977. 
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FILTER  PAPER  WATER  CONTENT  (wf) 


filter  paper  (next  to  the  top  of  the  sample)  was  used  as  the  measure  of  soil 
suction.  The  second  filter  paper  measurement  was  used  as  a  quality  control 
check.  If  the  upper  paper  was  slightly  drier  than  the  bottom  paper,  normal 
equil ibration  occurred,  and  the  measurement  was  good.  Drastic  differences 
(>0.5  pF)  or  a  wetter  top  paper  indicated  an  abnormal  condition.  Abnormal 
results  were  discarded. 

Several  things  created  problems  in  using  the  filter  paper  method.  Occa¬ 
sionally,  condensation  occurs  inside  the  container  as  a  result  of  poor  temper¬ 
ature  control  or  very  wet  soil.  If  the  filter  paper  contacts  it,  this  mois¬ 
ture  will  be  absorbed  rapidly  by  capillary  action.  Great  care  in  transferring 
filter  papers  is  required.  Dry  laboratory  air  can  drastically  alter  the  mois¬ 
ture  content  of  the  filter  paper  sensor  if  exposed  for  a  long  time.  Experi¬ 
ence  in  this  study  indicates  that  technician  training  is  very  important.  In 
most  cases,  errors  result  in  extremely  high  or  low  water  content  of  the  sen¬ 
sors  which  yield  improbable  data.  Thus  an  erroneous  result  is  obvious  when 
the  data  are  examined. 

Field  measurements--In  the  field  experiment,  suction  measurements  were 
made  regularly  on  soil  samples  removed  near  the  locations  of  the  in  situ  mea¬ 
surements.  The  samples  were  taken  by  pushing  a  smal 1 -diameter  tube  into  the 
ground.  The  soil  was  removed  and  sealed  with  two  filter  papers  in  a  moisture 
can.  The  containers  were  returned  to  the  laboratory  and  equilibrated.  The 
samples  were  obtained  at  152-mm  (6-in)  intervals  to  a  depth  of  1.5  m  (5  ft). 
When  soil  borings  were  taken  with  a  truck-mounted  drill  rig,  samples  were  also 
removed  at  intervals  for  suction  determination.  One  advantage  of  the  filter 
paper  method  is  that  large  numbers  of  samples  may  be  taken  at  very  low  cost. 

In  addition  to  the  field  samples,  several  sets  of  data  were  measured  on 
soil  samples  in  the  laboratory  for  evaluation  of  the  method  and  for  comparison 
of  procedures.  Results  are  discussed  later  in  this  section  of  the  report. 

Hydraulic  Tensiometers 

Hydraulic  tensiometers  were  installed  to  measure  the  low-suction  (wet 
season)  range  of  the  field  soils.  Tensiometers  were  installed  at  the  DFW  and 
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JSN  sites.  It  was  anticipated  that  the  weather  conditions  at  the  GAL  site 
would  be  too  dry  for  satisfactory  operation  of  the  instruments.  Figure  4-7 
illustrates  the  design  of  the  instruments  used.  The  ceramic  bulbs  were  placed 
at  depths  of  0.3,  0.6,  0.9,  1.2,  and  1.5  m  (1.0,  2.0,  3.0,  4.0,  and  5.0  ft)  at 
the  sites  studied.  A  water-alcohol  solution  was  used  to  provide  a  continuous 
liquid  column  from  the  bulb  to  the  mercury  container.  The  soil  surrounding 
the  ceramic  bulb  attracted  water  through  the  bulb  by  capillary  action,  pulling 
the  mercury  column  up.  With  appropriate  corrections  for  the  liquid  column 
lengths,  a  measurement  of  the  suction  in  the  soil  was  obtained.  Because  the 
ceramic  bulb  allows  salts  to  pass  through  freely,  this  device  measures  only  the 
matrix  component  of  suction* 

COMPARISON  OF  FP  AND  TCP  MEASUREMENTS 
Laboratory  Comparison 

The  TCP  and  filter  paper  (FP)  methods  of  suction  measurement  were  compared 
and  evaluated.  Two  comparisons  were  made  using  different  procedures.  First, 
TCPs  and  FP  were  exposed  to  the  same  soil  sample  in  the  same  closed  container. 
Four  different  groups  of  tests  were  performed,  totaling  thirty-one  samples. 

The  results  are  shown  in  Figure  4-8.  The  computation  of  the  linear  regression 
confirms  that,  under  the  conditions  of  this  evaluation,  the  TCP  and  the  FP 
methods  yield  the  same  results.  The  second  comparison  was  made  in  an  attempt 
to  verify  the  procedures  reported  by  Snethen  and  Johnson  (1980)  of  Waterways 
Experiment  Station  (WES). 

The  results  of  the  WES  experiment  indicate  that  FP  provides  lower  values 
of  suction  than  the  TCP.  This  difference  in  results  has  been  attributed  to 
two  specific  reasons:  the  equilibrium  period  and  the  sampling  technique.  The 
procedure  used  by  Snethen  and  Johnson  assumes  that  the  TCP  equilibrates  with 
the  air  above  the  soil  sample  in  two  days.  Data  to  substantiate  this  time 
period  selection  are  not  evident  in  the  literature.  In  the  present  study,  the 
TCPs  were  allowed  to  equilibrate  for  seven  days,  the  same  period  of  time  as 
the  FP.  Data  from  the  WES  report  are  plotted  in  Figure  4-9,  along  with  the 
NMERI  regression  curve.  It  is  apparent  that,  for  low-range  measurements,  the 


Snethen,  D.  R.,  and  Johnson,  L.  0.,  (1980)  Evaluation  of  .'oil  u..k  ft\  - 

Filter  Paper ,  Miscellaneous  Paper  GL-80-4,  U.S.  Army  Waterways  Experiment 
Station,  Vicksburg,  Mississippi. 
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MEASURING  TUBE 


FIGURE  4-7.  HYDRAULIC  TENSIOMETER  CONFIGURATION. 
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R2  =  COEFFICIENT  OF  DETERMINATION 
a  =  STANDARD  ERROR 


FIGURE  4-8.  LABORATORY  COMPARISON  -  SAME  SAMPLE. 
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TCP  instruments  in  the  WES  tests  did  not  achieve  vapor  equilibrium  with  the 
soil  samples.  The  porous  ceramic  protective  cover  on  the  TCP  may  not  have  been 
sufficiently  permeable,  at  low  pressure  gradients,  to  allow  equilibration  in 
two  days;  or  the  covers  may  have  been  contaminated,  preventing  free  diffusion 
of  water  vapor. 

The  interference  of  the  ceramic  cover  with  equil ibration  is  called  lag. 
Data  were  gathered  in  the  laboratory  to  obtain  an  estimate  of  this  lag.  Fig¬ 
ure  4-10a  shows  that  equilibration  is  rapid  in  calibration  solutions.  In  Fig¬ 
ure  4-10b,  equilibration  is  shown  with  a  soil  sample.  The  minimum  equilibra¬ 
tion  time  is  about  6  days. 

The  WES  sampling  technique  also  differed  from  the  NMERI  procedures.  In 
the  WES  technique,  an  undisturbed  soil  sample  was  divided  into  two  pieces,  one 
for  each  type  of  measurement.  In  the  NMERI  verification  tests,  the  WES  pro¬ 
cedures  of  time  and  sampling  were  followed.  Shelby  tube  samples  were  cut  in 
slices  approximately  25  mm  (1  in)  thick.  These  slices  were  then  divided  into 
two  parts,  one  for  suction  measurement  using  FP,  one  for  measurement  using 
TCP.  The  samples  were  allowed  to  equilibrate  for  two  days.  The  linear 
regression  was  computed  and  compared  with  the  regression  obtained  with  the 
NMERI  procedures.  The  comparison  is  shown  in  Table  4-5.  These  data  clearly 
show  the  dramatic  influence  of  equilibration  time  period  and  sample  variation, 
even  though  the  samples  were  adjacent  in  the  soil  profile. 

TABLE  4-5.  EFFECTS  OF  SAMPLE  VARIATION  AND 
EQUILIBRATION  TIME 


Samp! ing 

Number 

Coefficient  of 

Standard 

method 

of  tests 

determination 

error 

Same  sample 
for  7  days 

31 

0.865 

0.19 

Different 

samples 
for  2  days 

53 

0.624 

0.45 
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a.  CALIBRATION  WITH  SALT  SOLUTION 


b.  CALIBRATION  WITH  SOIL. 


FIGURE  4-10.  ILLUSTRATION  OF  TCP  EQUILIBRIUM  TIME 


From  the  present  study  it  is  concluded  that,  given  full  equil ibration,  no 
sample  variance,  and  laboratory  conditions,  the  FP  and  TCP  methods  are  equiva¬ 
lent  for  the  determination  of  soil  moisture  suction  in  the  range  3.2  to  4.8  pF. 
Because  this  range  covers  most  soils  of  interest  in  expansive  soil  studies, 
both  methods  can  be  applied  to  these  studies.  The  choice  of  method  should  be 
based  on  other  considerations;  the  precision  appears  to  be  equivalent.  It  is 
of  primary  importance  that  the  TCPs  be  allowed  to  reach  equilibrium. 

Comparison  of  In  Situ  Samples 

Borings  to  a  depth  of  6  m  (20  ft)  were  made  at  each  site  to  obtain  samples 
for  laboratory  characterization  studies.  Suction  profiles  were  measured  with 
FP  and  TCP  in  each  case.  A  total  of  six  borings  were  made:  three  at  DFW, 
two  at  JSN,  and  one  at  GAL. 

The  results  of  the  FP  and  TCP  measurements  of  the  in  situ  suction  profile 
from  these  borings  are  shown  in  Figure  4-11  for  the  DFW  site.  These  data 
require  some  explanation.  The  ideas  presented  here  evolved  throughout  the 
field  experiment  and  should  be  placed  in  chronological  order  for  clarity. 

Boring  DFW  4  was  made  in  October  1978  to  determine  the  suitability  of  the  site 
for  the  field  experiment.  Two  suction  profiles  were  measured  in  the  labora¬ 
tory  on  separate  samples,  one  with  TCP  and  one  with  FP.  At  that  time  the  FP 
was  being  used  on  a  trial  basis.  The  TCP  was  considered  a  reliable  method 
for  suction  determination.  As  shown  in  Figure  4-11,  dramatic  differences 
were  noted  between  the  two  methods  at  depths  below  1.8  m  (6  ft).  It  was  ini¬ 
tially  assumed  that  the  FP  data  were  erroneous.  The  TCP  data  indicated  suc¬ 
tions  above  3.4  pF,  indicating  no  possibility  of  free  water  or  a  water  table 
in  the  profile.  The  site  was  accepted  as  suitable  for  the  field  experiment. 

Within  about  30  days,  borings  were  made  for  TCP  installation.  Before  the 
installation  could  be  completed,  free  water  entered  the  borings  and  rose  to 
within  1.8  to  2.4  m  (6  to  8  ft)  of  the  surface.  The  undisturbed  tube  samples 
from  boring  4,  as  well  as  those  for  TCP  installation  (DFW  5  and  6),  were  thor- 
ougly  studied.  The  study  revealed  the  presence  of  numerous  sand  seams  in  the 
profile  (Figure  4-12).  It  was  found  that  the  FP  suction  measurements  made  on 
samples  from  DFW  4  had  more  accurately  depicted  the  actual  suction  distribu¬ 
tion  in  the  profile.  The  TCPs  evidently  did  not  reach  equilibrium. 
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The  results  on  samples  from  DFW  10  (Nov.  1979)  confirm  the  sharp  gradients 
in  suction  at  several  sand  seams.  During  the  sampling  for  DFW  10  there  was  no 
problem  locating  the  sand  seams.  Suction  samples  were  taken  above,  below,  and 
within  the  seams  at  several  points.  It  is  clear  that  the  water  content  in  the 
profile  gives  virtually  no  information  about  moisture  distribution  and  equi¬ 
librium  in  the  soil  profile,  except  in  conjunction  with  suction  measurements. 

The  water  content  profile  from  DFW  10  indicates  that  the  sand  seam  at 
2.1  m  (7  ft)  prevents  water  from  below  from  moving  up  in  the  profile.  The 
apparent  dryness  of  DFW  10  data  compared  to  DFW  4  is  believed  to  be  normal 
profile  variation.  The  two  borings  were  located  about  15.2  m  (50  ft)  apart. 

The  difference  of  0.2  pF  (3.7  to  3.5  pF)  is  well  within  the  variation  found  for 
sample-to-sample  differences. 

A  bentonite  layer  occurred  at  3  m  (9  ft)  below  the  surface.  Inspection 
of  suction  and  water  content  at  this  depth  clearly  show  that  suction  is  an 
indicator  of  moisture  condition  whereas  water  content  is  not.  The  greater 
clay  content  (98  percent  <  2  pm)  of  this  layer  requires  a  very  high  water  con¬ 
tent  relative  to  surrounding  clay. 

Soil  samples  were  removed  at  the  time  in  situ  TCP  measurements  were  made 
at  each  site.  These  samples  were  sealed  in  the  field  into  containers  with  FP 
sensors,  transported  to  the  laboratory,  and  equilibrated.  Data  showing  in  situ 
TCP  measurements  and  FP  measurements  on  samples  from  the  same  depth  are  plotted 
in  Figure  4-13.  A  regression  computed  for  these  field-lab  data  and  the  regres¬ 
sions  for  the  previous  laboratory  comparison  are  shown.  Observations  made 
.from  these  data  are  as  follows: 

1.  The  standard  deviation  is  less  for  these  tests  than  for  the  laboratory 
tests  on  different  samples  (using  WES  procedure). 

2.  The  deviation  from  a  one-to-one  slope  again  points  to  lower  suctions 
for  FP  on  the  low  end  of  the  suction  range  measured. 

3.  FP  measurements  of  suctions  made  in  the  laboratory  on  field  samples 
are  at  least  as  accurate  as  in  situ  TCP  measurements. 

The  slope  of  the  field-lab  regression  line  indicates  a  lower  FP  result 
as  compared  to  the  TCP  on  the  low  end  of  the  range.  Two  causes  are  respon¬ 
sible  for  the  discrepancy:  lag  caused  by  the  ceramic  cover,  and  installation 
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methods.  The  instruments  were  installed  in  the  field  by  boring  a  hole,  plac¬ 
ing  the  instruments,  and  recompacting  the  soil.  Thus,  the  field  permeability 
was  reduced  at  the  point  of  measurement,  resulting  in  increased  lag  time. 

RESULTS 

Field  Profiles 

In  situ  thermocouple  psychrometers--TCPs  were  installed  in  boreholes  at 
each  site  to  measure  in  situ  suction  profiles.  Two  boreholes  were  instrumented 
at  each  of  the  GAL  and  DFW  sites;  a  single  installation  was  made  at  the  JSN 
site.  These  TCPs  were  monitored  periodically  to  obtain  suction  estimates 
throughout  the  year.  The  profiles  are  illustrated  in  Figures  4-14  to  4-16. 
Figure  4-14  shows  that  the  suction  at  the  GAL  site  is  very  stable  below  1.5  m 
(5  ft).  Water  apparently  seeped  down  the  side  of  borehole  2,  affecting 
measurements  at  1.1  m  (3.5  ft)  and  1.7  m  (5.5  ft). 

The  DFW  data  in  Figure  4-15  shows  evidence  of  seasonal  effects  to  the 
1.8-m  (6-ft)  depth.  The  presence  of  sand  seams  which  periodically  contained 
water  at  depths  below  1.8  m  (6  ft)  produced  an  erratic  suction  profile.  How¬ 
ever,  data  obtained  from  this  site  with  the  FP  method  corroborated  the  TCP 
measurements.  Most  of  the  instrumentation  at  DFW  was  infiltrated  when  the 
sand  seams  filled  with  water.  At  DFW  boring  5,  for  example,  several  instru¬ 
ments  were  operational,  but  conditions  were  too  wet  to  allow  a  measurement. 
Apparently  the  relative  humidity  was  too  high,  preventing  evaporative  cooling. 
In  addition,  the  TCPs  were  too  widely  spaced  to  depict  suction  in  a  soil  pro¬ 
file  as  discontinuous  and  complex  as  the  DFW  site. 

The  data  for  the  JSN  site,  shown  in  Figure  4-16,  were  limited  because  of 
the  range  restrictions  of  the  TCP.  Although  the  instruments  did  not  fail,  the 
soil  was  too  wet  for  reliable  measurements  by  TCPs.  The  suitability  of  TCP 
for  in  situ  measurements  depends  on  the  expected  range  of  suction.  Profiles 
that  remain  below  3.5  pF  should  not  be  instrumented.  Profiles  between  3.5  and 
4.0  pF  are  satisfactory;  those  above  4.0  pF  appear  to  be  the  optimum  for  TCP 
measurements. 
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Laboratory  filter  paper  tests--Further  descriptions  of  suction  variation 
at  OFW  are  shown  in  Figure  4-17.  These  profiles  were  determined  in  the  labora¬ 
tory  using  the  FP  method.  The  dashed  lines  represent  the  profiles  determined 
from  the  TCPs  as  shown  in  Figure  4-15.  The  data  for  January  indicate  that  the 
TCPs  at  0.3  m  (1  ft)  and  0.9  m  (3  ft)  were  too  wet  to  respond.  The  results  in 
April,  September,  and  November  show  transition  periods;  May  and  July  are  at 
end  points  of  the  wet  and  dry  cycle  respectively.  In  the  transition  periods, 
water  is  either  entering  or  leaving  the  profile.  Because  of  variations  in 
drainage  of  the  water,  suction  variations  are  present;  and  TCP  lag  causes  the 
discrepancy  between  TCP  and  FP  measurements.  As  the  wetting  or  drying  influ¬ 
ence  persists,  the  soil  approaches  a  uniform  suction  condition.  In  the  May 
and  July  profiles  this  condition  is  clearly  demonstrated. 

Data  in  Figure  4-18  indicate  the  variety  of  suction  profiles  obtained 
during  the  transition  period  from  dry  to  wet,  the  recharge  period  in  the 
Thornthwaite  system.  The  four  plots  represent  suction  versus  depth  at  four 
different  points  horizontally  separated  by  less  than  15  m  (50  ft).  The  rate 
of  change  of  suction  with  elevation  change  is  a  good  indicator  of  the  rate  of 
moisture  movement  into  or  out  of  the  profile.  Thus,  the  following  observa¬ 
tions  can  be  made  from  Figure  4-18. 

1.  DFW  sites  1,  2,  and  4  were  wetting  at  about  the  same  rate. 

2.  At  the  zone  for  0.6  to  0.9  m  (2  to  3  ft)  to  1.2  m  (4  ft),  all  profiles 
show  the  effect  of  a  prolonged  wetting  period.  Weather  data  indicate  that,  20 
to  30  days  before  sampling,  8.6  cm  (3.4  in)  of  rain  occurred  over  a  period  of 
15  days. 

3.  The  field  permeability  and  surface  drainage  characteristics  are  cap¬ 
able  of  producing  dramatically  different  suction  profiles  over  short  horizontal 
distances  during  wetting  or  drying  phases.  Prolonged  wetting  or  drying  periods 
will  produce  more  uniform  suction  conditions. 

FP  suction  profiles  for  GAL  are  shown  in  Figure  4-19.  The  dashed  line 
shows  the  TCP  measurements  (from  Figure  4-14).  As  before,  the  data  from  May 
and  July  show  good  agreement,  but  the  data  for  November  '78,  March  '79,  and 
November  '79  do  not  agree  well.  Thus  the  data  show  the  inability  of  the 
in  situ  TCP  to  track  the  changing  suction  profile. 
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FIGURE  4-17.  COMPARISON  OF  TCP  AND  FP  SUCTION  PROFILES  AT  DFW. 
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FIGURE  4-18.  VARIATION  OF  SUCTION  DURING  RECHARGE  AT  DFW 


69 


SUCTION.  kPa 


FIGURE  4-19.  COMPARISON  OF  TCP  AND  FP  SUCTION  PROFILES  AT  GAL. 


Results  of  measurements  at  the  JSN  site  are  shown  in  Figure  4-20.  Note 
that  most  of  the  FP  suction  values  are  less  than  3.0  pF  in  the  top  0.76  m 
(2.5  ft)  of  the  soil.  Volume  changes  in  soil  at  these  levels  of  suction  are 
not  significant.  It  is  also  clear  that  the  clay  soil  from  about  0.76  m 
(2.5  ft)  down  remained  at  a  stable  suction  between  3.0  and  4.0  pF  throughout 
the  observations. 

Figure  4-21  shows  FP  data  obtained  from  samples  taken  to  the  6.1-m  (20-ft) 
depth  at  JSN.  The  first  boring  (Jan.  *79)  was  at  a  location  about  91.5  m 
(100  yd)  from  the  site  because  the  site  was  inaccessible  during  rains.  The 
profile  shows  ponding  water  at  the  surface  resulting  in  wetting  of  the  gray 
clay  layer.  However,  once  the  wetting  reached  the  interface  with  the  next 
clay  layer,  the  wetting  effect  disappears.  The  other  6.1-m  (20-ft)  boring 
(Dec.  ’79)  was  made  near  the  site  of  the  instrumentation.  These  data  show  no 
change  with  depth.  The  surface  differences  were  caused  by  the  elevation  and 
drainage  at  the  site. 

It  appears  that  the  interface  between  clay  layers  affected  the  movement 
of  water  in  the  profile.  Such  features  of  the  soil  profile  may  provide  access 
for  water  to  infiltrate  beneath  pavements  even  when  the  pavement  surface  is 
properly  maintained  and  provides  a  seal  to  moisture.  The  identification  of 
these  avenues  of  water  movement  must  be  a  part  of  the  investigation  of  the 
soil  condition  at  a  potential  site. 

Moisture  Characteristics 

The  soil  moisture  characteristic  defines  the  relationship  between  suction 
(or  free  energy  of  the  soil  water)  and  water  content.  Its  value  is  largely  a 
function  of  clay  type  and  amount.  Figures  4-22  to  4-24  show  suction  versus 
water  content  data  for  samples  obtained  in  the  field  experiments  at  DFW,  GAL, 
and  JSN.  These  data  are  for  samples  taken  in  the  top  1.5  m  (5  ft)  of  the  soil 
profile.  In  Figure  4-22,  DFW  data  provide  a  moisture  characteristic  (m)  of 
0.082  for  the  central  line  shown.  An  abrupt  change  is  apparent,  beginning 
about  3.5  pF.  At  this  point  the  soil  response  to  water  content  changes  is 
altered,  and  small  increases  produce  relatively  large  reductions  in  suction. 
This  aggregation  behavior  was  discussed  previously  in  McKeen  and  Nielsen 
(1978). 
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4-20.  COMPARISON  OF  TCP  AND  FP  SUCTION  PROFILES  AT  JSN. 
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FIGURE  4-23.  MOISTURE  SUCTION  DATA  FOR  GAL. 
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In  Figure  4-23  the  GAL  data  indicate  a  similar  behavior,  except  that  the  sur¬ 
face  material  [top  15  cm  (6  in)]  is  clearly  different  (m  =  0.092)  from  the 
bulk  of  the  soil  (m  =  0.042). 

At  JSN,  a  similar  result  was  obtained  because  of  the  alluvial  material 
at  the  surface.  For  the  surface  material,  m  =  0.135;  the  remainder  of  the 
profile  is  clay  (m  =  0.085).  Transitions  of  samples  from  one  group  of  data  to 
the  other  were  noted.  The  explanation  for  such  occurrences  is  that  the  wetting 
phenomena  reduces  the  suction.  As  suction  is  reduced,  the  soil  structure  is 
remolded  and  altered.  Because  of  the  structural  alterations,  the  soil  loses 
bearing  capacity,  exhibits  lower  density,  and  is  not  very  responsive  to  further 
suction  reduction.  The  soil  will  not  swell  under  loads  and  may  collapse  or 
consolidate,  driving  out  the  capillary  water  which  is  at  a  much  lower  suction 
level.  It  is  evident  that  field  soils  under  overburden  loads  will  not  reach 
suctions  as  high  as  those  of  laboratory  tests. 

It  is  also  apparent  the  field  soils  do  not  exhibit  any  significant  hyster¬ 
esis  caused  by  the  cycle  of  wetting  and  drying.  At  a  given  level  of  suction, 
there  is  little  difference  between  water  contents  of  the  wetting  and  the  dry¬ 
ing  soil.  This  lack  of  hysteresis  contradicts  a  wealth  of  data  in  the  lit¬ 
erature.  In  all  of  the  studies  reporting  hysteresis,  the  pressure  membrane 
device  and  axis  translation  technique  were  used.  All  measurements  made  for 
the  present  study  were  at  atmospheric  pressure  in  contrast  to  the  axis  trans¬ 
lational  technique  which  raises  the  sample  ambient  pressure  significantly  in 
order  to  obtain  measurements.  Data  in  McKeen  and  Nielsen  (1978)  show  no 
hysteresis.  Therefore  it  is  assumed  that  the  suction-water  content  relation¬ 
ship  is  linear  and  nonhysteretic. 

Suction  Compressibility 

knowledge  of  suction  changes  must  be  combined  with  a  measure  of  the  soil 
response  to  unit  suction  change  in  order  to  solve  engineering  problems.  As 
previously  discussed,  the  suction  compression  index  (the  ratio  of  volume  change 
to  suction  change)  has  been  selected  for  use  in  this  study.  Some  features  of 
expansive  soil  behavior  must  be  understood  before  proceeding  further.  Most 
expansive  field  soils  exhibit  in  situ  suctions  between  2.5  pF  and  5.0  pF.  As 
noted  in  earlier  work  on  this  project  (McKeen  and  Nielsen,  1978),  soil  volume 
response  to  suction  changes  in  this  range  may  be  modeled  as  a  linear  relation¬ 
ship.  For  purposes  of  engineering  analysis  it  may  be  assumed  that  swel1  ceases 


at  the  capillary  limit  (2.53  pF)  and  at  the  shrinkage  limit  (5.5  to  6.0  pF). 
Several  tests  are  discussed  below  which  provide  a  measure  of  the  slope  of  the 
volume-suction  curve  in  this  range. 


COLE  tests- -The  U.S.  Soil  Conservation  Service  uses  the  COLE  test  to  rate 
soils  as  to  their  shrink-swell  potential.  The  test  was  thoroughly  evaluated 
in  previous  work  (McKeen  and  Nielsen,  1978).  The  COLE  is  defined  as 


COLE  = 


AL 


V 

|  l/ 3 

-  1  = 

+  ll 
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y 

Yw 
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•  1 

1/  3 


-  1 


(4-2) 


AL/L^  =  linear  strain  relative  to  dry  dimensions 
y  ,  =  soil  sample  dry  unit  weight  after  oven  drying 
y  =  soil  sample  dry  unit  weight  at  2.53  pF 
AV/V  ,  =  volume  change  measured  in  the  test 

Several  features  of  the  COLE  test  should  be  noted: 

1.  A  specific  relationship  between  length  change  (AL/Ld)  and  volume 
change  (AV/Vd)  is  explicitly  assumed.  The  actual  length- volume  relationship 
in  the  analysis  problem  must  be  considered  independently. 

2.  The  volume  change  expression  quantifies  the  response  (A V)  in  relation 
to  the  oven-dry  volume  (V^).  Because  real  soils  are  wetter,  usually  between 
3.0  and  4.0  pF,  they  fluctuate  about  a  larger  volume.  Thus,  COLE  is  probably 
a  high  estimate  of  volume  response. 

3.  Standard  COLE  procedures  report  that  volume  change  occurs  between 
2.53  pF  and  oven  dry.  If  the  shrinkage  limit  is  assumed  to  occur  at  5.5  pF, 

7^  is  obtained  by  dividing  COLE  volume  change  by  3.0  pF.  This  assumes  that 
COLE  is  the  change  occurring  over  three  orders  of  magnitude  of  suction  change. 

4.  COLE  procedures  require  coating  the  sample  in  its  natural  state  with 
plastic  resin.  The  sample  is  then  allowed  to  absorb  water  and  swell.  The 
plastic  coating  restricts  this  volume  change.  For  very  wet  samples  little 
restriction  occurs  whereas  considerable  restriction  is  expected  for  dry 
samples.  Previous  work  indicated  that  the  plastic  coating  is  equivalent  to  an 
average  surcharge  pressure  of  10.3  kPa  (1.5  1 b/ in2). 


A  number  of  COLE  tests  were  performed  on  samples  taken  from  each  site 
under  study.  These  data  are  shown  in  Figure  4-25.  In  each  profile  the  soil 
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FIGURE  4-25.  SUCTION  COMPRESSIBILITY  FROM  COLE  TESTS. 


response  varies  dramatically  from  point  to  point.  It  is  clear  that  any  attempt 
to  characterize  these  soils  as  isotropic  homogeneous  materials  will  require 
the  determination  of  macroscopic  soil  parameters  from  a  statistically  sound 
sampling  and  testing  program. 

Clod  tests--The  suction  compression  index  was  also  determined  using  the 
clod  test,  a  modification  of  the  COLE  (McKeen  and  Nielsen,  1978).  A  variety 
of  initial  moisture  conditions  was  used  rather  than  a  single  starting  point 
(2.53  pF  in  the  COLE  test).  This  modification  allows  a  faster  test  and  better 
description  of  the  soil  volume  response.  The  values  of  the  suction  compression 
index  determined  from  the  clod  tests  are  shown  in  Figure  4-26.  The  solid 
lines  are  the  trends  identified  in  the  COLE  data  (Figure  4-25).  The  effect 
of  the  restrictions  involved  in  the  COLE  test  are  evident  in  the  lower  y^ 
values  from  COLE  as  compared  to  those  determined  from  clod  tests.  The  clod 
test  better  depicts  the  unrestrained  response  of  the  soil  to  suction  variation. 
The  COLE  is  perhaps  better  in  depicting  soil  response  under  light  loads;  how¬ 
ever,  the  actual  load  imposed  in  the  test  is  unknown  and  cannot  be  measured. 
Appendix  E  gives  the  detailed  procedures  used  in  the  clod  test. 

Chart  method— The  y^  may  also  be  determined  using  a  chart  method  developed 
from  data  in  the  literature.  Data  published  by  the  Soil  Conservation  Service 
were  used  as  a  source  for  COLE,  percent  clay,  and  cation  exchange  capacity  for 
a  large  number  of  samples  (Soil  Surveys,  1973,  1974,  1976).  Using  a  clay 
mineralogical  classification  system  reported  by  Pearring  (1968)  and  Holt  (1969), 

Holt,  J.  H.,  (1969)  "A  Study  of  Physico-Chemical,  Mineralogical,  and  Engineer¬ 
ing  Index  Properties  of  Fine-Grained  Soils  in  Relation  to  Their  Expan¬ 
sive  Characteristics,"  Dissertation,  Texas  A  &  M  University,  College 
Station,  Texas. 

Pearring,  J.  R.,  (1968)  "A  Study  of  Basic  Mineralogical  Physical-Chemical,  and 
Engineering  Properties  of  Laterite  Soils,"  Dissertation,  Texas  A  &  M 
University,  College  Station. 

"Soil  Survey  Laboratory  Data  and  Descriptions  for  Some  Soils  of  Arizona,"  Soil 
oM't.'ey  investigations  Report  No.  28,  (1974)  Soil  Conservation  Service, 

U.S.  Department  of  Agriculture,  Washington,  DC. 

"Soil  Survey  Laboratory  Data  and  Descriptions  for  Some  Soils  of  California, '• 

Soil  Survey  investigations  Report  No.  24,  (June  1973)  Soil  Conservation 
Service,  U.S.  Department  of  Agriculture,  Washington,  DC. 

"Soil  Survey  Laboratory  Data  and  Description  of  Some  Soils  of  Texas,"  Soil 
Survey  Investigations  Report  No.  30,  (January  1976)  Soil  Conservation 
Service,  U.  S.  Department  of  Agriculture,  Washington,  DC. 
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NMERI  produced  a  chart  which  gives  values  of  without  requiring  suction 
tests.  The  Pearring-Holt  system  classified  soils  into  mineralogical  groups 
according  to  empirical  relationships  between  mineralogy,  plasticity  index  (PI), 
cation  exchange  capacity  (CEC),  and  clay  content  (%C).  CEC  is  measured  in 
milliequivalents  (me).  Pearring  and  Holt  segregated  mineralogical  groups 
according  to  clay  activity  (Ac)  and  cation  exchange  activity  (CEAc) 
where 

rtC  %C 
and 

CEAc  =  (4-3) 

/oU 


(4-3) 


The  NMERI  chart  carried  this  idea  a  step  further  by  correlating  COLE  and 
mineralogical  groups.  About  two  hundred  samples  were  used  to  develop  the 
necessary  regressions.  Data  were  plotted  on  the  CEAc  -  Ac  chart.  Boundaries 
were  then  established  using  mineralogy  percentages,  particularly  the  amount 
of  smectite.  Five  regions  were  established,  ranging  from  pure  smectite  to 
none.  Table  4-6  shows  the  mineralogical  makeup  of  the  regions. 

TABLE  4-6.  REGION  MINERALOGICAL  COMPOSITION 


Region 

Smectite 

11 1 ite 

1 

Kaolinite 

Vermiculite 

I 

>50 

N 

N 

N 

II 

>50 

Tr  -  25 

Tr  -  25 

N 

IIIA 

25  -  50 

10  -  25 

5  -  10 

N 

1 1  IB 

5  -  50 

5-25 

Tr  -  25 

N 

IVA 

Tr  -  10 

5  -  25 

5  -  50 

N 

IVB 

Tr 

10  -  25 

5  -  50 

N 

VA 

N 

Tr  -  25 

5-50 

Tr  -  25 

VB 

N 

N 

10  -  25 

Tr 

All  values  expressed  as  percent  of  the  clay  fraction, 
N  =  none 

Tr  =  trace,  <5  percent 
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With  the  regions  identified,  the  y^  for  soils  in  each  region  were  studied. 
Linear  regressions  were  calculated  for  yh  and  percent  clay  in  each  region. 
All  coefficients  of  determination  were  above  0.9.  The  y^  values  for  a  pure 
clay  are  shown  on  the  chart  Figure  4-27.  The  five  regions  are  indicated  by 
the  Roman  numerals  given  in  Table  4-7.  To  obtain  a  y^  value  for  a  soil, 
the  chart  value  is  determined  by  plotting  CEAc  versus  Ac  and  multiplying  the 
number  obtained  by  the  clay  content  as  a  decimal. 

Example:  %C  =52  Ac  =  0.981 

PI  =  51  CEAc  =  0.610 

CEC  =  31.7  me/100  g  yh  =  0.52  (0.163)  =  0.085 


SUCTION  CONCLUSIONS 


The  relationship  of  suction  to  water  content  depicts  the  most  important 
feature  of  a  clay:  its  ability  to  attract  and  hold  water  on  the  clay  particle 
surfaces.  It  is  this  struggle  for  water  between  the  soil  and  its  surroundings 
that  produces  moisture  movement  which  invariably  is  accompanied  by  volume 
change.  As  the  suction  goes  below  about  2.5  pF  the  principles  of  saturated 


ACTIVITY,  Ac 


FIGURE  4-27.  CHART  FOR  SUCTION  COMPRESSIBILITY  PREDICTION. 
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soil  mechanics  become  increasingly  accurate  in  modeling  the  soil  behavior. 
Above  about  5.5  pF  the  cohesion  is  so  great  that  volume  changes  are  of  little 
consequence.  The  activity  of  importance  in  this  study  is  between  2.5  and 
5.5  pF,  the  range  of  activity  of  clay  soils.  The  more  water  active  in  this 
range,  the  more  potential  for  movement. 

Both  the  TCP  and  FP  methods  of  suction  measurement  can  be  used  to  measure 
soil  suction.  In  the  range  of  suction  encountered  in  this  study,  laboratory 
measurements  were  good  estimators  of  field  suction.  The  FP  has  the  advantage 
of  a  wider  range  than  the  TCP  method.  It  is  generally  assumed  that  the  TCP 
equilibrates  faster  than  the  FP;  however,  this  study  indicates  that  this  is 
not  true  for  ceramic-covered  TCPs.  Suction  determinations  should  be  made 
routinely  on  clay  soils  for  any  pavement  project. 

The  Yh  varies  dramatically  within  expansive  soils.  The  variation  is  a 
natural  result  of  the  heterogeneous  nature  of  the  deposits.  The  characteriza¬ 
tion  of  such  profiles  must  involve  a  large  number  of  measurements  and  the 
determination  of  properties  that  will  model  the  macroscopic  response.  The  use 
of  a  few  widely  spaced  samples  is  clearly  a  dangerous  approach  to  soil  charac¬ 
terization  in  expansive  clays. 

The  determined  from  clod  tests  is  a  good  measure  of  the  unrestrained 
response  of  a  soil  to  suction  variation.  Clod  tests  yield  a  better  estimate 
of  the  soil  behavior  than  the  COLE  test.  In  expansive  soil  profiles  with  wet 
densities  of  1442  to  1762  kg/m3  (90  to  110  lb/ft3),  the  COLE  models  overburden 
loads  to  0.61  to  0.73  m  (2.0  to  2.4  ft)  of  overburden.  Because  swell  reduc¬ 
tion  by  loads  is  linear  in  the  early  stages,  it  is  reasonable  to  assume  that 
the  COLE  test  can  be  used  to  characterize  the  top  1.5  to  1.8  m  (5.0  to  6.0  ft) 
of  a  soil  profile  without  making  overburden  load  corrections  to  the  swell 
predictions. 

Clays  which  are  cracked  extensively  are  expected  to  swell  three- 
dimensionally,  approaching  the  COLE  assumption.  When  the  clay  is  wetter  and 
does  not  contain  cracks,  it  becomes  confined  by  lateral  soil  pressure.  Then 
virtually  all  volume  change  occurs  in  a  vertical  direction.  Because  of  this 
single  degree  of  freedom,  internal  shear  stresses  must  increase  significantly, 
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producing  the  slickenside  feature  found  in  active  clays.  The  creation  of 
these  shear  failures  reduces  the  vertical  swell  below  the  total  volume  change 
that  would  have  occurred  in  an  unrestrained  soil.  This  is  an  important  consi 
deration  and  virtually  no  information  on  it  exists.  Further  consideration  is 
given  to  this  aspect  in  later  parts  of  this  report. 

The  presence  of  discontinuities  in  clay  soil  profiles  has  a  profound 
influence  on  the  in  situ  behavior.  In  the  present  study  it  was  apparent  that 
the  interface  between  two  clay  layers  may  provide  an  avenue  for  water  move¬ 
ment.  Because  of  the  very  low  permeability  of  intact  clay,  virtually  all 
significant  variation  in  moisture  occurs  by  means  of  seams  of  more  permeable 
material  embedded  in  the  clay  or  through  fissures  resulting  from  clay  volume 
changes.  The  determination  of  the  patterns  of  such  features  is  essential  to 
soil  characterization. 
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V.  ELEVATION  STUDIES 


ELEVATION  DATA  GATHERING 

Soil  and  pavement  profile  characteristics  were  periodically  observed  by 
conventional  rod  and  level  surveying.  Elevations  were  surveyed  along  three 
lines:  Line  1,  uncovered  soil;  Line  2,  pavement  edge;  Line  3,  pavement 
interior.  Elevations  were  recorded  at  0.6-m  (2-ft)  intervals,  in  most  cases, 
over  a  horizontal  distance  of  40.6  m  (130  ft).  Several  observations  along 
these  lines  were  made  at  0.2-m  (0.5-ft)  intervals  to  study  the  suitability 
of  the  data  for  analysis  using  transform  techniques.  In  addition  to  these 
observations,  some  measurements  on  rough  pavements  were  made  in  June  1980  to 
provide  a  data  base  for  the  analysis.  These  measurements  are  described  in 
Table  5-1.  At  each  of  the  regular  sites,  a  suitable  bench  mark  was  identified 
or  constructed  to  serve  as  an  elevation  reference.  Site  layouts  are  illustra¬ 
ted  in  Figures  5-1  to  5-3. 


TABLE  5-1.  ROUGH  PAVEMENTS  SURVEYED 


Location 

Designation 

Feature  surveyed 

Gal lup  Airport 

GAL  4,  5,  6 

Runway  centerline 

Jackson  Airport 

JSN  5,  6,  7,  8,  9 

Taxiway 

Jackson  Airport 

JSN  11,  12 

Gate  apron 

Jackson  Airport 

JSN  4,  15,  16 

General  aviation  apron 

Cl  inton,  Mississippi 

JSN  10 

Highway  pavement 

Dallas/Fort  Worth  Airport 

DFW  7  to  12 

Access  roadway 

Dallas  Naval  Air  Station 

NAS  1,  2 

Apron/taxiway 

Dallas  Naval  Air  Station 

NAS  3,  4 

Taxi  way 

Greater  Fort  Worth 
(Amon  Carter  Field) 

Visual  survey  and 
inspection  of  all 
pavements 
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TAXIWAY  E31  PAVED  SHOULDER 


38 

cm 

(15  in)  PORTLAND  CEMENT  CONCRETE 

10  cm 

(4  ini 

ASPHALT  CONCRETE 

23 

cm 

(9  in)  CEMENT  STAB  BASE 

15  cm 

(6  in) 

LIME  STAB  SOIL 

23 

cm 

(9  in)  LIME  STAB  SOIL 

ACCESS  ROAD 


4  cm  (1.5  in)  ASPHALT  CONCRETE 
11  cm  (4.5  in)  ASPHALT  CONCRETE  DENSE  GRADED 
15  cm  (6  in)  LIME  STAB  AT  7  PERCENT 
30  cm  (12  in)  LIME  STAB  AT  4  PERCENT 


O  BORINGS 

O  INSTRUMENTATION  SITE 


o  4  (OCT)  78) 


FIGURE  5-1.  DALLAS/FORT  WORTH  REGIONAL  AIRPORT  SITE  LAYOUT. 
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TAXI1JAY 


SHOULDER 


8  cm  (3  in)  BITUMINOUS  SURFACE 
28  cm  11  in)  CEMENT  TREATED  BASE 
15  cm  (6  in)  LIME  TREATED  SUBGRADE 
22  cm  (9  in)  COMPACTED  SUBGRADE 


DOUBLE  SURFACE  TREATMENT 
WITH  COVER  AGGREGATE 
13  cm  (5  in)  SELECT  BORROW 
(CRUSHER  FINES) 


FIGURE  5-2.  GALLUP  MUNICIPAL  AIRPORT  SITE  LAYOUT. 
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AD-A104  660 
UNCLASSIFIED 


NEW  MEXICO  ENGINEERING  RESEARCH  INST  ALBUQUERQUE  F/G  8. 

DESIGN  OF  AIRPORT  PAVEMENTS  FOR  EXPANSIVE  SOILS. (U) 

JAN  81  R  6  MCKEEN  F29601-76-C-0015 

NMERI-AP-37  FAA/RD-Rl/25  NL 


FAA/RD-Rl/25 


EAST  PARALLEL  TAXIWAY 


EAST  PARALLEL  TAXIWAY 

2.4  cm  \ 6  in)  ASPHALT  CONCRETE  OVERLAY 
29.2  cm  (11.5  in)  PORTLAND  CEMENT  CONCRETE  PAVEMENT 
13  cm  (5  in)  SOIL  CEMENT 


TURNOUT 


FIGURE  5-3.  JACKSON  MUNICIPAL  AIRPORT  SITE  LAYOUT. 
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GENERAL  AVIATION  PARKING  APRON 


SOIL  SURFACE  ELEVATIONS 


Undisturbed  soil  profiles  rise  and  fall  with  wetting  and  drying  influ¬ 
ences  imposed  by  the  environment.  It  was  recognized  early  in  the  literature 
survey  that  understanding  this  in  situ  behavior  was  important.  Therefore  the 
natural  soil  surfaces  adjacent  to  the  pavements  were  observed.  Figures  5-4 
to  5-6  show  typical  data  from  the  soil  surface  at  each  site.  The  GAL  site  is 
virtually  flat;  drainage  is  across  Line  1  shown  in  Figure  5-2.  At  the  DFW 
site  the  elevations  are  flat  from  0  to  18  m  (0  to  60  ft),  then  slope  toward  a 
culvert  at  about  a  1-percent  slope.  Drainage  is  perpendicular  to  Line  1  shown 
in  Figure  5-1  and  toward  the  higher  station  numbers.  The  JSN  site  slopes  for 
its  entire  length  at  about  0.5  percent.  These  areas  were  not  prepared  in  any 
way. 

Variations  of  the  surface  elevation  were  compared  in  several  ways.  In 
one  method,  the  maximum  elevation  difference  was  determined  at  each  point 
along  a  line.  Results  for  DFW  and  GAL  are  illustrated  in  Figure  5-7.  The  DFW 
results  reflect  variation  in  the  range  of  swell  at  regular  intervals  of  between 
3  and  6  m  (10  and  26  ft).  The  dominant  spacing  is  5.5  m  (18  ft).  These  obser¬ 
vations  are  consistent  with  the  characteristics  of  gilgai  surfaces  on  clay 
soils.  Data  for  the  GAL  site  show  about  half  as  much  vertical  movement  and  no 
distinct  pattern  of  response  variation. 

Another  method  of  investigating  the  variation  of  soil  response  to  envi¬ 
ronmental  changes  was  inspection  of  individual  profile  plots  for  each  observa¬ 
tion  date.  The  GAL  data  for  June  1980  clearly  indicate  the  varying  response 
of  the  soil  along  the  line  (Figure  5-8).  The  elevations  reflect  differential 
soil  response  exhibiting  a  regular  mound-depression  pattern  consistent  with 
the  gilgai  model.  The  spacing  of  the  mound-depression  pattern  is  nearly  6  m 
(20  ft).  Inspection  of  Figure  5-4  reveals  no  such  distinct  pattern  for  March 
1979  (drying  soil).  This  discrepancy  leads  to  the  conclusion  that  soil  rate 
of  drying  varies,  producing  depressions  where  water  moves  out  and  evaporates 
faster.  Likewise,  mounds  appear  where  water  is  moving  out  of  the  soil  more 
slowly.  Thus  the  development  of  differential  surface  features,  as  shown  in 
Figure  5-8,  is  attributed  to  varying  field  soil  permeabi 1 i ty . 

A  third  representation  of  the  soil  surface  elevation  was  obtained  with  a 
fast  Fourier  transform  to  determine  the  amplitudes  of  discrete  frequency  com¬ 
ponents  in  the  profile.  The  procedure  used  is  explained  in  the  next  section. 
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FIGURE  5-4.  ELEVATIONS  ON  SOIL  SURFACE  AT  GAL. 
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FIGURE  5-5.  ELEVATIONS  ON  SOIL  SURFACE  AT  DFW. 
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FIGURE  5-6.  ELEVATIONS  ON  SOIL  SURFACE 
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FIGURE  5-7.  RANGE  OF  ELEVATION  CHANGES  AT  DFW  AMD  GAL. 
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FIGURE  5-3.  EVIDENCE  OF  SOIL  DIFFERENTIAL  RESPONSE  AT  GAL. 


The  fundamental  idea  presented  in  this  report  is  that  soil  volume  responds 
to  soil  moisture  suction  changes.  If  this  is  true,  the  field  elevations  should 
reflect  this  relationship.  Figure  5-9  is  a  plot  of  change  in  elevation  versus 
change  in  suction  at  the  DFW  and  GAL  field  sites.  Thirteen  of  the  observed 
points  were  used  to  calculate  a  statistical  relationship  between  elevation 
(measured  at  several  points  near  the  suction  measurement)  and  the  average  suc¬ 
tion  in  the  top  1.5  m  (5  ft)  of  the  soil.  Results  yielded 

AH  =  0.0164  -  0.2202  (Ah  )  (5-1) 

a 

where 

AH  =  elevation  change,  ft 

Ail  =  average  suction  in  the  top  5  ft  (1.5  m),  pF 

a 

number  of  samples  =  13 

coefficient  of  determination  (R2)  =  0.922 

sample  standard  deviation  =  0.032  ft 

The  two  deviant  points  shown  in  Figure  5-9  indicate  that  the  soil  surface 
is  subsiding  as  the  soil  is  wetting.  This  behavior  is  contrary  to  the  concept 
of  the  suction  compression  index  previously  presented.  The  complete  model 
recommended,  however,  requires  evaluation  of  the  lateral  restraint  on  the  soil 
in  the  profile.  This  idea  is  illustrated  in  Figure  5-10.  Figure  5-10a  shows 
an  uncracked,  wet  soil  profile.  When  the  soil  dries,  it  cracks  as  shown  in 
5-10b  at  an  interval  H  ,  to  depth  1^.  Crack  width  at  the  surface  is  W£;  and, 
assuming  a  triangular  shape,  the  average  crack  width  is  W^/2 .  It  is  proposed 
in  5-10c  that,  as  the  cracks  fill  with  surface  water,  a  process  of  horizontal 
swelling  fills  the  cracks  in  the  soil.  This  stress  relaxation  at  the  cracks 
also  results  in  surface  subsidence.  Assuming  that  the  volume  of  surface  sub¬ 
sidence  is  equal  to  the  volume  of  cracks  filled,  estimates  of  crack  spacing 
and  depth  can  be  made.  The  problem  is  illustrated  in  Figure  5-11.  From  the 
information  shown, 

Vj  =  (H  -  W  )  AH  per  unit  of  crack  length  (5-2) 

Wc 

Vz  =  ■£-  (Zc  -  AH)  per  unit  of  crack  length  (5-3) 

where 

AH  =  surface  subsidence 
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O  DFW 

A  A  GAL  NO. 

0.04-  O  GAL  NO. 


FIGURE  5-9.  CORRELATION  OF  ELEVATION  AND  SUCTION  CHANGES 


c.  WETTING  SOIL 

Hc  =  CRACK  SPACING  z  =  ORIGINAL  DEPTH 

Zc  =  CRACK  DEP™  l  =  DEPTH  AFTER  DRYING 

Uc  =  !!?5 TuUr1  CRACK  Z  =  DEPTH  AFTER  WETTING 

W I  u  1  H  W 

FIGURE  5-10.  SOIL  SUBSIDENCE  CAUSED  BY  WETTING. 
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Zc  =  CRACK  DEPTH 

Wc  =  MAXIMUM  CRACK  WIDTH 

AH  =  SURFACE  SUBSIDENCE 

Vi  =  VOLUME  OF  SOIL  SUBSIDENCE 

V  =  VOLUME  OF  SOIL  TO  FILL  CRACKS 

FIGURE  5-11.  CRACK  SUBSIDENCE  PROBLEM. 


Equating  these  volumes  and  solving  for  Zc  yields 


Zc  -  AH 


f2Hc 


(5-4) 


A  comparison  with  observations  at  the  sites  studied  is  provided  in 

Table  5-2.  The  results  indicate  that  the  calculated  Z  values  are  about  twice 

c 

those  measured.  These  are  interpreted  as  good  estimates  because  the  average 
crack  depths  were  measured  with  a  steel  tape  from  the  surface.  Obviously  the 
measurement  could  only  be  made  to  a  depth  at  which  the  crack  would  not  permit 
further  penetration  of  the  tape.  The  depths  measured  could  easily  be  about 
half  the  actual  distance. 


These  results  show  that  a  cracked  expansive  soil  profile  will  undergo  sur¬ 
face  subsidence  upon  wetting.  When  the  surface  cracks  are  filled  by  slumping, 
swelling  soil,  then  surface  heaving  occurs.  The  nature  of  cracking  in  an 
expansive  soil  profile  very  strongly  influences  its  behavior.  A  method  of 
accounting  for  this  influence  is  needed. 
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TABLE  5-2.  CRACKING  CHARACTERISTICS  COMPARED 


Observed  values 

Calculated 

w  ,  m 

(ft) 

— 

Zc’ 

(ft) 

Zc,  m 
(ft) 

0.9 

0.023 

Pi 

0.46 

0.79 

Sep  79 

(3.0) 

(0.08) 

IHI 

0.5) 

(2.6) 

GAL 

0.2 

0.003 

Nov  79 

(0.8) 

(0.01) 

HI 

(0.13) 

(0.23) 

Notes:  Hc  =  interval  between  cracks;  W  =  crack  width  at 
c  surface;  Z  =  depth  of  cra£k;  .1H  =  surface 
subsi dence; 


PAVEMENT  SURFACE  DATA  ANALYSIS 

Root-Mean-Square  Profile  Elevation 

Observations  of  pavement  surface  elevations  were  made  several  times  during 
the  field  experiment  and  coincided  with  determinations  of  soil  moisture  and 
suction.  A  variety  of  tools  have  been  developed  for  analyzing  roughness.  In 
the  present  work,  several  specific  methods  were  used.  The  RMS  (root-mean- 
square)  profile  elevation  and  ARMS  (root-mean-square  elevation  changes)  were 
computed  in  accordance  with  Air  Force  roughness  procedures  (Seeman  and 
Nielsen,  1976). 

The  following  group  of  equations  describe  profile  characteristics  to  be 
used  and  their  interpretation. 


N 


where 


y  =  average  elevation 
y.j  =  individual  elevation  observation 
N  =  total  number  of  observations 


Seeman,  Douglas  R.,  and  Nielsen,  John  P.,  (1976)  "Suggested  Procedure  to  Iso¬ 
late  and  Repair  Runway  Roughness,"  Letter  Report  C ERF- 1 1 ,  University  of 
New  Mexico  Civil  Engineering  Research  Facility,  Albuquerque. 
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(5-6) 


N-l 

ARMS  =  — 


(yi+l  "  yi'2 
'T^~2 


where 

ARMS  =  root-mean-square  elevation  change  from  point  to  point 
y.+1,  y-  =  consecutive  elevation  observations 
N  =  total  number  of  observations 

If  the  ARMS  is  divided  by  the  observation  spacing  of  0.61  m  (2.0  ft),  it 
then  represents  the  RMS  slope  of  the  pavement  profile. 

Transform  Method 

Fast  Fourier  transforms  are  used  widely  in  the  analysis  of  signals.  Their 
value  is  derived  from  the  capability  of  obtaining  the  frequency  content  of  a 
time-domain  signal.  Knowing  the  frequency  content,  the  electrical  engineer  is 
in  a  much  better  position  for  solving  signal  analysis  problems.  In  the  present 
study,  it  has  been  proposed  that  expansive  soil  profiles  exhibit  microrelief 
of  a  periodic  nature.  It  has  also  been  stated  that  roughness  of  pavement 
structures  founded  on  expansive  soils  is  related  to  the  soil  microrelief.  If 
these  assertions  are  true,  a  frequency  domain  plot  of  the  surface  elevation 
data  should  exhibit  the  dominant  frequency  components  characteristic  of  the 
soil -pavement  system.  Thus  the  fast  Fourier  transformation  appears  to  be  a 
viable  analysis  technique. 

Pavement  surface  analysis  involved  calculating,  from  the  space  domain 
data,  a  discrete  frequency  domain  function  using  a  fast  Fourier  transform. 
Literally  this  transform  expresses  the  data  in  terms  of  a  complex  Fourier 
series  in  the  space  domain.  The  Fourier  coefficients  are  calculated  using  a 
least  squares  technique.  In  the  frequency  domain,  half  of  the  amplitude  com¬ 
ponents  are  plotted  versus  frequency  for  use  in  the  analysis.  Each  component 
also  has  a  phase  angle  associated  with  it. 

The  first  step  in  the  procedure  is  to  obtain  a  space  domain  function 
representing  the  observed  data  set  as  a  complex  Fourier  series.  The  observed 
data  can  be  fitted  exactly  using  a  Fourier  series. 
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y(xn )  =1^+2 

n  n=0 


N-l 


_ 27rmn  A  u  2Tmin 

a„COS  - TT—  +  b  Sin  - n- 

n  N  n  N 


(5-7) 


where 


y(x  )  =  the  observed  elevation  at  x„ 
n  n 

x  =  horizontal  station  of  elevation  measurement 
n 

m,  n  =  0,  1 ,  2. . .N 

N  =  number  of  observations 

a  ,  b  =  Fourier  coefficients 
n  n 

Equivalent  forms  for  this  equation  are  given  by  Stearns  (1975): 


Ac 


N-l 


y(xn)  =  r  +  n:Q  Ancos  (nfoXn  +  an] 


(5-8) 


+N/2 

y(x  )  =  E  c  exp(jnf0x  ) 
n  n=  -N/2  n  n 


(5-9) 


Frequency  components  of  the  function  are  expressed  using  a  discrete 
Fourier  transform  (DFT): 


N-l 

F(jw)  =  E  y(x  )e 
n=0  n 


■j2mn/N 


(5-10) 


N-l 

s  y(x  )  cos 
n=0  n 


(2mn ) 


N-l 

-  j  E  y(x  )  sin 
n=0 


2Trmn 


(5-11) 


where 

F(jw)  =  complex  Fourier  transform  of  y(*n) 
y(xn)  =  observed  elevation  at  x 
m,  n  =  same  as  above 

j  =  /-“T 

Some  features  of  this  technique  are  illustrated  in  the  following  example. 
Figure  5- 1 2a  is  a  space  domain  plot  of  64  elevation  observations  taken  on  the 
surface  of  an  aircraft  parking  apron  at  the  JSN  site.  The  data  are  properly 


Stearns,  Samuel  D. ,  (1975)  Digital  Signal  Analysis,  Hayden  Book  Company, 
Rochelle  Park,  N.J. 
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HORIZONTAL  DISTANCE 


FIGURE  5-12.  ILLUSTRATION  OF  FOURIER  TRANSFORMATION  OF  A  PROFILE 


represented  here  as  discrete  points  rather  than  a  continuous  function.  In 
Figure  5- 1 2b  the  amplitude  versus  frequency  data  calculated  using  the  DFT  are 
shown.  Because  the  Fourier  series  is  complex,  each  frequency  component  also 
has  a  phase  angle  associated  with  it  which  is  not  shown.  Notice  that  the  data 
in  Figure  5- 1 2a  appear  to  be  periodic  with  a  half-amplitude  of  0.09  m  (0.3  ft) 
and  a  wavelength  (or  period)  of  about  30.5  m  (100  ft).  In  the  frequency  domain, 
the  surface  is  expressed  at  discrete  frequency  components.  The  frequencies 
are  determined  by  the  sampling  interval  and  overall  length.  Thus, 

f  =  spatial  frequency  of  the  nth  component,  cycles 
n  per  unit  length 

n  =  1,  2,  3  ...  N 

L  =  overall  sampling  length 

A  =  wavelength  of  the  nth  component 

Therefore,  the  frequencies  with  estimates  are 
fi  =  ]28^ft  =  0.0078  ft'1  (0.26  m"1) 
f2  --  =  01^6  ft"1  (0.0513  m"1) 


fja  =  y^Yt  =  °*250  ft_1  <°-820  m_1) 

These  have  corresponding  wavelengths,  A  ,  as  follows: 

Ax  =  i-  =  128  ft  (39.0  m) 

*  i 

A 7  =  -l-  =  64.0  ft  (19.5  m) 

*  2 

•  •  ■ 

•  •  • 

A37  =  r-L  =  4.0  ft  (1.2  m) 

T  32 
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Only  half  of  the  computed  frequency  components  are  unique.  Thus,  for 
64  data  points  in  the  space  domain,  only  32  points  in  the  frequency  domain 
are  determined.  The  data  folds  at  the  midpoint  (called  the  foldovar  or  Nyquist 
frequency)  and  subsequent  points  are  a  mirror  image  of  the  first  half.  The 
sampling  procedure  used  in  the  present  study  yielded  frequency  components  from 
0.026  m"1  (0.0078  ft-1)  to  0.82  m"1  (0.25  ft'1}.  The  computed  frequency  com¬ 
ponents  are  at  discrete  frequencies  determined  by  the  sampling  procedure. 
Because  the  real  profile  contains  components  at  different  frequencies  from  the 
calculated  components,  the  real  frequency  components  are  distributed  to  those 
calculated.  In  the  above  example  (Figure  5-12a),  the  wavelength  was  about 
30.5  m  (100  ft).  The  frequency  domain  represents  this  as  two  components  at 
wavelengths  of  39.0  m  (128  ft)  and  19.5  m  (64  ft).  Note  that  the  sum  of  the 
amplitudes  of  these  components  is  about  equal  to  the  space  domain  half¬ 
amplitude;  i.e.,  0.055  m  +  0.043  m  =  0.097  m  (0.18  ft  +  0.16  ft  =  0.34  ft). 

The  overestimate  of  amplitude  results  from  phase  differences  in  the  components. 

Appendix  B  contains  frequency  domain  data  for  the  pavements  studied. 
Regression  lines  were  fitted  through  the  data.  Equations  are  of  the  following 
form : 


an  =  c(f)b  (5-13) 

an  =  half-amplitude,  ft 

f  =  spatial  frequency,  ft"1 
c,  b  =  regression  constants 

Velasco  (1980)  has  computed  statistical  relationships  between  these 
regression  coefficients  and  a  number  of  factors  involved  in  the  swelling  soil 
problem  for  highway  pavements.  Table  5-3  shows  the  factors  considered  in  the 
regression.  The  results  of  his  study  yielded  regression  equations  for  the 
prediction  of  c  and  b  for  a  particular  profile.  These  could  then  be  related 
to  the  Present  Serviceability  Index  of  a  highway  pavement  surface  built  on 
the  soil . 


Velasco,  Manuel  0.,  (1980)  "Pavement  Roughness  on  Expansive  Clays,"  Thesis  at 
Texas  A  &  M  University,  College  Station,  Texas. 
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TABLE  5-3.  VELASCO'S  INDEPENDENT  VARIABLES 


OEPTH  =  effective  depth  of  pavement,  inches 


TH  =  mean  value  of  Thornthwaite  moisture  index  for  a  20-year 
period 

RANGE  =  range  of  values  of  Thornthwaite  moisture  index  for  a  20- 
year  period 

TIME  =  time  since  construction  or  last  rehabilitation  before  the 
roughness  was  measured,  years 


CLAY  =  percent  clay  (grain  size  less  than  0.002  mm) 
AC  -  activity  j?Jaiticit^Jndexj 


CEC  =  cation  exchange  capacity 
CEAC  =  cation  exchange  activity 


me 

’  100  g 

CECl 

CLAY 


COLE  =  coefficient  of  linear  extensibility 


ESP  =  exchange  sodium  percentage 


Velasco,  1980 
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Other  Methods 

Stone  and  Dugundji  (1965)  studied  microrelief  of  a  variety  of  surfaces 
for  the  purpose  of  developing  a  quantitative  method  of  expressing  the  surface 
characteristics.  Their  approach  involved  using  a  cosine  series  representation 
of  the  profile  and  determining  a  number  of  quantitative  characteristics  of  the 
profile.  A  group  of  similar  characteristics  were  obtained  in  the  present  study 
as  shown  in  Table  5-4.  Data  computed  for  the  pavements  in  this  study  are  shown 
in  Appendix  C. 

An  important  pavement  characteristic  for  the  study  of  roughness  is  the 
stiffness  of  the  pavement  structure  itself.  The  approach  taken  by  Velasco 
(1980)  used  an  effective  depth  of  the  pavement  structure  to  characterize  it 
for  study.  The  same  approach  was  followed  to  characterize  the  pavements 
studied  in  the  present  work. 


Stone,  R.  0.,  and  Dugundji,  T. ,  (1965)  "A  Study  of  Microrel i ef —  Its  Mapping, 
Classification,  and  Quantification  by  Means  of  a  Courier  Analysis," 
iu.jincev'rii  'icolojy,  Vol.  1,  No.  2:39-187. 


Symbol 


M 


P 


K 


A 
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VI.  AIRCRAFT  SIMULATIONS 


TAXI  COMPUTER  CODE 

The  TAXI  computer  code  was  developed  to  determine  the  dynamic  response 
of  an  aircraft  to  runway  roughness  during  takeoff  and  constant-speed  taxi. 

The  mathematical  model  has  been  programmed  for  a  CDC  6600  digital  computer, 
and  results  to  date  indicate  good  agreement  with  results  of  runs  using  instru¬ 
mented  aircraft.  The  procedures  involved  are  fully  described  by  Gerardi  and 
Lohwasser  (1973  and  1974).  A  summary  of  this  description  is  provided  here 
from  Gerardi  (1978). 

The  general  mathematical  model  is  represented  as  an  asym¬ 
metrical  body  with  a  nose  and  a  right  and  a  left  main  landing 
gear.  Each  landing-gear  strut  is  assumed  to  have  point  contact 
with  the  profile,  and  it  is  assumed  that  each  landing  gear 
traverses  a  different  profile.  The  model  has  aerodynamic  lift 
and  drag  and  thrust  applied  at  the  center  of  gravity  of  the 
aircraft. 

The  aircraft  is  free  to  roll,  pitch,  plunge,  and  translate 
horizontally  down  the  runway,  and  each  landing-gear  unsprung 
mass  is  free  to  translate  vertically.  In  addition  to  these 
rigid-body  degrees  of  freedom,  there  are  up  to  30  flexible 
modes  of  vibration.  This  aircraft  motion  is  controlled  by  the 
landing-gear  strut  forces  and  the  lift,  drag,  thrust,  and 
resisting  parameters  of  aircraft  mass  and  inertia. 

The  landing-gear  struts  are  nonlinear,  single-  or  double¬ 
acting,  oleopneumatic  energy-absorbing  devices  and  are  repre¬ 
sented  in  the  model  as  the  sum  of  the  three  forces--pneumatic, 
hydraulic,  and  strut  bearing  friction.... 

The  runway  elevation  data  are  input  into  the  model  in 
0.61-m  (2-ft)  long  increments.  The  profile  is  made  continuous 
by  fitting  the  following  polynomial  through  three  points  and 
the  slope  at  the  end  of  the  previous  profile  segment: 

y(x)  =  a!  +  a2x  +  a3xz  +  aux3  (6-1 ) 


Gerardi,  Anthony  G. ,  (1978)  "Dynamic  Response  of  Aircraft  to  Pavement 

Uneveness,"  in  Kuccaivk  in  Airport  !  u:\vicn  tr.,  Special  Report  175,  Trans¬ 
portation  Research  Board,  Washington,  DC,  pp.  91-96. 

Gerardi,  Anthony  G.  and  Lohwasser,  Adolph  K. ,  (1973)  Am ;:<u:r  iXunK  Hu 
Prediction  of  Aircraft  Response  to  Runway  Rou.jhr.css,  Volume  i,  :  rojrari 
1/coclopmcnt,  U.  S.  Air  Force  Flight  Dynamics  Laboratory,  Technical  Report 
No.  AFWL-TR-73-109,  Vol.  I,  Wri ght-Patterson  AFB,  Ohio. 

Gerardi,  Anthony  G.  and  Lohwasser,  Adolph  K. ,  (1974)  I'onpat.  >■  !ve.:r,v-;  for  ; 
Prediction  of  Aircraft  Hoot  otter  to  Runwuy  Roujnn  roc,  'r 

Manual,  U.S.  Air  Force  Flight  Dynamics  Laboratory,  Technical  Report  No. 
AFWL-TR-73-109,  Vol.  II,  Wright-Patterson  AFB,  Ohio. 
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where  ax,  a2,  a3,  and  a*,  =  coefficients  derived  from  elevation 

and  slope  data.  This  is  done  for  each  of  the  three  lines  of 

runway-profile  data. 

Several  factors  involved  in  the  present  study  affect  the  use  of  this  simu¬ 
lation.  The  data  were  gathered  on  the  basis  of  anticipated  roughness  patterns 
caused  by  swelling  soil.  Wavelengths  expected  were  in  the  range  3.0  to  15.0  m 
(10  to  50  ft).  The  longest  wavelength  obtained  was  the  sampling  interval  (the 
total  length  sampled):  39.0  m  (128  ft). 

Previous  studies  investigated  the  effect  of  filtering  out  long  wavelength 
components  (Seeman  and  Nielsen,  1976).  These  results  indicate  that  at  the 
39.0  m  (128  ft)  cut-off  the  accelerations  are  about  96  percent  of  those  for 
the  unfiltered  data  for  the  B727-100  aircraft.  It  was  concluded  that  the  pre¬ 
dicted  accelerations  were  satisfactory  for  roughness  evaluation  in  this  study. 
Results  reported  are  averaged  over  four  runs  on  the  same  pavement  segment, 
where  the  segments  were  placed  head  to  toe. 

A  single  profile  was  available  for  each  observation  made.  Thus,  all  of 
the  landing  gear  was  assumed  to  traverse  the  same  profile. 

Another  factor  not  simulated  is  the  flexural  movement  of  the  pavement 
structure.  The  TAXI  code  considers  the  surface  of  the  pavement  to  be  rigid; 
the  pavement  is  not  allowed  to  deflect  under  aircraft  loads.  Some  of  the 
roughness  encountered  on  airport  pavements  results  from  the  elastic  deflection 
of  the  pavement  under  aircraft  loads.  These  deflections  are  particularly 
important  for  large  aircraft  and  light  pavements.  In  this  study  all  pavements 
are  treated  as  rigid  structures  in  the  aircraft  simulations. 

Four  aircraft  were  used  in  the  simulation:  the  B727-100,  B727-200, 
DC-9-40,  and  the  TriJet  Composite.  Input  aircraft  characteristics  and  initial 
conditions  are  in  Appendix  D.  All  aircraft  simulations  consisted  of  constant- 
speed  taxi  runs  over  the  pavement  profiles  under  study. 

RESULTS  OF  SIMULATIONS 

Results  of  aircraft  simulations  are  presented  in  Table  6-1.  The  results 
shown  are  RMS  and  peak  accelerations  at  the  pilot  station.  The  date  column  in 


TABLE  6-la.  RESULTS  OF  AIRCRAFT  SIMULATIONS 


Line 

Da  te 

Ai rcraft 

Pilot  station  acceleration,  g 

RMS 

Peak 

Gal  2 

Mar  79 

727-100 

0.0126 

0.046 

Mar  79 

DC-9-40 

0.0136 

0.049 

Mar  79 

727-200 

0.0103 

0.042 

Mar  79 

Wide  body 

0.0151 

0.063 

Jun  80 

727-100 

0.0111 

0.046 

Gal  3 

Mar  79 

727-100 

0.0219 

0.059 

Mar  79 

DC-9-40 

0.0136 

0.045 

Mar  79 

727-200 

0.0080 

0.026 

Mar  79 

Wide  body 

0.0203 

0.056 

Jun  80 

727-100 

0.0174 

0.047 

Apr  79 

727-100 

0.0147 

0.049 

May  79 

727-100 

0.0164 

0.050 

Jul  79 

727-100 

0.0161 

0.057 

Sep  79 

727-100 

0.0157 

0.052 

Nov  79 

727-100 

0.0149 

0.054 

Mar  80 

727-100 

0.0162 

0.057 

Jun  80 

727-100 

0.0161 

0.058 

Apr  79 

DC-9-40 

0.0144 

0.053 

Apr  79 

727-200 

0.0198 

0.074 

Apr  79 

Wide  body 

0.0213 

0.081 

DFW  3 

Apr  79 

727-100 

0.0113 

0.035 

Apr  79 

DC-9-40 

0.0082 

0.030 

Apr  79 

727-200 

0.0134 

0.040 

Apr  79 

Wide  body 

0.0143 

0.049 

Apr  79 

727-100 

0.1270 

0.387 

May  79 

727-100 

0.1255 

0.386 

Jul  79 

727-100 

0.1263 

0.381 

Sep  79 

727-100 

0.1243 

0.370 

Nov  79 

727-100 

0.1262 

0.390 

Mar  80 

727-100 

0.1183 

0.327 

Jun  80 

727-100 

0.1185 

0.335 

Apr  79 

0C-9-40 

0.1201 

0.821 

727-200 

0.0905 

0.344 

Apr  79 

Wide  body 

0.1407 

0.452 

DFW  6 

May  79 

727-100 

0.1056 

0.300 

Jul  79 

727-100 

0.1050 

0.290 
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TABLE  6- lb.  RESULTS  OF  AIRCRAFT  SIMULATIONS 


Line 

Date 

Ai rcraft 

Pilot  station  acceleration,  q 

RMS 

Peak 

DFW  7 

Jun  80 

727-100 

0.0146 

0.057 

QFW  8 

Jun  80 

727-100 

0.0324 

0.108 

DFW  9 

Jun  80 

727-100 

0.0711 

0.192 

DFW  10 

Jun  80 

727-100 

0.0557 

0.147 

DFW  11 

Jun  80 

727-100 

0.1296 

0.372 

DFW  12 

Jun  80 

727-100 

0.1330 

0.390 

Jun  80 

727-200 

0.0974 

0.397 

Jun  80 

DC-9-40 

0.0973 

0.441 

Jun  80 

Wide  body 

0.1452 

0.561 

DFW  13 

Jun  80 

727-100 

0.0317 

0.089 

DFW  14 

Jun  80 

727-100 

0.0243 

0.068 

Jun  80 

727-200 

0.0179 

0.060 

Jun  80 

DC-9-40 

0.0317 

0.096 

Jun  80 

Wide  body 

0.0261 

0.076 

NAS  1 

Jun  80 

727-100 

0.1326 

0.393 

Jun  80 

727-200 

0.1243 

0.468 

Jun  80 

0C-9-40 

0.1071 

0.450 

Jun  80 

Wide  body 

0.1508 

0.588 

NAS  2 

Jun  80 

727-100 

0.1495 

0.431 

Jun  80 

727-200 

0.1069 

0.422 

Jun  80 

DC-9-40 

0.1126 

0.474 

Jun  80 

Wide  body 

0.1651 

0.642 

NAS  3 

Jun  80 

727-100 

0.0805 

0.224 

Jun  80 

727-200 

0.0620 

0.230 

Jun  80 

DC-9- 40 

0.0844 

0.564 

Jun  80 

Wide  body 

0.1011 

0.331 

NAS  4 

Jun  80 

727-100 

0.0772 

0.199 

GAL  4 

Jun  80 

727-100 

0.0268 

0.085 

Jun  80 

727-200 

0.0339 

0.114 

Jun  80 

DC-9-40 

0.0398 

0.131 

GAL  5 

Jun  80 

727-100 

0.0556 

0.143 

Jun  80 

727-200 

0.0396 

0.148 

Jun  80 

DC-9-40 

0.0412 

0.215 

GAL  6 

Jun  80 

727-100 

0.0587 

0.158 

Jun  80 

727-200 

0.0430 

0.188 

Jun  80 

DC-9-40 

0.0512 

; _ 

0.374 
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TABLE  6-1 c.  RESULTS  OF  AIRCRAFT  SIMULATIONS 


Line 

Date 

Ai  rcraft 

Pilot  station 

acceleration,  g 

RMS 

Peak 

JSN  2 

l 

727-100 

0.0789 

0.220 

BW  SSfl 

727-200 

0.0560 

0.205 

iW 

DC-9-40 

0.0707 

0.468 

■ji  sis 

Wide  body 

0.0878 

0.324 

Jun  80 

727-100 

0.0776 

0.215 

JSN  3 

Apr  79 

727-100 

0.0825 

0.232 

Apr  79 

727-200 

0.0615 

0.236 

Apr  79 

DC-9-40 

0.0757 

0.487 

Apr  79 

Wide  body 

0.0927 

0.352 

Jun  80 

727-100 

0.0814 

0.242 

JSN  4 

Dec  79 

727-100 

0.1017 

0.277 

Dec  79 

727-200 

0.0537 

0.152 

JSN  5 

Jun  80 

727-100 

0.0934 

0.249 

JSN  6 

Jun  80 

727-100 

0.0575 

0.158 

727-200 

0.0575 

0.114 

DC-9-40 

0.0427 

0.226 

Wide  body 

0.0582 

0.175 

JSN  7 

Jun  80 

727-100 

0.1007 

0.244 

727-200 

0.0772 

0.286 

DC-9- 40 

0.0955 

0.627 

Wide  body 

0.0962 

0.449 

JSN  8 

Jun  80 

727-100 

0.0962 

0.218 

JSN  9 

Jun  80 

727-100 

0.0440 

0.102 

JSN  10 

Jun  80 

727-100 

0.0254 

0.086 

727-200 

0.0316 

0.110 

DC-9-40 

0.0392 

0.212 

Wide  body 

0.0367 

0.143 

JSN  11 

Jun  80 

727-100 

0.1465 

0.387 

JSN  12 

Jun  80 

727-100 

0.1288 

0.300 

JSN  15 

Jun  80 

727-100 

0.0991 

0.281 

727-200 

0.0698 

0.220 

DC-9-40 

0.1118 

0.517 

Wide  body 

0.1041 

0.285 

JSN  16 

Jun  80 

727-100 

0.0944 

0.220 

3 


the  table  gives  the  date  of  the  pavement  profile  observation.  The  pilot 
station  acceleration  was  used  for  evaluation  of  roughness,  and  0.4  g  was  the 
acceptable  acceleration  criterion  initially  used. 

The  acceptabil ity  of  the  pavement  surface  is  clearly  dependent  on  the 
aircraft  under  consideration.  The  same  pavement  may  produce  different 
responses  in  different  aircraft.  For  example,  at  NAS  3,  JSN  7,  and  JSN  15, 
the  smaller  aircraft,  DC-9-40,  exhibits  higher  peak  acceleration  response  in 
the  simulation. 

Reports  from  users  indicate  that  the  stretched  B727-200  has  more  trouble 
with  roughness  than  the  B727-100  model.  No  such  consistent  trend  is  evident 
in  the  data.  It  is  believed  the  TAXI  assumption  of  a  rigid  profile  may  con¬ 
tribute  to  this  discrepancy.  The  TAXI  simulation  estimates  are  not  as  good 
for  the  B727-200  as  for  other  aircraft.  Peculiarities  in  the  strut  for  this 
aircraft  are  not  modeled  properly  in  the  TAXI  program. 
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VII.  ANALYSIS  AND  CONCLUSIONS 


The  current  state  of  the  art  for  characterizing  expansive  soils  involves 
the  mathematical  description  of  a  three-dimensional  response  surface.  In 
this  scheme,  the  soil  volume  is  related  to  mechanical  load  and  suction  applied 
to  the  soil.  Natural  near-surface  expansive  clays  contain  other  features, 
such  as  cracks  and  seams,  that  affect  soil  behavior  significantly  and  must 
enter  into  the  characterization  process. 

HEAVE  PREDICTION 

In  Figure  5-9,  surface  elevation  chanqe  is  plotted  against  suction 
change.  Several  important  facts  were  observed:  (1)  At  two  points  the  actual 
soil  behavior  was  opposite  to  the  behavior  predicted  by  a  suction  model. 

(2)  The  error  variance  was  apparently  a  function  of  the  magnitude  of  suction 
change.  (3)  The  elevation  changes  are  serially  correlated  with  the  suction 
change,  as  indicated  by  the  Durbin-Watson  Test  (Durbin  and  Watson,  1951). 

These  observations  clearly  indicate  that  no  simple  model  can  describe  field 
oehavior  of  expansive  soils  adequately  for  pavement  design. 

The  two  deviant  points  in  Figure  5-9  indicate  that  surface  subsidence  may 
accompany  initial  wetting  in  cracked  expansive  clays.  This  surface  subsidence 
is  attributed  to  heave  caused  by  lack  of  lateral  restraint  on  the  soil  at  the 
cracks.  Although  it  is  a  principal  factor  in  soil  behavior,  lateral  restraint 
is  not  accounted  for  in  current  state-of-the-art  behavior  models.  These 
models  either  assume  a  fully  restrained  condition  or  arbitrarily  assign  a 
value  of  lateral  restraint. 

Lateral  restraint  may  be  thought  of  in  terms  of  a  factor,  f: 

(AL/L) 

i7-') 

where 

f  =  lateral  restraint  factor 

(AL/L)  =  vertical  linear  strain 
\V/V.j  =  volumetric  strain 

Durbin,  J.,  and  Watson,  G.  S.,  (1951)  "Testing  for  Serial  Correlation  in 
Least  Squares  Regression  II,"  i,  Vol  .  38:159-7?.!. 


The  f  value  will  vary  between  1,  for  fully  restrained,  and  1/3,  for  unre¬ 
strained  conditions. 

However,  roughness  is  the  primary  concern  of  the  present  research,  and 
more  accurate  prediction  of  heave  does  not  predict  roughness.  Therefore,  it 
is  reasonable  to  proceed  toward  a  characterization  of  wavelength  and  amplitude 
rather  than  the  prediction  of  heave. 

PREDICTION  OF  WEIGHTED  AMPLITUDE 

Weighted  amplitude,  A,  represents  a  weighted  value  of  the  amplitudes  of 
the  frequency  components  calculated  from  the  elevation  data.  The  weights  put 
more  value  on  longer  wavelength  components.  The  calculation  of  the  weighted 
amplitude  was  described  in  Table  5-4.  The  next  step  is  to  determine  a  method 
of  predicting  A  at  a  site. 

Figure  7-1  shows  the  weighted  amplitudes  for  the  uncovered  soil  studied 
in  the  field  experiment.  The  complexity  of  the  behavior  results  from  the  soil 
cracking.  Also,  the  suction  values  used  in  this  plot  were  determined  at  dis¬ 
crete  points,  whereas  the  A  values  are  a  function  of  data  taken  over  a  39-m 
(128-ft)  distance. 

If  the  soil  equilibrates  normally  beneath  a  newly  constructed  pavement, 
the  soil  volume  changes  only  once.  The  associated  suction  change  would  repre¬ 
sent  the  profile  moving  from  a  dry  condition  during  construction  to  a  condi¬ 
tion  of  equilibrium  beneath  the  pavement.  For  example,  DFW  ‘-eached  4.0  pF 
(103  kPa)  at  its  driest  point.  At  depth  the  suction  is  at  3.5  pF.  Normal 
equilibrium  would  involve  a  change  in  suction  of  0.5  pF.  Some  cases  warrant 
a  design  based  on  saturated  conditions.  Previous  work  (McKeen  and  Nielsen, 
1978)  indicated  an  end  suction  of  2.5  pF  would  be  appropriate  for  this  case. 
The  suction  change  would  then  be  1.5  pF. 

The  point  is  that  initial  and  final  suction  profiles  must  be  selected  by 
the  designer.  With  the  profiles  determined,  the  data  in  Figure  7-1  suggest 
the  following  relationships  for  predicting  A: 

Site  Slope  (AA) 

DFW  0.021  m  per  pF  (0.068  ft  per  pF) 

JSN  0.023  m  per  pF  (0.075  ft  per  pF) 

GAL  0.015  m  per  pF  (0.050  ft  per  pF) 
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3  4  5 

AVERAGE  SUCTION  IN  ACTIVE  ZONE,  pF 


FIGURE  7-1.  CHARACTERIZATION  OF  SOIL  BEHAVIOR  FOR  PAVEMENT  DESIGN. 
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The  relationship  between  A  and  suction  change  depends  on  several  factors. 
These  include  the  depth  of  the  active  zone  and  the  soil  response  to  suction 
change  as  i nd i ca t<  i  by  The  cracking  fabric  within  the  active  zone  as 
indicated  by  f  is  also  very  important. 

The  data  shown  in  Figure  7-1  are  the  result  of  the  field  experiment  con¬ 
ducted  in  this  study.  Clearly  each  of  these  sites  exhibits  distinctive  fea¬ 
tures.  Future  work  should  be  directed  to  obtaining  more  data  of  this  type  at 
a  variety  of  sites. 

WAVELENGTH  AND  ROUGHNESS 

The  amplitudes  must  be  accompanied  by  a  horizontal  distance  in  order  to 
be  useful  in  the  study  of  roughness.  To  characterize  the  wavelength,  the 
autocorrelation  function  for  the  elevation  data  was  calculated.  The  horizon¬ 
tal  distance  at  which  the  autocorrelation  function  goes  to  zero  represents  a 
characteristic  of  the  cyclic  nature  of  the  surface.  With  this  distance,  it  is 
possible  to  quantify  the  profile  in  terms  of  roughness  parameters,  an  ampli¬ 
tude  and  a  wavelength. 

Wavelengths  and  amplitudes  for  the  uncovered  soil  surfaces  studied  are 
given  in  Table  7-1.  Note  that  the  sites  with  deeper  active  zones  (DFW,  JSN) 
exhibit  a  more  constant  wavelength.  The  shallow  active  zone  (GAL)  shows  more 
variation  in  wavelength.  When  a  pavement  is  placed  on  these  soils,  this 
natural  wavelength  is  changed.  This  change  must  be  understood  in  order  to 
proceed  with  the  proposed  approach  to  pavement  design.  Wavelength  and  ampli¬ 
tude  data  for  the  paved  surfaces  studied  are  presented  in  Tables  7-2  through 
7-4.  The  layers  in  the  pavements  were  converted  into  equivalent  depths  of 
asphalt  concrete  for  ready  comparison.  The  equivalent  depth  of  each  of  these 
pavements  is  shown  in  Table  7-5. 

The  equivalent  depth  was  calculated  by  first  assuming  modulus  values  as 
follows:  Portland  cement  concrete,  2.8  x  107  k Pa  (4  x  106  lb/in2);  asphaltic 
concrete,  3.5  x  10f'  kPa  (5  x  10s  lb/in2);  cement-treated  base,  1.4  x  10s  kPa 
(2  x  10s  lb/in2);  lime-stabilized  subgrade,  2.1  x  10"  kPa  (3  x  103  lb/in2). 
Each  layer  thickness  was  then  converted  to  an  equivalent  thickness  of 
asphaltic  concrete  using  the  modular  ratios.  The  sum  of  these  thicknesses  in 
the  pavement  is  the  equivalent  depth. 


TABLE  AMPLITUDE-WAVELENGTH  FOR  SOIL  SURFACES 


Pavement 

Ampl i tude , 
m  (  f  t ) 

Wavelength 

m 

DFW  1 

JAN 

0.0421 

(0. 138) 

APR 

0.0454 

(3.149) 

8.2 

MAY 

0.0649 

(3.213) 

7.0 

JUL 

0.0454 

(3.149) 

8.2 

SEP 

0.0527 

(0.173) 

8.2 

NOV 

0.0369 

(3.121) 

8.2 

MAR 

0.0552 

(0.181) 

JUN 

0.0427 

(3.  140) 

JSN  1 

JAN 

0.0399 

(9.131) 

APR 

0.0363 

(3.119) 

11.0 

MAY 

0.0405 

(3.  133) 

10.7 

JUL 

0.0430 

(9.141) 

10.7 

SEP 

0.0408 

(9.134) 

9.5 

DEC 

(9.134) 

10.4 

JUN 

■ 

,  ’  .  |  j  -  -  . 

(9.154) 

GAL  1 

NOV 

0.0146 

(3.048) 

MAR 

0.0204 

(9.067) 

7.0 

MAY 

0.0183 

(3.060) 

4.6 

JUL 

0.0110 

(0.036) 

7.0 

NOV 

0.0046 

(9.015) 

3.1 

JUN 

0.0189 

(9.062) 

(27) 

(23) 

(27) 

(27) 

(27) 


(36) 

(35) 

(35) 

(31) 

(34) 


(23) 

05) 

(23) 

(10) 
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TABLE  7-2.  AMPLITUDE-WAVELENGTH  DATA  FOR  DFW  PAVEMENTS 


Pavement 


Amplitude, 
m  (ft) 


Wavelength, 
m  (ft) 


0.014 

0.020 

0.013 

0.018 

0.012 

0.013 

0.013 

0.014 

0.015 


0. 0070 
0.0067 
0.0085 
0.0073 
0.0076 
0.0082 
0.0070 
0.0058 
0.0079 


0.027 

0.020 

0.023 

0.018 

0.022 

0.025 

0.024 

0.024 

0.022 


0.015 

0.014 

0.014 

0.014 

0.014 

0.018 

0.014 

0.015 

0.014 


0.046 

0.064 

0.044 

0.058 

(0.040) 

(0.043) 

(0.042) 

(0.046) 

(0.050) 


(0. 023) 
(0.022) 
(0.028) 
(0.024) 
(0.  025) 
(0.027) 
(0.023) 
(0.019) 
(0.026) 


(0.087) 

(0.064) 

(0.076) 

(0.059) 

(0.073) 

(0.083) 

(0.079) 

(0.077) 

(0.073) 


(0.049) 
(0.  045) 
(0.047) 
(0.045) 
(0.046) 
(0.059) 
(0.  047) 
(0.048) 
(0.046) 


TABLE  7-4.  AMPLITUDE-WAVELENGTH  DATA  FOR  ROUGH  PAVEMENTS 


Pavement 

Amplitude, 
m  (ft) 

Wavelength, 
m  (ft) 

DFW  7 

0.0314 

(0.103) 

2. 7-3.7 

(9-12) 

8 

0.0250 

(0.082) 

9.1 

(30) 

9 

0.0454 

(0.149) 

7.9 

(26) 

TO 

0.0469 

(0.154) 

5.8 

(19) 

11 

0.0600 

(0.197) 

8.2 

(27) 

12 

0.0393 

(0.129) 

7.3 

(24) 

13 

0.0351 

(0.115) 

5.9 

(19.5) 

14 

0.0344 

(0.113) 

7.3 

(24) 

NAS  1 

0.0433 

(0.142) 

3.7-12.2 

(12-40) 

2 

0.  0323 

(0.106) 

4.9 

(16) 

3 

0.0122 

(0.040) 

3. 2-9.8 

(10.5-32) 

4 

0.0140 

(0.046) 

5.6 

(18.5) 

GAL  4 

0.0527 

(0.173) 

4.0 

(13) 

5 

0.0223 

(0.073) 

6.2 

(20.5) 

6 

0.0143 

(0.047) 

11.3 

(37) 

JSN  4 

0.1067 

(0.350) 

6.4 

(21) 

5 

0.0518 

(0.170) 

6.2 

(20.5) 

6 

0. 0326 

(0.107) 

5.2 

(17) 

7 

0.0469 

(0.154) 

8.2 

(27) 

8 

0.0649 

(0.213) 

8.5 

(28) 

9 

0.0271 

(0.089) 

7.3 

(24) 

10 

0.0311 

(0.102) 

5.5 

(18) 

11 

0.0588 

(0.193) 

5.8 

(19) 

12 

0.1006 

(0.330) 

6.1 

(20) 

15 

0.1070 

(0. 351) 

6.4 

(21) 

16 

0.0866 

(0.284) 

7.0 

(23) 

TA8LE  7-5.  PAVEMENT  EFFECTIVE  DEPTH 


Pavanent 

Surface 

material 

Equivalent  depth,  1 
m  ( <  n ) 

i 

DFW  2 

AC 

0.115 

(4.5)  | 

DFW  3 

PCC 

0.820 

(32.3)  : 

DFW  5  through 

14 

AC 

0.200 

17.9) 

NAS  1,2 

PCC 

0.549 

(31 .6) 

NAS  3,4 

PCC 

0.511 

(2C.1) 

GAL  2  through 

6 

AC 

0.400 

<15.7) 

JSN  2,3 

AC/ PCC 

0 . 600 

•  •  r  \ 

' .  0  < 

JSN  5-9 

PCC 

0.518 

f.:u.4) 

JSN  10 

AC 

0.857 

'ic.D 

JSN  11,12 

PCC 

0.625 

<2:  .  ! 

JSN  4,15,16 

L 

PCC 

!  n.,101. 

i...  -  --  -  . 

(12.0; 

Notes:  AC  is  asphalt  concrete;  PCC  is  pot  1 1  arid 
cement  concrete. 


A  plot  of  equivalent  depth  versus  wavelength  for  the  pavements  at  GAL, 
DFW,  and  NAS  is  shown  in  f  igure  7-2.  Figure  7-3  is  a  similar  plot  tor  the  JSN 
area.  These  data  reveal  that  the  wavelengths  on  the  pavement  edges  and  those 
in  the  interior  differ  markedly.  In  Figure  7-2,  the  edges  contain  snorter 
wavelengths  whereas  the  interior  lines  exhibit  wavelengths  equal  to  or  greater 
than  the  uncovered  soil  wavelength  (\<>).  The  data  m  F inure  7-3  show  that  the 
reverse  is  true  of  the  pavements  studied  in  the  J:>N  are.',  'ins  difference  is 
interpreted  as  evidence  of  a  fundamental  difference  in  so i 1 -pavement  behavior 
at  these  test  sites. 

The  shorter  wavelengths  exhibited  at  the  pavement  edges  in  figure  7-2 
result  from  climatic  changes.  The  soil  wet-dry  sequence  produces  more  move¬ 
ment  near  the  edge  than  in  the  interior,  resulting  in  shorter  wavelengths.  It 
is  clear  that  .the  design  condition  should  he  the  pavement  edge,  where  cyclic 
movements  occur  and  produce  short-wavelength  roughness. 

In  contrast,  the  JSN  data  reveal  cyclic  movements  under  the  interior 
pavement  areas.  Three  possible  mechanisms  are  postulated:  (a)  moisture  move¬ 
ment  through  the  pavement  surface,  i.e.,  pavement  permeability  greater  than 
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FIGURE  7-3.  DEPTH-WAVELENGTH  RELATIONSHIPS  AT  JSN. 
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subgrade  permeability;  (b)  cyclic  suction  changes  occurring  via  avenues  within 
the  profile  such  as  seams  and  fissures;  (c)  changes  in  the  osmotic  suction  in 
the  profile  which  are  not  related  to  climatic  variation. 

Because  shorter  wavelengths  produce  more  severe  roughness,  they  are  an 
attractive  design  parameter.  The  data  in  Figures  7-2  and  7-3  indicate  that 
the  design  wavelength  (A^)  is  about  equal  to  half  the  soil  wavelength  (0.5  ■  ) 
for  airport  pavements.  This  approximation  will  be  applied  in  tne  next 
section. 

ACCEPTABILITY  CRITERIA 

The  use  of  wavelength  and  amplitude  characteristics  for  the  design  of 
pavements  requires  the  establishment  of  acce.  Uhil  i  t>  •;.!  :  tori  a.  Acceptable 
pavement  roughness  depends  on  acceptable  levels  of  aircraft  rosoonse. 

Although  no  widely  accepted  aircraft  criterion  is  available,  a  pilot  stat 
vertical  acceleration  of  0.4  g  has  been  used  most  often,  i he  TAXI  computer 
code  was  first  used  to  simulate  aircraft  response  for  the  puiyose  of  estab¬ 
lishing  acceptability  quantitatively.  The  application  of  t.l,e  0.4-g  criteria!’, 
was  not  entirely  satisfactory  for  the  following  reason..  (1)  Face  aircraft 
responds  differently  to  a  pavement  surface;  only  four  commercial  aircraft  were 
used  in  the  simulation  in  this  study.  (2)  As  the  ornss  weight  of  the  aircraft 
varies,  its  response  to  the  pavement  surface  varies.  (3)  The  pavement 
deflection  under  aircraft  loads  contribute  to  the  roughness  encountered  by 
real  airplanes.  TAXI  does  not  model  this  behavior.  (4/  rne  B727--20P  has 
not  been  modeled  as  well  as  other  aircraft  because  the  TAXI  simulations  do 
not  properly  model  its  specific  strut  characteristics. 

In  addition,  several  pavements  clearly  unacceptable  to  users  and  airport 
operators  did  not  produce  unacceptable  accelerations  in  the  aircraft  simula¬ 
tions.  This  discrepancy  may  result  from  variations  n  the  actual  operating 
conditions,  such  as  the  gross  weight  and  speed  of  the  aircraft.  Another  fac¬ 
tor  involved  in  complaints  by  users  is  the  appearance  of  the  pavement  profile. 

The  acceptability  was  therefore  determined  in  two  ways.  First,  an 
acceleration  criterion  of  0.3  g  was  used.  It  was  reduced  from  0.4  g  because 
it  is  intended  for  design  use.  The  reduction  was  made  based  on  an  examination 
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of  the  data  in  this  study  and  on  experience.  The  second  criterion  was  subjec¬ 
tively  determined  and  involved  rating  the  acceptability  of  the  pavement  as 
expressed  by  users  and  operators.  Visual  inspection  of  the  pavement  and  study 
of  actual  elevations  were  also  used. 

Acceptability  data  are  shown  in  Figure  7-4  for  B727-100  and  -200  and  in 
Figure  7-5  for  the  DC-9-40  and  TriJet  Composite.  The  lines  on  the  graphs 
represent  criteria  for  acceptability  of  the  amplitude  and  wavelength  combina¬ 
tions.  The  acceptability  for  the  DC-9-40  and  TriJet  Composite  was  more 
conservative  and  was  selected  for  further  use  in  the  present  study. 
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FIGURE  7-4.  ACCEPTABILITY  OF  PAVEMENTS  FOR  THE  B727-200  AND 
B727-1 00. 
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b'lOMMf.  'J[v*T  IONS 


The  findings  of  the  field  st-cu  ;•*•••  vic<>  .1  basis  tor  i-p-.-owenciat ion1’  '• 
the  design  of  airport  pavements  is  expansive  soil  .  A  ; roposed  desim  pro¬ 
cedure  was  shown  ’n  Figui e  1-2.  in  ihi  section,  ice  findings  of  the  field 
experiment,  will  be  applied  to  this  ;  rnceduro. 


INDICATOR  TEST 

An  indicator  tost  may  be  used  so  determine  whether  t. f i ■  -  s. c ■  i  1  re  :u i  1 
special  expansive  soils  design  procedure?..  In  -an.  cast  -,  thi ?  d-t'n-::nr  - 
will  not  require  a  oh.-sita!  test,  bat  will  be  isr‘e  on  sue  :iar-is  of  e-.,.,  ■” 
or  available  data.  In  tit  absent  e  of  other  rvi  !e;e’.  sm  out.-'  ■- '  u,-.*  ?  ■  r 
undisturbed  soil  response  te  '  action  cnanges  ;  re" i.:  *■  be  *  mnii. 

of  potential  soil  behavior.  Suction  tests  pro'.  1  i<-  iirf>M,-:  vul-.es  dope' 
on  the  lateral  restraint  imposed  in  the  tei  •  : 

( 'V/V,  If 

h  log-.,  11  f 

h . 

7 


where 

=  suction  compression  index 
.•‘,V/V.  =  volume  change  in  terms  of  origin.)  1  volume 
hf,  h.  -  f  i  na  1  and  ini tia !  -  u'-'ion  in  ari  thm*w  «  .civ 

T  1 

f  -  lateral  restraint  factor 


The  clod  test  and  the  COLL  tes4 ,  foe  example,  p*"vid(-  v.<l  .u\.  of  ;  ■ .  i- .  no 

restraint  (f  =  1/ 3) . 

The  values  of  .  r,  shown  in  f  in  me  )  are  for  ;  I  ?  •  raint  factor 

f  =  1/3.  Either  <  va  I  ue  mav  he  used  11:  ca  Iculat  :orts  as  Iona  as  the  lateral 

restraint  factor  is  expressed  explicitly.  1  he  values  obtained  frost  the  chart 
in  Figure  4-2’/  are  for  r  -  1.0. 

SOIL  CHARACTER  I /AT IDN 

The  two  -teir,  shown  in  i  igum  I  •  ..  for  swell  eve  I  uat  inn  and  heave  predic¬ 
tion  should  be  <  muni  n -o  mV)  .  single  1  ej  <  ill.-1  .  i’re- 

cedures  for  obta’n-'n  ;  1  . na  ■  .1 1  t  ■  1  i  s  t.  i .  .  a  re  .  r  1  i--d  P'lcnv. 


1 1'.') 


Predict  i.it.  -oil  airp'i  itude  arid  wave lenytn  is  the  fir-i  -tep.  Ir.c- 
wavelenqtri  . ,  a  j  i  rec  t  f  uric  liar  of  tne  active  zone  depth  tor  s  it-  silt, 
studied  (Hyjre  8-  i , .  Active  zone  deptn  fur  .jSN  was  inferred  from  tie  tuCt 
that  interior  gavement  areas,  rutncr  Inan  the  edges,  were  undergo  in-.  move  sent . 
Interior  mover  rent  suggt*-  ts  that  tne  principci  tuu..*;  if  distortion  ;<  not  sur¬ 
face  climatic  variation,  it  is  the  movement  of  water  via  discontinuities  i" 
the  soil.  Data  in  1  .  j.n>-  a-  •  •  .  ■•-’i  .  t-  ■  !' ..cum  inn  .  ties  at  3./  :  'la  -t;  d 
the  test  site.  The  data  F  i..i.re  b-  ;  wan  oe  useo  for  estimating  for  too 

soil  at  a  sue.  Vneir  use  ;s  net  recommenuoa  outside  of  tne  given  data  range 
or  in  the  presence  a  „.:ler  t»i  i ■: .  Wavelength  may  .J  so  oe  estima tea  by 
surveying  the  surface  ,  ccr:ie  and  calculating  tne  docurre  i tier,  otai.ce  for 
the  data.  Trie  '  -  S  t-n.:  distance  at  wmci.  tire  auto 

correlation  fun. lion  got;.  to  zer  ■ 

The  o  trier  oil  erooert  v  ;  equi  red  it  trie  wei  giited  amplitude  T  'if-  slope 
of  A  versus  act  race  ..  f :  v  c  cone  sue  lion  ['  ic  are  7-1 )  i  s  the  test  estii  alcr  for 
the  sites  st.udie-j.  Thus,  an  esti-ate  of  in  *  and  final  .octun  t.rof  ilcs 
is  required.  bits  t'-.i  s  e-,t  --ite.  mo  .mi  ion  ten  m  <  an  :e  n  .eu  to  :.-'edict 
the  A.  Tfi  i  s  prediction  :  •  fur  a  oire-ti.ie  mange  ;  roi ,  initial  to  tim: ,  If 
the  suction  profile  f  1  ,.t. Luaies,  the  A  mi;!  nv'ohab1;.  ouilG  up  from  c.'.ele  to 
cycle,  result!'  '  i>  a  '.an,  >•  A. 

For  site-  otuer  than  studied,  a  meidon  *‘or  prediction  ot  A  u 

desired.  '-Vasu  reiser  it  and  m  l ;  mat  ion  are  .lternative  prediction  procedures. 

Measure'nen  t 

Witn  tins  proceuui  •.  ,  elevition  data  are  gathered  and  A  c  .imputed  at  two 
or  three  times  of  the  year.  These  data  are  plotted  against  suction,  os  shown 
in  Figure  7-1,  and  tne  slope  of  tne  least,  square  f't  is  determined.  The  fol¬ 
lowing  must  be  considered  : n  using  this  method.  Daring  the  initial  swell  of 
a  cracked  soil  profile,  dm  reuses  in  A  are  likely.  The  area  studied  must  be 
undisturbed  soil  in  order  to  reflect  the  benavior  ot  the  active  zone  as  it 
interacts  with  the  climaie.  Preferably,  the  area  should  be  relatively  flat. 
These  factors  will  combine  to  produce  a  worst-case  or  greatest  amplitude  soil 
for  the  study. 
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FIGURE  8-1.  ESTIMATION  OF  SOIL  WAVELENGTH. 
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WAVEI  ENG7W 


Estimation 


Soil  properties  may  be  used  to  predict  the  slope  of  the  amplitude-suction 
relationship.  On  the  basis  of  the  factors  involved  in  the  swelling  of  a  clay 
profile,  A  is  a  function  of 


•  suction  compression  index  (,.  ) 

•  active  zone  depth  (Z  ) 

d 

•  variation  of  suction  compression  index  from  point  to  point 
■  variation  of  suction  profile  from  point  to  point  (h^,  h^ ) 

The  development  of  an  estimate  for  a  soil  follows  these  steps: 


Vnax  'hj°^'h. 

!  1 


(hf  i 

'“fid 


J  ‘ 

— 

__ 

^  Kax  -  iHmin]  *  1/2 

IV 

-  .hr'09 

hfl 

T~\ 

h .  | 

k. 

1 

(8-2) 


Now  it  is  necessary  to  obtain  a  measure  of  the  variability  of  vh  and 
hf/h.  in  the  profile.  So  fa.  it  has  been  assumed  that  soil  beneath  pavement 
areas  is  in  equilibrium  after  swell.  Therefore,  h^/h-  will  not  vary,  although 
the  suction  will  vary  at  intermediate  times.  The  case  for  no  variation  in 
suction  was  calculated  and  is  illustrated  in  Figure  8-2. 


SOIL-PAVEMENT  INTERACTION  MODEL 

When  a  pavement  is  placed  on  an  expansive  soil,  pavement  stiffness  can  be 
expected  to  reduce  the  A.  In  order  to  account  for  this  effect,  the  following 
model  is  proposed. 


The  effect  of  pavement  stiffness  in  smoothing  expansive  clay  roughness 
is  modeled  using  a  beam  on  an  elastic  foundation  (Figure  8-3).  By  summing 
vertical  forces  on  the  element,  one  obtains 

Q  -  (Q  *■  dQ)  +  kb(w  -  w')dx  -  qbdx  =  0  (8-3) 

^  =  kb(w  -  w' )  -  qb  (8-4) 

where 

Q  =  shear 
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COEFFICIENT  OF  VARIATION  OF  y.  PERCENT 


FIGURE  8-2.  PREDICTION  OF  A  FROM  y.  AND  ITS  VARIATION 


I*  -  - - x - i 


q  -  unit  weight  of  beam  material 
k  =  subgrade  modulus 
b  =  width  of  beam 

a  =  hal f-ampl i tude  of  deformed  soil 
h  =  thickness  of  beam 
w  =  deflection  of  beam 
w'  =  shape  of  soil  surface 
Q  -  shear 
a  =  wavelength 
dx  =  incremental  length 


FIGURE  8-3.  BEAM  ON  DEFORMED  ELASTIC  FOUNDATION  PROBLEM 


w  =  deflection  of  beam 
w'  =  shape  of  soil  surface 
b  =  width  of  beam 
dx  =  incremental  length 
k  =  subgrade  modulus  (explained  below) 
q  =  unit  weight  of  beam  material 


Because 


d([  =  dj  M 

dx  , 
dx 


M  =  -El  ^ 


where 

E  =  modulus  of  elasticity  of  beam 
I  =  moment  of  inertia  of  the  section 
M  =  moment 

we  can  rewrite  the  expression  as 

^  =  -El  d--  =  kb(w  -  w '  )  -  qb 
dx1’ 

Now  an  expression  for  w  is  assumed: 


(8-5) 


w  =  Ae~'xcos  dx  +  Be_4,xsin  dx  +  Ce^’xcos  fc>x  +  DehXsin  fix  +  ^  +  w' 


where 


A,  B,  C,  D  are  constants 


=  a  characteristic  stiffness 


(8-6) 


Taking  the  necessary  derivatives  and  substituting  into  Equation  8-5 
yields 

— -  -  -4Ak')’’e  1  xcos  rx  -  4Bd4e  1  xsin  (3x  -  4Cis4e,xcos  Bx 

dx4 


-40.-  'e'  xs  i  n  ,-.x  +  d  w- 


(8-7) 
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Then  the  following  boundary  conditions  are  imposed  on  the  function  w': 


At  x  =  0,  w‘  =  0,  =  0 

At  x  =  w'  =  a 


At  x 


2a 


dV 

dx4 


=  0 


a  =  the  half-amplitude  of  the  soil 


Beginning  with  the  last  expression  and  integrating  gives 
w‘  =  +  r*-  +  Gx  +  H 

With  the  boundary  conditions  above,  the  constants  (E,  F,  G,  H)  may  be 
evaluated. 

H  =  G  =  0 

(8-8) 
(8-9) 

Solving  these  equations  simultaneously  yields 


E  f  0  3 

6  4. 


=  a 


EfA. 

6|2j 


3  +  I 

'  r\ 


=  2a 


w'  =  -32a  [*] 3  +  24a  W  " 

(8-10) 

In  solving  the  equation  for  the  beam  deflection, 

the  following  boundary 

conditions  are  imposed: 

n  dw  n 

"'•■dr0 

(8-11) 

X  -  0.  ^  =  0 

dx3 

(8-12) 

A  dw  n 

x  2’  dx  ~  0 

(8-13) 

A  d3w  „ 
x  =  -o,  -  =  0 

2  dx3 

(8-14) 
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The  boundary  conditions  8-11  through  8-14  produce  a  system  of  four 
equations  and  four  unknowns: 


dw' 


(-A  +  B)8  +  (C  +  D)6  =  - 
2[A  +  B]63  +[-C  +  D]83  =  - 


x=0 
d3w‘ 


dx; 


x=0 


(-A  +  B)Bbi  +  (-A  -  B)8b2  +  (r  +  D)8b3  +  (-C  +  D)Bb4  =  -  jjj' 
(A  +  B)2(33bi  +  (-A  +  B)203b2  +  (-C  +  D)2B3b3  +  (-C  -  D)233b4 
d3w' 


dx3 


x  = 


(8-15) 

(8-16) 

(8-17) 

x  =  X 
2 

(8-18) 


where  _X 

"Pjf  x 

bj  =  e  cos 

fiA 

PT  \ 

b2  =  e  sin  3 j 

A 

b3  =  e  cos  fe- 

oA 

p?  X 

b4  =  e  sin  Bo 


This  system  of  four  equations  and  four  unknowns  may  be  solved  by  the  use 
of  matrix  methods  which  provide  values  for  A,  B,  C,  and  D.  Then  the  maximum 
pavement  elevation  difference  is 


Aw 

2a 


Aw  =  w 
max 


=  A 
max  a 


A  "  I*  =  0 

x  =  2 


B 

[bA 

1  +  i 

[b3  -  V, 

+ 

fb4l 

a 

; 

1  a 

l  2  J 

a  I 

1.2  J 

+  1 


(8-19) 


AW 


For  ease  of  use,  a  plot  of  versus  BA  was  prepared  for  a  variety  of 
values  of  stiffness  (Figure  8-4). 
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EXPANSIVE  SOIL  DESIGN  PROCEDURE 

Figure  8-5  outlines  the  revised  design  procedure  for  airport  pavements 
on  expansive  soils.  It  may  be  viewed  as  a  three-step  procedure.  Step  1  is 
the  determination  of  the  need  for  special  expansive  soil  design.  Step  2  is 
soil  characterization  in  terms  of  parameters  necessary  for  expansive  soil 
design.  Step  3  is  a  study  of  soil -pavement-aircraft  interaction  to  arrive 
at  a  structurally  adequate  thickness  for  the  pavement.  These  steps  are 
described  in  detail  below.  Examples  of  numerical  values  for  the  sites  studied 
are  also  provided.  Important  assumptions  and  limitations  are  listed. 

Indicator  Test 

The  suction  compression  index  (y^)  is  recommended  for  use  as  an  indicator 
of  soil  swell  potential.  In  addition,  the  use  of  judgement  and  experience  is 
necessary.  In  many  cases,  judgement  and  experience  will  eliminate  the  need 
for  an  objective  system  such  as  the  rating  criterion. 

The  for  a  particular  site  may  be  measured  by  any  one  of  several  tests, 
including  the  COLE,  clod  test,  and  standard  consol idometer- swell  test.  A 
chart  for  indirect  determination  of  yh  has  been  developed  and  is  presented  in 
Figure  4-27.  The  data  in  Figures  4-25  (COLE)  and  4-26  (clod)  indicate  very 
few  values  below  0.06,  placing  all  of  sites  studied  in  the  high  swell 
category  (Figure  8-5). 

Soil  Characterization 

Once  a  soil  is  placed  in  the  high  category,  special  procedures  are 
required.  First,  a  number  of  soil  parameters  which  are  compatible  with  the 
interaction  model  are  determined.  Then  these  parameters  are  used  to  evaluate 
the  soi 1-pavement  interaction. 

Active  zone  depth  (Z  )--The  depth  of  soil  that  interacts  with  the  surface 

_  a 

climatic  conditions  is  called  the  active  zone  depth.  Its  determination  is  not 
simple.  The  measurement  of  suction  versus  depth  at  several  times  of  the  year 
is  probably  the  best  indicator  of  Z=.  However,  it  should  be  clear  that 

a 

periods  of  extreme  weather  (wet  or  dry)  may  drive  the  seasonal  influence  to 
greater  depths. 
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FIGURE  8-5.  REVISED  DESIGN  PROCEDURE. 
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From  the  suction  data  in  Figures  4-17  and  4-18,  the  active  zone  depth  at 
DFW  was  estimated  to  be  2.1  m  (7  ft).  Data  in  Figure  4-19  for  GAL  indicate 
that  the  Z,  was  1.2  m  (4  ft).  The  JSN  site  did  not  undergo  a  normal  climatic 

a 

cycle  during  this  study;  the  soil  profile  remained  wet  for  the  duration  of 
the  observations.  The  approach  taken  was  to  determine  the  depth  at  which  the 
suction  profile  reached  a  constant  value.  The  data  in  Figures  4-20  and  4-21 
indicate  that  a  depth  of  3.7  m  (12  ft)  Is  appropriate  for  the  JSN  active  zone 
depth. 

Subgrade  modulus  (k)— This  modulus  represents  the  reduction  of  swell  by 
loads  applied  at  the  soil  surface.  Requirements  for  determining  k  are  a  swell 
value  at  little  or  no  load  and  a  swell  value  at  loads  approximately  equal  to 
the  pavement  load.  It  is  the  pavement  dead  load  that  will  tend  to  reduce  the 
swell  of  the  underlying  soil.  Values  of  this  modulus  are  not  comparable  to  the 
conventional  modulus  of  subgrade  reaction.  In  the  present  study,  k  values 
were  obtained  by  conducting  consolidometer-swel 1  tests  at  a  token  load  and 
then  consolidating  the  sample  in  normal  increments.  The  modulus  was  calcula¬ 
ted  from  the  percent  swell,  sample  height,  and  the  load  difference.  Values 
obtained  were  as  follows: 

Site  k,  kPa/m  ( 1 b/i n 3 ) 

GAL  1.65  x  105  (608) 

DFW  3.56  x  105  (1310) 

JSN  2.26  x  105  (537) 

Suction  change  (Ah)--When  expansive  soils  are  covered,  they  no  longer  can 
interact  with  the  climate  at  their  surface.  The  normal  result  is  moisture 
equilibration  with  the  soil  below  the  active  zone.  A  water  table  or  some 
other  feature  of  the  profile  which  influences  water  movement  (e.g.,  seams  and 
fissures)  will  alter  the  equilibration.  The  selection  of  the  suction  change 
during  equilibration  is  not  a  well-defined  procedure.  Some  workers  recommend 
using  saturated  final  conditions  while  others  propose  numerous  alternatives, 
depending  primarily  on  their  particular  climate  and  drainage  conditions.  In 
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(he  present,  work,  the  driest  suction  profiles  are  about  1.5  pr,  and  the 
'ipf  -c^t  about  2.5  |jf  .  These  extremes  were  selected  as  initial  and  final 
profile:,  because  there  is,  at  present,  no  rational  method  for  the  selection 
of  these  values.  Use  of  these  extremes  results  in  the  greatest  reasonable 
Ah,  Therefore,  for  the  present  examples.  Ah  of  2.0  pF  is  used  for  design. 

boil  amplitude  (A)- -Data  from  the  field  experiment  were  used  to  compute 
the  weighted  amplitude  of  the  frequency  components  of  the  surface  profile. 
These  computations  were  discussed  in  Section  V.  The  values  of  A  are  plotted 
versus  the  change  in  suction  in  Figure  7-1.  In  some  cases,  because  of  crack¬ 
ing,  the  results  do  not  follow  a  linear  relationship.  In  each  case,  however, 
an  estimate  of  the  A  versus  Ah  slope  (AA)  was  made.  These  values,  together 
with  estimated  suction  changes,  were  used  to  predict  values  of  A  for  each  of 
the  sites  studied.  Results  of  these  computations  are  0.020  in  (0.100  ft)  at 
GAL,  0.042  m  (0.136  ft)  at  DFW,  and  0.046  m  (0.150  ft)  at  JSN. 

Soil  wavelength  (")--A  characteristic  wavelength  of  the  soil  is  required 
to  complete  the  roughness  characterization.  The  decorrelation  distance  of  the 
autocorrelation  function  of  the  elevations  measured  was  used.  This  distance 
may  not  be  the  actual  fundamental  wavelength  of  the  soil,  but  it  is  some 
multiple  of  that  wavelength  and  therefore  serves  the  present  purpose.  Data 
were  plotted  versus  active  zone  depth  in  Figure  8-1,  and  may  bt  used  to 
select  a  value  of  ac  for  the  soil.  In  Figures  7-2  and  7-3,  values  of  \0  for 
pavement  sections  were  plotted  versus  the  equivalent  depth  of  each  pavement. 

It  was  recognized  that  pavement  structures  caused  shorter  wavelengths  than 
those  that  occur  on  the  uncovered  soil.  The  design  wavelength  was  determined 
to  be  half  of  the  uncovered  soil  wavelength.  Values  of  the  design  wavelength 
(A^j)  are  therefore  2.7  m  (9  ft)  at  GAL,  4.0  in  (13  ft)  at  DFW,  and  5.2  m 
(17  ft)  at  JSN.  Results  of  the  soil  characterization  for  all  sites  in  shown 
in  Table  8-1 . 

interaction  Study 

The  next  step  is  to  determine  the  amount  of  stiffness  required  in  a  pave¬ 
ment  section  to  resist  the  development  of  excessive  pavement  roughness  caused 
by  soil  swelling.  The  pavement  was  modeled  as  a  beam-on-elastic  foundation 
with  the  foundation  deformed.  Through  a  combination  of  beam  bending  and  soil 
compression,  the  amplitude  of  the  soil  deformation  was  reduced.  The  objective 
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TABLE  8-1.  RESULTS  OF  SOIL  CHARACTERIZATION 


Parameter 

GAL 

,  DFW 

JSN 

Indicator  test 

High 

High 

High 

Active  zone  depth  (Z  ), 

1.2 

2.1 

3.7 

m  (ft)  a 

(4) 

(7) 

02) 

Modulus  (k),  kPa/m 

1.65  x  105 

3.56  x  105 

1.45  x  105 

(lb/in3) 

(608) 

(1310) 

(537) 

Suction  change  (A^) ,  pF 

2.0 

2.0 

2.0 

Soil  amplitude  (A  ), 

30.5 

41.5 

45.7 

mm  (ft) 

(0.100) 

(0.136) 

(0.150) 

Soil  wavelength  (A0), 

5.5 

7.9 

10.4 

m  (ft) 

08) 

(26) 

(34) 

Design  wavelength  (A,), 

2.7 

4.0 

5.2 

m  (ft)  d 

(9) 

03) 

07) 

in  design  is,  then,  to  reduce  the  amplitude  in  the  pavement  surface  to  a  toler¬ 
able  amount.  Acceptable  values  of  amplitude  were  determined  from  the  TAXI 
computer  code  and  are  shown  in  Figure  7-5  for  four  aircraft  (B727-100, 

B727-200,  DC-9-40,  TriJet  Composite). 

Results  of  the  differential  equation  are  plotted  in  Figure  8-4  in  non- 
dimensional  quantities  BA  and  Aw/2a.  In  order  to  simplify  the  problem,  the 
pavement  is  considered  in  terms  of  its  equivalent  depth  D,  where  all  pavement 
layers  are  converted  to  an  equivalent  thickness  of  asphaltic  concrete  with  a 
modulus  (E)  of  1.95  x  106  kPa  (5  x  10s  lb/in2). 

Acceptable  values  of  amplitude  were  selected  from  the  solid  line  in  Fig¬ 
ure  7-5  according  to  the  design  wavelength  for  each  site.  Results  are  in 
Table  8-2.  Then  the  acceptable  amplitude,  A.  was  divided  by  the  estimated 

a 

soil  amplitude,  Ag,  to  provide  the  required  Aw/2 a.  Using  the  model  results 
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TABLE  8-2.  RESULTS  OF  INTERACTION  STUDY 


Parameter 

GAL 

DFW 

JSN 

Design  wavelength  (A.), 

2.7 

4.0 

5.2 

m  (ft)  d 

(9) 

(13) 

07) 

Acceptable  amplitude,  (A  ), 

9.1 

15.2 

22.9 

mm  (ft) 

(0.03) 

(0.05) 

(0.075) 

Soil  amplitude,  A 

30.5 

41.5 

45.7 

mm  (ft)  e 

(0.100) 

(0.136) 

(0.150) 

Aw/ 2a  =  Aa/A~ 

0.30 

0.37 

0.50 

B>  (Figure  8-4) 

3.6 

3.9 

4.5 

Equivalent  depth  (D), 

384 

699 

643 

mm  (in) 

05.1) 

(27.5) 

(25.3) 

Existing  pavement  equivalent 

400 

820 

518* 

depth  mm  (in) 

05.7) 

(32.3) 

(20.4)* 

600** 

(23.6)** 

*OriginaI  pavement. 

**Incl tides  recent  overlay. 


plotted  in  Figure  8-4,  values  of  BA  were  obtained.  Now  using  the  k  value  and 
Xd  from  Table  8-1  along  with  E,  B  and  then  D  may  be  obtained.  The  results  are 
shown  in  Table  8-2  along  with  the  estimated  equivalent  depth  of  the  existing 
pavement  at  the  sites  studied. 

Table  8-2  shows  the  equivalent  depth  required  by  the  expansive  soils 
design  and,  for  comparison,  the  existing  depth  at  each  site.  At  the  GAL  site, 
the  existing  pavement  equivalent  depth  slightly  exceeds  that  required.  The 
original  pavement  at  the  site  was  structurally  inadequate,  and  an  overlay  was 
placed  in  1972  to  upgrade  the  pavement.  In  June  1980,  another  overlay  was 
constructed  to  upgrade  the  pavement.  Distress  in  the  pavement  was  attributed 
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to  edge  effects,  utility  locations  improperly  backfilled,  and  isolated  areas 
of  water  infiltration.  The  design  procedure  provides  a  means  of  evaluating 
the  pavement-soil  interaction  as  the  soil  responds  to  pavement  construction. 

It  is  assumed  that  the  surface  does  not  leak  (e.g.,  through  joints  or  cracks) 
and  that  the  edges  are  isolated  from  wetting  effects.  From  evaluation  of 
the  design  and  performance  of  the  GAL  pavement,  it  is  concluded  that  these 
requirements  were  not  adequately  fulfilled. 

Values  in  Table  8-2  for  the  DFW  site  indicate  that  the  existing  pavement 
is  thicker  than  the  requirement  for  expansive  soils  design.  Pavement  perform¬ 
ance  at  DFW  has  been  excellent.  In  this  particular  case,  the  final  design 
resulted  from  extensive  on-site  study  and  monitoring  of  test  sections.  These 
results  indicate  that  the  expansive  soils  design  may  not  be  as  demanding  as 
the  high  volume,  high  gross  weight  traffic  analysis  used  for  the  DFW  pavement 
design. 

Thickness  requirements  for  the  JSN  site  exceed  the  equivalent  depth  of 
the  original  pavement  as  well  as  the  total  depth  after  a  recent  overlay  on  the 
taxiway  studied.  This  discrepancy  indicates  that  the  pavement  originally  con¬ 
structed  was  inadequate  for  the  expansive  soil  at  the  JSN  site.  Pavement  per¬ 
formance  has  been  unsatisfactory  at  JSN,  and  serious  problems  have  been  appar¬ 
ent.  Existing  pavements  at  the  gates  were  originally  thicker  than  the  taxiwav 
and  runway  structures.  For  example,  the  equivalent  thickness  for  JSN  11  and 
12  (Table  7-5)  was  638  mm  (25.1  in)  which  is  about  equal  to  the  required  thick¬ 
ness  shown  in  Table  8-2.  One  would  anticipate  distinctly  different  pavement 
performance  from  an  adequate  pavement  (JSN  11,  12)  when  compared  to  an  inade¬ 
quate  pavement  (JSN  5,  6).  Study  of  the  elevation  profiles  in  Appendix  A 
reveals  that  the  inadequate  pavement  shows  a  continuous  pattern  of  mounds  and 
depressions.  In  contrast,  the  distortions  are  centered  around  a  single  joint 
on  JSN  11  and  12,  indicating  surface  moisture  infiltration  through  the  joint. 
Thus,  JSN  5  and  6  evidence  general  structural  distress  because  of  insufficient 
pavement  stiffness  whereas  JSN  11  and  12  show  isolated  distress  caused  by 
infiltration  via  surface  discontinuities.  Proper  design  must  insure  that  pro¬ 
visions  are  made  against  both  forms  of  distress. 
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Alternative  designs  were  considered  for  adequate  structural  stiffness  for 
the  soil  conditions  at  JSN.  Two  alternatives  were  computed  following  the  pro¬ 
cedures  outlined  previously  and  are  shown  in  Table  8-3.  Alternative  1  incor¬ 
porated  additional  base  and  stabilized  soil  in  the  design.  The  required  thick 
ness  was  reduced  by  reducing  the  active  zone  depth  and,  thereby,  the  suction 
change.  Thus,  the  required  depth  was  599  mm  instead  of  643  mm  as  shown  in 
Table  8-2.  Alternative  2  increased  only  the  Portland  cement  concrete  thick¬ 
ness.  Both  alternatives  produced  equivalent  depths  which  exceeded  those 
required,  resulting  in  structurally  adequate  pavements. 

TABLE  8-3.  ALTERNATIVE  PAVEMENT  DESIGNS  FOR  JSN 


Material 

Alternative  1 
mm  (in) 

Alternative  2 
mm  (in) 

Portland  cement  concrete 

241  (  9.5) 

305  (12.0) 

Cement-treated  base 

305  (12.0) 

152  (  6.0) 

Lime-stabilized  soil 

254  (10.0) 

— 

Total  equivalent  depth 
of  AC 

615  (24.2) 

663  (26.1) 

Required  equivalent  depth 

599  (23.6) 

643  (25.3) 

LIMITATIONS  AND  RECOMMENDATIONS 

The  method  presented  above  provides  a  rational  procedure  for  selecting  a 
pavement  thickness  sufficient  to  mitigate  the  effects  of  swelling  soil  behav¬ 
ior.  Several  aspects  of  this  method  are  important  and  cannot  be  ignored  in 
applications  to  real  problems. 

Localized  areas  of  soil  saturation  cannot  be  overcome  by  the  design 
thickness.  If  the  pavement  leaks  or  seepage  beneath  the  pavement  occurs  (e.g. 
along  a  utility  installation),  additional  swelling  and  pavement  distress  may 
be  expected.  Therefore,  maintenance  of  joints  and  cracks  and  the  protection 
of  shoulders  from  water  are  essential  design  features.  Drainage  must  be 
adequate  and  properly  maintained. 
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The  presence  of  a  near-surface  water  table  or  a  perched  water  table  com¬ 
plicates  the  evaluation  of  suction  conditions  in  the  soil  profile.  The  pres¬ 
ent  state-of-the-art  is  inadequate  for  establishing  the  initial  and  final 
suction  profiles.  Additional  work  is  needed  to  better  define  equilibrium 
conditions. 


An  additional  complication  is  the  presence  of  a  significant  osmotic  suc¬ 
tion  component.  This  portion  of  the  total  suction  may  vary  dramatically  while 
the  water  content  remains  fixed.  It  also  may  change  very  slowly  with  time. 

A  rapid  method  is  needed  to  evaluate  the  osmotic  suction  component  and  deter¬ 
mine  its  change  potential. 


Follow-on  work  is  needed  to  compile  a  data  base  for  the  comparison  of 
predicted  pavement  requirements  with  performance  data  of  existing  structures. 
This  type  of  validation  is  a  necessary  step  in  developing  the  concepts  pre¬ 


sented  into  useful  procedures.  The  results  shown  in  Figure  8-4  should  also  be 
converted  into  a  design  chart  format  for  various  values  of  (3,  X,  and  k. 


At  this  time  it  is  recommended  that  values  of  Aw/2 a  be  between  0.1  and 
0.8.  If  Aw/2a  is  larger  than  0.8,  special  procedures  for  expansive  soils 
are  probably  not  required.  If  Aw/2a  is  less  than  0.1,  the  pavement  will  prob¬ 
ably  bridge  or  lose  contact  with  the  soil.  This  situation  may  result  in 
serious  pavement  damage  and  should  be  avoided  in  design. 
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SYMBOLS  AND  ABBREVIATIONS 


A 

AC 

Ac 

C 

CEC 

CEAc 

COLE 

CV 

D 

DF 

DFW 

DP 

E 

FP 

GAL 

AH 

HC 

JSN 

L 

LL 

AL/Ld 

AL/LZ 

NAS 

NMERI 

PCC 

PE 

PI 

R 

R2 

RC 

S 

TAXI 


Weighted  amplitude  of  frequency  components 
Asphalt  concrete 
Clay  activity 
Percent  clay 

Cation  exchange  capacity 
Cation  exchange  activity 
Coefficient  of  linear  extensibility 
Coefficient  of  variation 

Equivalent  depth  of  pavement  in  terms  of  asphalt  concrete 

Deficit,  Thornthwaite  method 

Dallas/Fort  Worth  airport  field  experiment  site 

Depletion,  Thornthwaite  method 

Modulus  of  elasticity 

Filter  paper 

Gallup  airport  field  experiment  site 
Elevation  change 

Water  holding  capacity,  Thornthwaite  method 
Jackson  airport  field  experiment  site 
Overall  sampling  length 
Liquid  limit 

Linear  strain  relative  to  dry  dimensions 
Linear  strain  in  the  vertical  direction 
Dallas  Naval  Air  Station  survey  site 
New  Mexico  Engineering  Research  Institute 
Portland  cement  concrete 

Potential  evapotranspiration,  Thornthwaite  method 

Plasticity  index 

Rainfall,  Thornthwaite  method 

Coefficient  of  determination 

Recharge,  Thornthwaite  method 

Degree  of  saturation;  surplus,  Thornthwaite  method 
Computer  code  used  to  calculate  aircraft  accelerations 
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TCP  Thermocouple  psychrometer 

TI  Thornthwaite  moisture  index 

AV/V  Volumetric  strain 

AV/Vd  Volume  change  with  respect  to  dry  volume 

AV/V^  Change  of  volume  with  respect  to  initial  volume 

WES  U.  S.  Army  Waterways  Experiment  Station 

a^  Half-amplitude  of  the  nth  frequency  component 

e  Void  ratio 

emf  electromotive  force 

f  Lateral  restraint  factor 

f  Spatial  frequency  of  the  nth  component 

h  Suction 

Ah  Average  suction  in  the  top  1.5  m 

k  Subgrade  modulus 

m  Moisture  characteristic 

h^,hi  Final  and  initial  suction 

pF  Logarithmic  unit  for  suction 

w  Water  content  of  soil 

Aw/ 2a  Ratio  of  pavement  amplitude  to  uncovered  soil  amplitude 
p\  Nondimensional  stiffness  characteristic 

AA  Slope  of  the  A  versus  Ah  relationship  (Figure  7-1) 

Yd  Soil  sample  dry  unit  weight  after  oven  drying 

Y ^  Suction  compression  index 

>  Soil  sample  dry  unit  weight  at  2.53  pF 

i  Total  stress  compression  index 

A  Wavelength 

\d  Characteristic  wavelength  for  design 

A  Wavelength  of  nth  component 

\0  Characteristic  wavelength  of  uncovered  soil 
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APPENDIX  A 
ELEVATION  PROFILES 


6467.5 


gure  A-l.  Elevations  along  profiles  at  GAL. 


Figure  A-2.  Elevation  data  along  Gallup  runway 
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Horizontal  distance,  ft 


Horizontal  distance. 


Horizontal  distance,  ft 


Figure  A-5.  Elevation  profile  on  DFW  access  road 


Horizontal  distance. 


igure  A-8.  Actual  profile  of  Dallas  Naval  Air  Station. 


Horizontal  distance,  ft 


Figure  A-9.  Actual  profile,  Dallas  Naval  Air  Station 
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Figure  A-10.  Elevations  along  OSN  taxiway. 
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Figure  A-ll.  Elevation  profile  along  apron 
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Horizontal  distance,  ft 


Horizontal  distance 


Horizontal  distance,  ft 


Actual  profile,  Jackson  airport  General  Aviation  Apron. 


APPENDIX  B 


FREQUENCY  DOMAIN  DATA 

The  following  data  are  given  for  the  regression  coefficients  c, 
regression  equation 


where 

a  =  half  amplitude  of  frequency  components 
n 

f  =  frequency 

The  tabulated  coefficients  were  calculated  by  least  squares  for 
frequency  domain  data  generated  by  the  Fast  Fourier  Transform  used  ii 
study.  The  coefficients  were  determined  for  data  a ^  in  feet  and  f  i 
per  foot. 


Prof i le 

c  x  10'4 

b 

GAL  2  AVG 

1.240 

-1.067 

GAL  3  AVG 

0.305 

-1.368 

GAL  4 

1.051 

-1.337 

GAL  5 

0.834 

- 1 . 285 

GAL  6 

0.242 

-1.476 

DFW  2  AVG 

2.162 

-0.891 

QFW  3  AVG 

2.378 

-0.778 

DFW  5  AVG 

2.446 

-1.092 

DFW  6  AVG 

1.339 

-1.054 

DFW  7 

0.254 

-1.646 

DFW  8 

0.424 

-1.409 

DFW  9 

0.308 

-1.684 

DFW  10 

0.426 

-1.586 

DFW  11 

0.353 

-1.647 

DFW  12 

0.574 

-1.447 

DFW  13 

0.652 

-1 .516 

DFW  14 

0.806 

-1.379 

NAS  1 

0.120 

-2.007 

NAS  2 

0.242 

-1.760 

NAS  3 

0.180 

-0.917 

NAS  4 

1.822 

-1.061 

JSN  2  AVG 

1.926 

-0. 594 

JSN  3  AVG 

1.299 

-0.864 

JSN  4 

0.154 

- 1 . 882 

JSN  5 

0.450 

-1.585 

JSN  6 

1.340 

-1.087 

JSN  7 

0.243 

-1  662 

b  in  the 


the 

this 

cycles 


Profile 

c  x  10“ 

b 

JSN  8 

0.220 

-1.724 

JSN  9 

0.201 

-1.593 

JSN  10 

0.297 

-1.712 

JSN  11 

0.013 

-1.872 

JSN  12 

0.191 

-1.942 

JSN  15 

0.120 

-2.027 

JSN  16 

0.173 

-1.860 
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APPENDIX  C 

CALCULATED  PROFILE  CHARACTERISTICS 
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.29964951 E-06 


APA8 

•33430693E-01 
.?9P4155eE-ei 
.17608648E-01 
. 7 36 9 243 3E -02 
.2012 131 2 E-01 
. ?1 261430  F-01 
.14132159  E-n 
.21568S46F-"1 
.  14542539E-'>1 
.21569821 E-81 
. 1 3676965E-01 
.  22064^49  E-<11 


178 


LINS  -ATS  SJMAN 


LFV 

i 

A?0  1979 

.  155 12261 F +21 0 

.42>5021F“.,2 

r?*1 

l 

F  KAY  197S 

.2  ?4376F  i»v)M 

.  02527  c.  3  l’-'I 

LF^ 

l 

JUI  1979 

.246709496  >3 

.55 1021::  IF -O? 

LFW 

l 

SSF  li>79 

.26095^299+^7 

. 539567  4? •> -02 

I  FV 

i 

4  NCV  lo79 

.1995437,  2 +00 

,  2 r 6B0285.;  ~ :?? 

if* 

2 

36  HP.i  19  7 

.55r  effaets-n 

.1  8566,-237  -03 

I FW 

2 

3  KAY  19 

.586667? 

.2  64761 23 1 -03 

r.Fv. 

2 

2?  JUL  IF? 

.  n  191827  .T -01 

.1  1563-73 

IF* 

2 

27  3  IF  19  7 

.534 20699} -01 

.2111152?; -03 

I'FV 

2 

14  NOV  19 

.53150761 P-*)  1 

. 1 7662 3o j  -33 

I F* 

3 

At/  lave 

.2l9V'0ii  F— 01 

.91391462  >04 

r  j 

n 

3  i“  A Y  1979 

.  3 ,;7 24622 '  -  1 

.65625  f -04 

w  *  & 

•7 

u 

JUI  1979 

.  7  144052'3 1-  3 1 

.  82720963  >  > 

!.?;• 

St?  19-79 

.3112716*,:;-  »i 

.  va  500675  -34 

D  7  # 

T 

4  NCV  197 J 

.  V,3c257.5.-31 

.6  -'43  b  37/  "-04 

IFV 

5 

A  ?  .  1979 

.6533707  «■  3- ;1 

.  4 5 1 8 -x cj 4 3  >>;3 

?  Fa 

3 

*  iOY  1979 

.6668525-  *>01 

.3.»82lf-> 

E  J 

c. 

JUI  19? 9 

.53915729: -01 

..3b  ?2b937  23 

I  > 

5 

37.?  1379 

.665  5C  56 '>-01 

.477':  17  24:-  33 

r?« 

5 

4  0  V  1 3  7  9 

>  393P462:,-01 

.  -i  1  4 5  4  72  ;  -0 

IFV. 

o 

'•56  API  lb 7 

6J292:''-0l 

. 214625750-05 

IFir' 

7  '  MAY  19 

.0?76l62?E~Ol 

.  21 V46593  •  -  -'I. 

IF  A 

l' 

27  JUL  19? 

.4315502LP-01 

.1*  68V? 4c,  - 

DF* 

k3 

"7  19? 

.44  5  541?  It. -01 

.35  359  263  .-  ,3 

r-F» 

6 

14  N  (.•  V  19 

.271767501;  -41 

. 26221 1 02 v-73 

r.~\\ 

1.3 

AP  ,  1» 

.195  30971  >00 

.  24  9.63232--?? 

ru 

13 

JL'I  lo?o 

>?2cCol.,:r'+0  Y 

.  6  >4.0445  c,-02 

.  •34336839 E— 31 

.4916o 6381+60 

.  *6177142;  -  "l 

.?43i» 1-  'i 

.  123092062+  :•/• 

.53862264 » 

.7  '96fi?592-  '3 

.1  '6o7c67V- 

.3 4C575127+0Z 

.2.3924  773  8+  30 

. 19037?; 9. 

.74  3  ?  1 

.3  569S161E+0U 

.  223564592+  ,'0 

.  1  ,1917.1  .  -  P 

.  '*■  6  0 ' 1 57  4  f.  -  > 

.24146676  3  +  0? 

.2857i.  1687+02 

.GO  7  81  >  >7.3 

.6  06i.-462F-ri 

.121423063-61 

.  1-'*P37  36?  v+«.  ? 

.2432*' 5  22  •  -25 

.  2 1 7  J  1  5 1  ?  r  -  .5 1 

.14190O61L-/1 

.146742541 +J0 

.3^ 0  034  52  >05 

.•'5.-3  »•  -  ■: 

.  13602  1 67  .>'51 

.  15  080 1 7  OF  +!•>’ 

.''79>it?~.  -~5 

.201117ft 

.1 39(07771  T-fll 

.  12 422966 F +00 

.2912:  496F-  '6 

.213  J  9922- "'1 

.1343  >33271-01 

.11475228^+00 

.24  104  '•  5 

.7  ,t:  .07861-  ’1 

.  3  9r?13722F-02 

.255974275+00 

.  >2*221  wF-09 

. .?7t:55f.7F- ••' ' 

.42400665F-O2 

.14267556?-  7 

. ;  v  07  4 '’6>  > 

.11'  91  4  72  >31 

.44411954E-0V 

.80^551  -  -”•'  +  0? 

.  .I  :’i:  1. -> 

. ;  : .  o9.  }-  'i 

.  3949?  39P>-'2 

.  19?  '2597'  +92 

.3  r'8^7;-  vt.  *.  -/r. 

.130 "22 5 1  T-  5.  t 

.43381240F-0? 

.  15  J9996'1  -  t 

.  '  851  15  21-  ’6 

•  11  1  >84F-’l 

.11  141175  c-Ji 

.  15L  85  2  '4  l  +t'9 

•  '  2 1 '  0  »  —  1 0 

.  1  127  4  >  51 

.192744523- Jl 

. 17  +  8785 7 F  +  5  0 

.?40 

..’y!5 :  ■•>  ol  F- 

.5 ?65964  4F-02 

.16658149*+  0 

.  .244°?'’»'r 

~t f.  r_  •«  p  •  -  o  r  -  *•  ’ 

. 11757  873F-J1 

.  1  *’491624'.' +C0 

.5‘7>  A  •-  3 

.*:!  1  9464  t-  '1 

.901605179-0? 

.176417309+  0 

.''7371,  7L-06 

.;‘-c2:52i”-*',i 

.266761  sJE-ti  ? 

.3223624  1 F  +•  0 

.  •;726‘T.  52  -  ’-r 

.  3.2r.  fpf  y-  >1 

.32198142 F-02 

.324  83046  F  +  x.0 

.'  988119*5  2-  .6 

.824-  1767 

.  A  0862  894  E_,,2 

.26°!42229>80 

.  ?:'4  075O  JF-  Vj 

.2  '  " .V;  F-‘i 

.3233576 AF-0? 

.3592V 133V.  g 

.  1 1918.727  t  -Or. 

.:  .  ;?f  9  -"F —  1 

.294270525-02 

.  05  54636  0> 'HI 

.7716097.  --06 

.  ?  8 1 7 ' 1 7  r-  *■  1 

.16 174170F+20 

. 1 2 5 8 1 5 4  9  r  -90 

.4  136  Jl  131  -77 

.*,  •  , 

.c  454  643  V»P4 

.247- 5432? +  8 

. '  '*7  •  9: r  r-?r 

.71  '1 1 71  '  0-  1 

1  79 


iu  hi) 


1 

7 

5 


LINE 

LATE 

S  IM AN 

M 

NAS  1 

JUN 

1980 

.34 } 27 646 E -01 

. 15245746E-02 

NAS  2 

JUN 

1980 

.71972982E-01 

.10869810F-02 

NAS  3 

JUN 

1980 

.40478820E-01 

.13147522E-03 

NAS  4 

JUN 

1980 

.498513561-01 

.20230519E-03 

P 

K 

E 

A  BAH 

.12850250E-01 

. 17511117E+00 

. 19591164E-24 

.711568341-21 

.91310421E-02 

.  30807 089E+00 

.99252692E-05 

. 52845070E-01 

.61177849E-02 

.  11275971E+00 

. 8043371 IE-05 

.20119757E-01 

.90007798E-02 

.145S2777E+00 

.18209045E-05 

.231426141-71 

LINE 

EAT  i. 

SJhAN 

M 

D  J*  1 

JUN 

19Su 

59572998E  +k;0 

523948  i2e.-*l 

LI  a  2 

JUN 

1S36 

.594339596 -3 1 

.21967951: -23 

D Ft  3 

JUN 

13=0 

T 04 05312 F -01 

.74531073  -.'-'H 

r  n  5 

JUN 

1960 

.7  ,c30134£-71 

.5=2812416-03 

rr*  e 

J’JN 

1980 

.384 57728"-9l 

.21 4063??- -03 

JSN  1 

JUN 

1980 

.531927951 +UU 

2‘'b?  06 12  "-bl 

JSN  2 

J  JN 

lb  60 

154  55682E-91 

.1?  1892730-^4 

JSN  3 

JUN 

I960 

.  23783456  E-J 1 

.55677478  —  '4 

P 

tC 

4 

ALA  i 

.126840357+01 

.363893077+0* 

.6645769  lr—*l 

.  V*5*51 9  E+-'* 

.1 1396294 E-?l 

.  1261  '’38SE +09 

.250307651-05 

.25  871189  F--  1. 

.33784 9B5E-02 

.21 461906P+  ,:0 

.  *;P9 65 8 5^1  -  i6 

. 127632:91-01 

.1553735SE-01 

.21823276E  +  77/ 

.781237637-05 

.  6644.3'x31 T?- ‘*1 

.336237218-02 

.  3033753=v+00 

.71976322'. -06 

.  23  ='61  841  y-17! 

.76'- 10839  E+OU 

.28003 7B6F+00 

.  78085150 r  -  *1 

.5=066552F+33 

.18Q31145V-  2 

.71309260E-61 

.22765485 \-i?7 

.  ':5221^27  '-67 

.2^64439=8-92 

.  17’5’6o601T,J-'1*« 

.  1147363C  ,'-'>6 

. 1 14  c? "4  7  F-,',l 

:  lip 

£*  l  .  ry 

«  J  \ 

••  T  L‘ 

.1  * 

1  J  V' 

.6.1  -.1 

■  '■*•—*•11  __ 

T  • 

i-  - 

JU'1 

U';:  I 

.■  V  <75  ’  >■■-.  1 

T-  /  L 

t  ,  ;  • 

1  b>3  0 

.;^'’i::'.  -v 

.J  -V  \  '•  -  ^ 

rit  l 

T  ■  |  • 

i  « 

1  Jr 

-  •'  1 

•*  7  1  __  ) 

•  4.  4 

r  /  ■  ii 

t]  •  ?  \| 

3  i  r 

i.  • 

.  1  .  J.'  ,  >  '  * 

•  '  ’"l  •  l  ♦  •  ’• :  ~ . 

I  :•  ■.  1 2 

O  j  t  . 

1  ir  *  * 

.  bj  -:i 

+  >  r.  ■  l.  _ 

~  7  *  r 

C 

Jl  ’• 

.  ?r.  •  - .  1 

•  •  ~  7 1  '  •  - 

r  "  *;  2,  *1 

J  UN 

4  )  ' 

^  •*  J 

.*?:  -  -  . 1 

•  *  X  j 

.=  67  7*  773  r-  2 

.1  -  Yt  ■ 

:  ~  .—  4-i50"  *■  .  J 

*  r  * 

r.  .  1  ~  r-  *»■ 

•  ■  1  A  ■  '  A.  .  * 

.  :  •••’711 '  - 

.  j  17:  • 

. -  . 

1  •  r  "*  .>  t,  “  r*.  • 

.  1 ’  1  r  ’  - 

.  '  •;  ‘  1 

.  .  •  _  ■ 

.  1  j :  ■  J  5 r  ’  5j  5i "—  .1 

.  1 

*•  ’•  ^  '  ■ 

:  ■  -  •• 

.'  i .  244  !u  7—  0 

.  ^  .j  1  ’  t  '  32  •  *■  •. 

.  -  •  • 

.  1  ’3  : .  : 

.  '  j  .  •  ' 

•  * 

■  -  • 

« 1 1!*  b  1 1  t  o  ~~  1  k 

.■.'574  i:  * 

- 

•  1  ’  - 
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LINE 

DATE 

SUM  AN 

M 

JSN 

4 

06  DEC  197 

.  19846  089E  +00 

.138208537-01 

JSN 

1 

APR  1979 

.142467971+00 

. 189268691-02 

JSN 

2 

19  APR  197 

.17875657E-01 

.203172006-04 

JSN 

3 

APR  1979 

.  22407  545E--01 

.375599007-04 

JSN 

1 

MAY  1979 

.140346255+00 

.20200758E-02 

JSN 

2 

31  MAY  197 

.  169 10389 E -01 

.25115522E-04 

JSN 

3 

MAY  1979 

.234295805-01 

.432588407-04 

JSN 

1 

JUL  1979 

.  149  55439E  +00 

.224028537-0? 

JSN 

2 

27  JUL  197 

.209544421-01 

.552200087-04 

JSN 

3 

JUL  1979 

.18621460E-01 

.23126644E-04 

JSN 

1 

8  SEP  1979 

.14220493E+00 

.205241 16«-e? 

JSN 

2 

28  SEP  19 

.16327225E-01 

.1 64667 53E-04 

JSN 

3 

8  SEP  1979 

. 31453 024E-01 

.56597235E-04 

JSN 

1 

DEC  1979 

.1511906SE+00 

.230522877-02 

JSN 

2 

06  DEC  197 

.13608983E-01 

.84182689E-05 

JSN 

3 

DEC  1979 

.23452169E-01 

.521316067-04 

p 

.378869483-01 
,3436761 1F-01 
.24015198E-02 
. 27381 297 E-02 
.73056398E-01 
.  1 8245302E-02 
.23976148E-02 
.75 8301671-01 
.18818905E-02 
.176908503-32 
.76184055E-01 
.190835117-02 
.534725583-02 
.80085619E-01 
.16433379E-02 
.19108419E-02 


K 

.25538817E+00 
.24263042E+00 
.17413622E+00 
.  11784093E+00 
.24964165E+00 
.27425319E+00 
•11834539E+00 
.  265412655+00 
.36550579E+00 
. 1 178081PE+00 
.28070799E+00 
.15161569E+00 
.763376657-01 
.27508995E+00 
.49211165E-01 
.172093693 +00 


v 

.  36479221 E +00 
.22433833E-01 
.846014247-02 
. 13717356 7-01 
. 2765091 2  v-01 
. 13765479E-01 
.  18042448!''.  "-01 
.295551501-01 
. 293479785-01 
. 130726587-01 
.269401737-01 
.962878583-02 
. 10564 351E-01 
.28784552E-/.1 
.512266453-02 
.272820103-01 


ABAS 

.175169893  '  ‘?'f 
.  59328451  E-  9. 1 
.413074737-0:' 
.  88698354 E-U? 
.666  360601-21 
.  48656613E-P2 
,572735347-0? 
.705115367-01 
.867495757-0? 
.68497317  2-02 
.66910c97E-ri 
.411540557-02 
.1211 /036E-2 1 
.6.^951  ??3  7-01 
.22  ■■■74  <*}  57-0? 
.113368907-21 


1 


LINE 

DATE 

S  JNAit 

M 

JSN  5 

JUN 

1980 

.  13507027E +00 

.260l6571R-'»2 

JSN  6 

JUN 

1930 

.73272325E-01 

.103446136-02 

JSN  7 

JEN 

1980 

.96610486E-01 

. 27511351 £-02 

JSN  8 

JUN 

1980 

.  12306 096E +30 

.7  264  0783.5 -(’2 

JSN  9 

JUN 

1930 

.  5635345 IE-01 

.949215437-33 

JSN  13 

JUN 

1980 

.75531455E-01 

. 81397933 K— S3 

JSN  11 

JUN 

1960 

.  11585075E  +00 

.33300288E-f’2 

JSN  12 

JUN 

1980 

.1P436389E+03 

.99341126^-02 

JSN  15 

JUN 

198J 

. 199 1407  4E  +  30 

.11501396 '-'-01 

JSN  16 

JUN 

1960 

. 161B9273E+ 00 

.£14288533-02 

p 

.1163138315-31 
.74994418 15-02 
.82153653 ?-02 
.1124332615-01 
.3139678715-32 
.9617850?  f-02 
.1 207382515-01 
.26555911F-21 
.34487655E-01 
.212  J6097E-U1 


K 

.222203622+06 
.248242772+60 
.3dlv51221E+02 
.47754066E+<6 
.-.1  j269?7r+?0 
.22087c74E+00 
.27089651-5+ 6  M 
.25943S71E  +  k,0 
.22932S22E+JW 
.3'--78631wS+00 


r 

.3326217t'E-',*4 
.7757L861E-05 
.225019S1E-?! 
.81708721E-«4 
.298e?3 147-06 
.77998753E-25 
.4  18298761.-34 
.26301274:5-^2 
.  3966734  lf.-'JO 
.  192^7? 74  2 


A  9  A  .<• 

.851 16264 E- 31 
.  CJ617  Jf;l  E-21 
.'’7188462  tl 
.  l-v  4  ^745 7 +0i’ 
.'/I424’~t8b¥-'>1 
. 637354 12E-01 
.  9657’585  i’.-/,l 
.  lbt  3922 P7  +  0'-‘ 
.  17r-:66  4l3y  +  -)i' 
.  l42ll?16!5  +  ,,0 
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FORMAT  FOR  PRINTOUT  OF  AIRCRAFT  DATA 


AIRCRAFT 

MODEL  GROSS 

WEIGHT 

W 

A 

B 

PSARM 

TAILRM 

SPEED 

THRUST 

TAKOFF 

CL 

AREA 

CD 

WM 

WN 

SXM 

SXN 

AHN 

AAN 

AHM 

AAM 

PAON 

PAOM 

VON 

VOM 

SLM 

SLN 

TSM 

TSN 

STROKE 

PIN  DIAMETER 

NOSE 


STROKE  PIN  DIAMETER 
MAIN 


SIMAIN  SINOSE  SICG  SITAIL  SIPS 


GEN  MASS  OMEGA 


OAM 


OAN 


184 


**********************************************************  TAXI 


READ  AND  PRINT  INPUT  CATA  TAXI 
W=V  EH  I CLE  WEIGHT  AT  CG  (POUNDS)  TAXI 
A=DISTANCE  MAIN  GEAR  TO  CG  (INCHES)  TAXI 
B=EI STANCE  NOSE  GEAR  TC  CG  (INCHES)  TAXI 
MMI=MASS  MOMENT  OF  INERTIA  (LB  IN  SEC  SQ)  TAXI 
PLANE=AI RPLA  NE  BEING  SIMULATED  AND  GROSS  WEIGHT  TAXI 
PSARM=DI STANCE  OF  PILOT  STATION  TO  CG  TAXI 
TAILRM=EISTANCE  CF  TAIL  STATION  TO  CG  TAXI 
TAKOFF=TAKE-OFF  SPEED  (FEET/SEC)  TAXI 
SPEED= INITIAL  VEL  OF  AIRPLANE  TAXI 
THRUST=TOTAL  AIRPLANE  THRUST  TAXI 


CL=LIFT  COEFF.  TAXI 

AREA  =  W I N G  AREA  TAXI 

CD=DRAG  COEFF.  TAXI 

WM*WEI GHT  OF  MAIN  GEAR  (EACH)  TAXI 

WN=WEIGHT  OF  NOSE  GEAR  TAXI 

SXM=NUMBER  OF  MAIN  GEAR  STRUSS  TAXI 

SXN=NUMBIR  OF  NOSE  GEAR  STRUTS  TAXI 

AHN=HYDRAULI C  PISTON  AREA  NOSE  SO  INCHES  TAXI 

AAN  PNEUMATIC  POO  TAXI 

AHN  HYDRAULIC  PISTON  AREA  NOSE  SQ  INCHES  TAXI 

AAN  PNEUMATIC  PISTON  AREA  NOSE  SQ  INCHES  TAXI 

AHM  HYDRAULIC  PISTON  MAIN  SQ  INCHIS  TAXI 

AAM  PNEUMATIC  PISTON  AREA  MAIN  SQ  INCHES  TAXI 

PAON  NOSE  STRUT  PRELOAD  PRESSURE  FSI  TAXI 

PAOM  MAIN  STRUT  PRELOAD  PRESSURE  PSI  TAXI 

VON  NOSE  STRUT  INITIAL  VOLUME  CU .  IN.  TAXI 

VCM  MAIN  STRUT  IN  IT  AL  VOLUME  CU.  IN.  TAXI 

OAM  ORIFACE  AREA  MAIN  TAXI 

OAN  ORIFACE  AREA  NOSE  TAXI 

SLM=MA I N  GEAR  STRUT  LENGTH  UNLOADED  INCHES  TAXI 

DISTANCE  FROM  CL  CF  AXLE  TO  CG  LINS  TAXI 

SLN=NOSE  GEAR  STRUT  LENGTH  UNLOADED  INCHES  TAXI 

DISTANCE  FROM  CL  OF  AXLE  TO  CG  LINE  TAXI 

TSM  MAIN  TIRE  SPRING  CONSTANT  PER  STRUT  TAXI 

TSN  NOSE  TIRE  SPRING  CONSTANT  PER  STRUT  TAXI 

DX=T  IME  STEP  SIZE  TAXI 

READ  METERING  PIN  DESCRIPTION  STARTING  AT  ZERO  STROKE  TAXI 

NS  CN  =  OF  METERING  PIN  CHANGES  NOSE  GEAR  TAXI 

NS CM=  OF  METERING  PIN  CHANGES  MAIN  GEAR  TAXI 

NFM=NrJM BER  OF  FLEXIBLE  MODES  TAXI 

SIXXXX(I)=MODE  SHAFS  DEFLECTION  (NON  DIM.)  TAXI 

GM( I )=GFNERALIZEL  MASS  (POUNDS-SE C/IN)  TAXI 

OMEGA ( I )=MODAL  FREQUENCIES  (RAD/SEC)  TAXI 


ft**********************************************************  TAXI 


J 
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WIT?  ?Q 

I ID  TPI-JT? 

C  OMPOS  IT? 

391500.0 

56.0 

767 .0 
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APPENDIX  E 
CLOD  TEST  PROCEDURE 
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Clod  Test  Procedures  and  Data  Reduction 


General  Description 

1.  Clods  weighing  120  ±  20  grams  are  separated  from  undisturbed  samples  and 
placed  in  4  ounce  moisture  cans;  usually  seven  for  the  initial  test. 

2.  Three  are  wetted,  one  remains  at  natural  moisture  and  three  are  dried  to 
achieve  a  variety  of  moisture  conditions.  For  example: 

18,  20,  22,  24  (natural),  26,  28,  30  percent 

3.  When  the  desired  moisture  condition  is  achieved,  a  sensor  is  placed  in 
the  container  and  allowed  to  equilibrate.  The  sensor  may  be  filter  paper 
or  a  thermocouple  psychrometer.  Separate  procedures  for  each  are  used. 

4.  The  containers  with  the  sensor  and  the  sample  are  placed  in  an  insulated 
chest  to  equilibrate.  After  equilibration,  the  suction  sensor  is  removed 
and  the  procedures  for  the  remainder  of  the  clod  test  are  independent  of 
sensor  type. 

5.  Clods  are  weighed  (Wi),  followed  by  preparation  of  the  clod  for  bulk  den¬ 
sity  measurement.  A  wire,  tag,  and  possibly  a  hair  net  are  attached  to 
provide  a  means  of  handling  the  sample.  These  tare  items  may  be  added 
before  adjusting  the  moisture  content  if  desired  in  order  to  prevent 
difficulty  in  the  wetter  region.  Experience  has  shown  that  many  clods 
crumble  when  wet  if  handled  at  all.  Following  this  installation,  a 
weight  is  taken  (W2). 

6.  The  clod  is  then  coated  with  saran  resin.  Solutions  of  1:7  or  1:4,  saran 
methyl  ethyl  ketone  are  used.  The  1:4  is  used  only  for  very  porous 
materials.  Sequence  for  coating  is  as  follows: 

a.  For  1:4  Solution--dip  in  liquid,  dry  five  minutes;  dip  in  again,  dry 
fifty-five  minutes. 

b.  For  1:7  Solution--dip  in  liquid,  dry  five  minutes;  dip  again,  dry 
eight  minutes;  dip  again,  dry  fifty-five  minutes.  (This  is  pat¬ 
terned  after  the  SCS,  NSSL  COLE  procedure.) 

7.  Immediately  weigh  the  clod  in  air  and  water  at  the  end  of  the  drying 
period.  These  weights  are  W3  and  W*  respecti vely.  In  this  test  proce¬ 
dure,  the  buoyant  force  on  the  clod  is  measured  rather  than  the  submerged 
weight. 

8.  Clods  are  then  dried  in  air  until  they  reach  a  constant  weight.  A  weight 
in  air  and  water  is  determined  as  before.  These  are  W5  and  respec¬ 
tively. 

9.  Clods  are  placed  in  a  cold  oven  which  is  started  and  raised  to  105°C; 
they  are  dried  for  48  hours. 
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10.  The  clods  are  removed  and  cooled  until  easy  to  handle.  Weights  are  taken 
in  air  and  water  again  (W7  and  We).  The  weights  taken  in  water  (W4,  W6, 
We)  ate  the  buoyant  force  on  the  clod  rather  than  clod  weight  while  sub¬ 
merged. 

11.  The  saran  plastic  coating  used  loses  10-20%  by  weight  during  oven  drying. 
The  Soil  Conservation  Service  uses  15%  weight  loss  in  COLE  calculations. 
NMERI  tests  using  oven  dried  alumina  blocks  as  samples  indicated  the  aver¬ 
age  weight  losses  as  follows: 


BLOCK 

COATING 

WEIGHT  LOSS 

q  % 

A 

None  ' 

0.01 

0.01 

B 

7:1 

0.14 

16.7 

C 

4:1 

0.26 

16.7 

In  NMERI  calculations,  a  15%  weight  loss  is  assumed. 

12.  Sensor  measurements  involve  the  determination  of  water  content  of  a  fil¬ 
ter  paper  or  output  measurement  of  a  thermocouple  psychrometer.  In  both 
cases,  the  quantity  determined  is  compared  to  a  calibration  curve  for 
the  particular  sensor  involved.  The  calibration  is  used  to  convert  the 
readings  to  suction  in  pF  units. 

SUMMARY  OF  MEASUREMENTS 


A. 


B. 


Clod  Samples 

Wj  =  W  +  W  (Weight  of  wet  clod) 

5  W 

W2  =  Ws  +  Ww  +  T  (Wet  clod  plus  tare) 

W3  =  W  +  W  +  T  +  W  (Wet  clod,  tare, 
s  w  r 

W4  =  Buoyant  force  on  submerged  clod 
W5  =  Ws  +  (Ww)a  +  T  +  Wr  (Air  dried  to 
W $  =  Buoyant  force  on  submerged  clod 
T6  =  Water  temperature  at  which  W6  was 
W7  =  Ws  +  T  +  0.85  (Wr)  (Oven  dried) 

W8  =  Buoyant  force  on  submerged  clod 
T8  =  Water  temperature  at  which  W8  was 


saran  coating) 
constant  weight) 
determined 


determined 


Suction  (Filter  Paper) 


See  also  suction  measurement  procedures. 

Ida  =  Wf  +  +  Tc  (Weight  of  filter  paper,  water,  tare-cold) 

W i 2  =  W^  +  T^  (Weight  of  filter  paper,  tare-hot) 
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W x 3  =  Tc  (Cold  tare) 

W i 4  =  Th  (Hot  tare) 

Temperature  at  equilibrium  (measured  in  the  insulated  chest). 

C.  Suction  (Thermocouple  Psychrometer) 

See  also  suction  measurement  procedures. 

E  =  Output  in  microvolts 
Tmv  =  TemPerature  in  millivolts 
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